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FOREWORD 


The great challenge of our time is the exploration and 
conquest of outer space. To meet this challenge success- 
fully, NASA must have help from educational institutions 
and from industrial organizations. Such a cooperative en- 
terprise known as the Goddard Summer Workshop was con- 
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ducted at this Center recently. The results obtained from 
xperimentare summarized in this volume, 

jfe joint efforts of university faculty teams 

t^ni^oddard Staff members. Certain problems relating to 
sinlilation of space environments in the laboratory were 
investigated for the purpose of advancing our knowledge and 
understanding of the phenomena involved. Our greatest 
benefits are derived from this type of study partnership 
between academic researchers and the NASA family. The 
relationships are important to our future development and 
growth. 





PREFACE 


The Goddard Summer Workshop was born of necessity; in twelve weeks it grew 
to be a source of hope and promise for a better understanding of space -environment 
simulation. 


In a time of rapid change and unprecedented technology, it is vitally important 
to increase fundamental knowledge which applies to the reliability of spacecraft. 
This need for innovations in environmental testing was further emphasized last Jan- 
uary by the prospects of Goddard’s new space simulators. At that time, the NASA 
administrator, Dr. James E. Webb, happened to visit Goddard, and the discussion 
turned to the Buildings 7-10 complex with the diversity of equipments to be housed 
in them. The thought that by the end of 1962 this complex structure would consti- 
tute one of the most advanced facilities of this type in the world was gratifying to 
all who heard it expressed. 

Then there was talk about utilizing these elaborate facilities with maximum 
effectiveness. What are the various approaches toward achieving such effective- 
ness? Does the Test and Evaluation Division staff have the time and the in-house 
capabilities to explore novel ways and means which would provide the applications 
of new scientific knowledge that is becoming available in this field? Who will spec- 
ify the test procedures when the new installation is ready for use? The answers to 
all these questions indicated generally that the projected workload for the newly in- 
stalled facility would require a large addition to the present staff, just to meet the 
test schedules. Hence, it appeared very unlikely that there would be available in 
the Division enough qualified personnel with sufficient time to undertake these 
studies in the near future. Nevertheless, it was fully recognized that the need for 
such investigations must be met, and, that prompt action was necessary. 

Usually one turns first to a university for help in scientific research and de- 
velopment; that is why this exigency created the idea of a summer workshop. Im- 
mediately, attention was focussed upon the possibility of procuring for the summer 
a number of university professors and graduate students who are specialists in 
certain disciplines. The Washington area campuses as well as more distant uni- 
versities were canvassed for available and suitable people. As a result of this 
search, the academic groups listed in this publication have contributed whole- 
heartedly to the efforts described herein. 
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SUMMER WORKSHOP EXPERIMENT 



PART I 


SUMMER WORKSHOP EXPERIMENT 


INTRODUCTION 

This final report is prepared to serve a twofold purpose. It provides a means for 
estimating the value of a summer workshop experiment conducted at Goddard within the 
confines of a test and evaluation activity. It also records the contributions to several 
aspects of simulation testing made by a group of university scientists during a twelve- 
week period. The substantive and tangible results obtained by the academic teams are 
adaptable to immediate uses and, therefore, may be assessed promptly as they supply 
needed information. 

It is far more difficult to evaluate the long-range assets and benefits which this ex- 
periment could eventually bring to Goddard and NASA. To a considerable extent, the 
future progress of space science and technology depends. upon the potential returns from 
such ventures. This was hopefully expressed by Dr. Goett, the Director of Goddard, at 
the opening session of the Summer Workshop: 

"I consider the Summer Workshop a very, very interesting 
experiment, and I re-emphasize that the main product and 
value of this experiment would hardly appear during the 
twelve weeks that you will be working here. The more im- 
portant benefits should come later, from the familiarity with 
our problems which you people can pick up, and the degree 
to which you can focus your specialized talents on these 
problems. We are hopeful that you will be in a position to 
help us when we come to you. We hope also that as univer- 
sity teachers you will guide your brightest graduates to 
come here and help us do the job at NASA.” 


In a sense, this experimental adventure yielded a measure of satisfaction to all who 
were interested in its undertaking. The Summer Workshop made available some cur- 
rently needed information for immediate application to environmental testing. It awak- 
ened, in the farseeing investigator of space-environment simulation, a keen desire to 
know why certain procedures and techniques are followed in the evaluation of payloads, 
and what new ideas could be developed which would prove more useful to test engineers. 

The Summer Workshop may have raised more questions than it answered, but it also 
pointed the way to possible growth and progress for test and evaluation in space activi- 
ties. Moreover, the Summer Workshop helped to establish fruitful relationships between 
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academic researchers and practical users of science and technology related to space ex- 
ploration. Finally, it may be said that a worthwhile mission was accomplished if the aca- 
demic groups have motivated a few of Goddard’s engineers and technicians to learn more 
about their job and do it better than before. 

The more tangible results of the Summer Workshop, produced by the several aca- 
demic teams, are presented in the six subsequent parts of this report. Each of these 
parts start with a preliminary statement made at the opening session of the Summer 
Workshop, suggesting to the different university groups some desirable approaches that 
might be made toward the solution of a particular problem, and indicating the technical 
goal to be achieved. These preludes represent a few of Goddard T s ’’voices of experience” 
requesting help in certain areas of research and development. The results contained and 
implied in the six articles may be considered as a cooperative response from university 
faculty members who for three months contributed to the analytical and experimental in- 
vestigations of the tasks set before them. 


GENERAL REMARKS 

Success in modern space exploration demands that the flight vehicle be first care- 
fully tested on the ground. Experience has indicated that such an evaluation tends to en- 
sure the performance reliability of the spacecraft. In the words of Dr. Goett, Director 
of Goddard, 

”If we have learned anything at all during the past three years, 
we have learned that a tremendously important factor in the 
reliability of spacecraft is the thoroughness of the testing that 
is provided on the ground— and this building with all its equip- 
ment indicates the importance we place upon this particular 
part of the job.” 

Thus, ground testing has become a keystone in the reliability of spacecraft opera- 
tion, and the need for preflight evaluation is now well recognized. But the procedures 
and techniques for simulating space environments are by no means fully established, be- 
cause the knowledge of conditions out there is undergoing rapid changes, and new methods 
must constantly be devised to meet different situations. A vast amount of outer space 
information is swiftly unfolding. Consequently, the state-of-the-art in simulation testing 
is being transformed at an unprecedented rate. By constantly introducing these applied- 
science innovations into the test procedures, it is possible to do a more thorough job of 
ensuring spacecraft -performance reliability. 

With the shortage of qualified technical manpower in government activities, e.g., 
Goddard, it is almost impossible for the same group of test and evaluation engineers to 
carry out all the phases involved in advancing the technology of simulation testing. The 
background experience and unique skills possessed by the Goddard test group are vital to 
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the accomplishment of today 1 s routine testing, and must be devoted to it. But the addi- 
tional responsibilities of keeping pace with new developments in this field, and at the 
same time putting novel ideas into practical use for ground testing of spacecraft, are 
beyond the capabilities of this or any other engineering group similarly situated. Present- 
day simulation involves many disciplines; to integrate the diverse applications of science 
into a comprehensive set of test procedures requires not only ample time, but also new 
skills and versatility in several fields. Therefore, to take advantage of the new tools 
which are becoming available and to utilize them in this effort, it becomes necessary to 
supplement the available manpower with other professional talents so as to enable GSFC 
to carry out the job effectively. 

There are several ways of procuring the required technical assistance in this effort. 
It is possible to increase the staff by a few carefully selected specialists who could sup- 
ply the immediate and future needs of this activity. Such an approach, assuming the pro- 
curement of qualified people, would benefit the overall spacecraft testing, and would 
serve the best interest of the government in its long-range objective. But considerations 
of manpower ceilings, space limitations, and other temporary handicaps to Federal in- 
house research preclude the immediate activation of such a program. For prompt action, 
therefore, other approaches to this .problem should be explored, even though the pros- 
pects for much greater support of Federal in-house research and development are favor- 
ably regarded. In time, the situation will be remedied, because it is recognized that in- 
house research is likely to add to the government activity’s experience and prestige with 
possible breakthroughs. It builds competence and scientific know-how among its people; 
it enriches in-house talents and capabilities, and makes for growth and progress of the 
entire organization. 

Still, the job at hand requires immediate technical assistance, and it must be found 
outside the Center. A common procedure is to contract out government research and 
development work, and it should be feasible in this case. Perhaps a competent industrial 
organization would undertake to carry out, jointly with the testing staff on the premises, 
the different phases of the needed investigations. Another possibility is to invite indivi- 
dual experts or specialists in the technical areas under consideration to work with the 
permanent staff on the existing facilities to obtain the desired results. To be sure, the 
effectiveness of such arrangements depends upon careful planning for the use of test 
equipment and upon precise timing of personnel availabilities; otherwise, the primary 
objective of the investigation might be vitiated. Another factor which may influence the 
value of the overall study is the time spent in discussing the various aspects of the inves- 
tigation and the conditions under which the interchange of ideas must take place. It is 
rather difficult to assimilate a relationship between theory and experiment while the pres- 
sure is on to test a satellite which has an early target date to meet. 

In other words, the hiring of consultants and experts to help expand certain phases of 
engineering and applied science in the test program can be very effective. In general, 
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however, the busy permanent staff gains little from this temporary relationship, in terms 
of research experience and training for the future. Unless there is a realistic inter- 
change of pertinent information, and both groups participate in the effort on a mutually 
contributing basis, the residual benefits to the permanent staff and the government are 
evanescent in nature. From the long-range viewpoint, the hiring of consultants and ex- 
perts would prove a huge success if, besides achieving high standards in the current 
evaluations, the permanent staff gained in its basic understanding of the phenomena in- 
volved so as to be encouraged and stimulated to investigate on its own the environmental 
simulation of outer space. 


ORIGIN OF THE SUMMER WORKSHOP 

During the past few years fresh information has become available on space - 
environment simulation. This information is based upon new scientific knowledge of 
outer space conditions, which has been accumulating rapidly. In attempting to adapt 
some of these novel ideas to existing test procedures, the Test and Evaluation Division 
engineers have found that, before certain technological innovations can be put to practi- 
cal use, they must be subjected to considerable scientific investigation. 

For example, in testing the thermal equilibrium of a satellite whose coating does not 
absorb spectral energy uniformly, it is important that the solar-simulator spectrum 
should match that of the sun. Hence, before a rational evaluation of this satellite's ther- 
mal balance can be made, a careful investigation is necessary to find whether the artifi- 
cial source is similar to the spectral distribution of the sun. This effort entails an ex- 
tensive research in spectroscopy, a large departure from the normal function of test and 
evaluation engineers; yet, many such problems have arisen during the past three years. 
Solutions to these problems would provide more effective procedures for assessing the 
reliability of spacecraft, and would also help establish criteria and standards for opera- 
ting the many varieties of test facilities. 

In their desire to update the requirements for reliable spacecraft in orbit, the test 
and evaluation groups have attempted in spite of obstacles to resolve some of these prob- 
lems to meet modern testing needs. Several advanced development projects were under- 
taken in-house, and enough exploratory work was done to indicate the phases of the prob- 
lems which require attention. Although these efforts, interrupted by high-priority 
commitments, produced few substantive results, they provided by virtue of experience a 
number of particular topics which should be investigated promptly. These are: 

• Solar Simulation and Albedo Effect on Spacecraft 

• Heat Transfer and Cryogenic Pumping in Vacuum Technology 

• Corona and Discharge Effects Inside Spacecraft at Low External Pressure 

• Magnetic Simulation 

• Response of Vehicle Structures to Acoustic Fields 
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• Relating Damage and Reliability to Environmental Data 

Last January the Test and Evaluation Division concluded, that, under existing condi- 
tions, there is little probability of conducting in-house one or more of these studies to 
achieve definitive results; that the solutions to these problems are essential; and that 
technical assistance for these investigations should be procured from outside Goddard as 
soon as convenient. A search for qualified scientists was instituted at that time. The 
competence and talents in the particular disciplines associated with the above topics were 
not then readily available to assume these studies. 

After due consideration of the space availability at Goddard and the timing of the 
various local personnel who might be involved in this effort, it was decided to tackle the 
job with university faculty and graduate students. These academic people are usually 
free for research work during the summer months, about June 15 to September 15, and 
the Test and Evaluation Division expected to have the facilities in Building 7-10 ready for 
this purpose by June. That is the immediate practical reason for the Summer Workshop 
experiment. There is another. It stems from the firm belief that the boundaries of outer 
space knowledge will be extended by the basic research of academic institutions, and that 
NASA's greatest benefits will come from its close relationships with universities. 


OBJECTIVES AND ORGANIZATION 

A primary objective of the Summer Workshop has been to develop methods and tech- 
niques for effective utilization of the new facilities in Buildings 7-10. In order to accom- 
plish all that is implied in this statement, full advantage must be taken of the progress 
which has been made in adapting basic science to test and evaluation engineering. Another 
aim has been to interest and indoctrinate faculty members from different universities in 
Goddard’s current and future simulation problems. It is believed that goals such as these 
are attained through a common understanding of the problems involved and the combined 
efforts of academic scientists and practical engineers. The Summer Workshop was set 
up with this purpose in mind. 

A plan to relate these objectives to the six topics listed on page 1-4 was begun last 
January, when the Test and Evaluation Division was still huddled in Building 4. From 
that viewpoint, the prospects of Building 7 becoming available in spring of 1962 seemed 
most promising, for in this fine structure the Summer Workshop could pursue its studies 
in a creative atmosphere with ample working space and adequate facilities. At that time, 
however, the Ad Hoc Committee in charge of this plan urgently needed competent person- 
nel who might use this exquisite laboratory for investigating the existing problems. The 
plan was to obtain for the summer five or six professors who would act as principal in- 
vestigators in the different technical areas mentioned, and who could team up with one or 
more graduate students to tackle a specific problem. This group of university professors, 
who are specialists in definite scientific disciplines and are experienced in the conduct 
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of research work, represented the vital part of the Summer Workshop program. It was 
Goddard's good fortune to be able to procure capable and cooperative professors as well 
as graduate students. Their wholesome attitudes and unstinting efforts together with the 
actual results they achieved this summer should stimulate greater accomplishments 
among those who undertake programs of this type in the future. 

Organizational planning and preparation were most important in making effective and 
optimum use of the time during which the academic groups were to be at Goddard . It was 
essential to anticipate many of their needs and wants in advance of their arrival in Build- 
ing 7. First of all, by the end of May, there was hardly a doubt in the mind of any pros- 
pective member of the Summer Workshop as to which technical area he would be investi- 
gating. The team groupings were arranged either by correspondence, by personal 
interviews, or by telephone conversation. And it is gratifying to state that the teams, by 
and large, remained intact and functioned very well as units, even though a large major- 
ity of the individuals had not met before coming to Goddard. Then, there were ready for 
the different teams as they arrived, stacks of pertinent technical material in each of the 
various subjects to serve as a start for background reading. Thus, every man upon arri- 
val was aware of the general scientific area with which he would be concerned, he was 
given several publications in his field of study to acquaint him with the problem; and he 
was furnished a comfortable desk at which to begin his summer's work without loss of 
time and without feeling at loose ends on a new job. 

Another means which proved most helpful to the personnel of the Summer Workshop 
was the procedure for obtaining urgently needed reference materials for their use. In 
order to overcome the usual delays in normal government inter-library loans (from one 
to eight weeks, or more, to obtain a book or a periodical), a more direct approach to the 
different libraries was pursued in the case of the Workshop project. Inquiry was made 
directly to the several libraries of the local research institutions as to the availability of 
the particular book or journal that was needed. If the required reference material was 
on hand and could be loaned for a couple of weeks, the book was on the desk of the person 
who asked for it within two days from the date of the request. Indeed, this exceptional 
service which the Summer Workshop personnel enjoyed was due to a rare set of circum- 
stances, and proper advantage was taken of them. The fact remains, however, that the 
Workshop efforts were tremendously expedited by this approach; it is estimated that 
about 15 man-weeks in the overall period were saved by handling the book situation in 
this manner. 

CONCLUSIONS AND RECOMMENDATIONS 

The Goddard Summer Workshop experiment has provided some interesting conclu- 
sions. To begin with, it has been clearly established that a coordinated interchange of 
technical information between academic scientists and Goddard staff members can be most 
valuable to the solution of immediate space problems as well as to the accomplishments 
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of long-range projects. This conclusion is based on the interest shown by Goddard per- 
sonnel in the Workshop program and in its technical achievements. Also the comments 
received, oral and written, indicate that a far better understanding with the university 
groups has been effected by virtue of the communication channels which have resulted 
from this experiment. In addition, there is reason to believe that the favorable appraisal 
of the workshop effort by the Goddard technical people as well as by the visiting partici- 
pants represents an integrated critical assessment of the situation as a whole. The close 
association between the academic groups and the permanent Goddard staff provided a 
procedural flexibility that yielded desired results in less time than ordinarily would be 
required by contractual means. It is believed that this general environment stimulated 
new ideas and original thinking. 

Now we come to the research findings of the university teams. Some of these are 
practical and easily recognized as being worthwhile contributions to our problems, while 
others have given us a good start and are most promising for the future. The workshop 
studies not only have produced results but they have also encouraged a research spirit 
and a higher morale factor in Goddard personnel toward scientific approaches and ef - 
forts. Because of the Summer Workshop activities here at Goddard, it may be concluded 
that: 


1 . At least three of the workshop reports have been or are expected to be submitted 
for publication in standard scientific journals, and will probably be issued as NASA tech- 
nical notes. 

2. Two local professors associated with the workshop are continuing one or two days 
per week at Goddard on certain phases of their respective projects. These investigations 
are conducted in cooperation with and under the control of the Goddard technical staff. 

3. Three of the workshop professors have already submitted, or are about to submit, 
unsolicited proposals to the Director of Grants and Research Contracts, NASA, on certain 
follow-up studies of the Summer Workshop which they wish to undertake. 

4. A number of the graduate students who participated in the various Summer Work- 
shop projects have expressed an interest in coming to Goddard as permanent employees, 
after completing their work for advanced degrees. 

5. One of the investigations carried through this summer was sufficiently definitive 
and practical to warrant a patent application by Goddard on the idea, product, and techni- 
que developed by Messrs. John and Hilliard. 

It is therefore recommended that: 

1. Plans should be initiated at once to institute a Summer Workshop for the 1963 
season. 
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2. The following improvements should be incorporated in the 1963 Workshop on the 
basis of experience gained in the 1962 Workshop: 

• Start planning and contacting candidates at an early date, preferably before 
December 1962 

• Hold preliminary conferences with team leaders in early spring to discuss 
programs and arrive at decisions on technical approaches 

• Determine what equipment is needed and take timely action to provide it 

• Arrange rapid library services to supply references with minimum delay 
after receipt of requests 

• Provide a Goddard permanent staff member to participate full-time with 
each team 

• Solicit and encourage full participation of the Directorates of Tracking and 
Data Systems and of Space Sciences and Satellite Applications in all scien- 
tific phases of the Summer Workshop. 
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SOLAR SIMULATION 


PRELIMINARY STATEMENT 

(By Milton Schach, Thermal Systems Branch, Spacecraft Technology Division, GSFC) 


The general objective of any kind of space simulation is to provide a means for ob- 
serving the performance of devices or systems on the ground in the laboratory. The 
space environment (i.e., the atmosphere, the particle flux, and the electromagnetic radi- 
ation) is extremely complex and is not readily simulated. Gas composition, particle- 
energy distribution, and spectral-energy distribution are extremely difficult to simulate, 
singly or simultaneously. 

What are the specific areas which generate an interest in solar simulation? First 
of all, radiation from the sun is a principal factor in the thermal equilibrium of satel- 
lites; solar radiation constitutes at least 75 percent of the average energy incident on the 
spacecraft. Second, solar radiation is a source of energy for a variety of solar convert- 
ers, photovoltaic devices, thermoelectric devices, and thermionic devices. Third, the 
sun is responsible for a variety of photophysical and photochemical effects, which are of 
interest in the study of materials properties, principally the degradation of materials 
under solar radiation. The simulation of the ultraviolet spectrum is attracting consider- 
able attention. In addition, solar energy may be used as a means of active temperature 
control by taking advantage of certain of these photochemical and photophysical properties. 

The Summer Workshop will be concerned chiefly with solar radiation and its influ- 
ence on spacecraft thermal equilibrium. Let us take a quick look at some relevant prop- 
erties of solar radiation which are reasonably well known. 


In the first place, if we talk about solar radiation in the vicinity of the earth, we are 
talking about an energy of the order of about 130 watts per square foot on a surface nor- 
mal to the sun's rays. This energy is spread over a spectrum that goes to infinity in 
both directions. Happily for those of us who are concerned with thermal equilibrium, 
little more than 1 percent of this energy lies below .3 of a micron in wavelength and 
slightly more than 3 percent lies above 3 microns in wavelength. 

I should point out that the regions both above and below the limits cited are of vital 
interest in understanding the structure of the sun, and they may be extremely important 
in understanding materials degradation. 

The solar radiation may be described by comparing it to that of a blackbody. In 
order to match the energy within the bandwidth-from 0.3 to 3 microns-the best fitting 
blackbody curve would correspond to a temperature of 5800° K. 
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This energy is remarkably stable, at least in the region with which we are concerned. 
There is relatively small variation in time. (In practice, there is a variation for those of 
us who are concerned with earth satellites; namely, the variation that results when the 
earth moves around in its elliptic orbit. The changing earth-sun distances produces a 
variation of ±3 percent in the 130 watts during the year.) 

Summing up, then, solar radiation at any given distance is uniform in intensity. It is 
almost parallel, having an angle of collimation of less than one-half degree, and its 
spectral-intensity distribution curve looks roughly like that of a blackbody at 5800°K. 


The external energy sources which contribute to the thermal balance of a spacecraft 
are two, in addition to direct sunlight: (1) the reflected sunshine that the earth directs 
towards a satellite, and (2) the earth-emitted radiation. The reflected radiation varies 
rather widely in intensity, as it is a function of the orbital position and height above the 
earth’s surface. As a result, there may be instantaneous inputs of reflected solar energy 
up to 50 percent of the direct solar energy. These inputs may be as small as zero when 
the satellite is over the dark side of the earth. The earth’s emitted radiation, too, is a 
variable quantity, which can achieve values up to 15 percent of the direct solar radiation. 
The reflected solar energy and earth-emitted energy drop to less than 2 percent of the 
incident solar energy at a height of 20,000 miles above the earth’s surface. 


Finally, there is the internal-energy dissipation which must be considered. Most of 
the spacecraft launched to date dissipate less than 100 watts internally. This source is 
small compared to solar radiation, but is significant in effect. 


What is required in the way of solar simulation for thermal design and for our ther- 
mal evaluation studies? First of all, we’d like a source with a good match to the solar 
spectrum. We would like the radiation to be parallel and uniform over the spacecraft 
volume. In addition, there are spurious radiation sources which arise as a result of the 
simulated space geometry, in particular the inter-reflection between the walls of the 
chamber, the satellite, and the optical elements of the simulator. These spurious sources 
need to be minimized. 

Let us turn briefly to spectral match. The degree of closeness that is required is a 
function of the optical properties and geometry of the surface of interest. Thus, for ex- 
ample, by using a homogeneous flat plate, the need for spectral match can be eliminated 
entirely. The objective of simulation is to make certain that the spacecraft will absorb- 
on every unit area-exactly the same amount of energy in the laboratory as it will in or- 
bit. This requirement can be met by a uniform collimated source of the proper intensity. 


It turns out, however, that as soon as simulation is attempted for any elaborate 
spacecraft whose surface properties vary, this kind of simulation proves quickly inade- 
quate, so that it becomes necessary to achieve an improved spectral match. 


The collimation requirement is important, again depending upon the complexity of 
the satellite. If there is poor collimation-that is to say, if there is divergent sunshine- 
then parts of the spacecraft which would not normally be illuminated are indeed illumi- 
nated, and, conversely, some parts of the satellite will receive less illumination than 
from a collimated source. 
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Uniformity is important because of the relationship between incident solar energy 
and the temperature. It may be worth while to take a moment to indicate what the effects 
of simulation errors are. All the factors-spectral match, collimation, uniformity, the 
inter-reflection between the satellites and the simulating source and the walls-act to 
produce a certain absorbed energy at the surface of the spacecraft. If, for any reason, 
the energy differs from the true absorbed solar energy, then there is an error in space- 
craft temperature. As an example, consider an element of area which is in equilibrium 
with solar radiation. The equation for equilibrium shows that the temperature varies as 
the one -fourth power of the absorbed energy. For a surface in equilibrium near room 
temperature, an error of 10 percent in absorbed energy would result in a maximum tem- 
perature uncertainty of 7.5° C. 


One general comment may illustrate a point made earlier. Solar simulation is an 
area that involves more than one discipline; partly because of this, it is an area in which 
progress is slow. It requires a knowledge of the spectrum of gases and other radiation 
sources. It requires a knowledge of optical systems. It requires a knowledge of a vari- 
ety of technologies that go along with these disciplines. 

It is therefore our hope that this group can in some measure contribute suggestions 
and ideas to implementing this technology. 


Thank you. 
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SPECTRAL ENERGY DISTRIBUTION 
OF A MERCURY-XENON LAMP 

by 

M. P. Thekaekara 

ABSTRACT 

This report presents the results o f the work done by the 
Summer Workshop Team A* Test and Evaluation Division, God- 
dard Space Flight Center. 

A critical analysis has been made of the problems of solar 
simulation for studying the thermal balance of spacecraft, in par- 
ticular of the spectral-energy distribution of the high-pressure 
mercury-xenon arc as compared to that of solar radiation in 
outer space. The results of previous investigators have been 
analyzed. A technique has been developed for precision meas- 
urement of intensity, using a high-dispersion, high-resolution 
monochromator. The spectral characteristics of a 2.5-kw 
mercury-xenon lamp have been studied in the visible region. 
The data are presented in a series of charts and tables. The 
applicability of mercury-xenon arcs for solar simulation has 
been examined with the aid of the relevant experimental data. 


INTRODUCTION 

The GSFC spacecraft programs at present are the earth-orbiting scientific payloads, 
in particular, the Nimbus, the Orbiting Geophysical Observatory, and the Orbiting Astro- 
nomical Observatory. To test them prior to launching, a space-environment simulator is 
being constructed. The radiation in the simulator will be that of 127 mercury-xenon 
lamps. Many interesting and urgent problems are connected with solar simulation, such 
as the requirements of energy balance, the dependence of thermal balance on the shape 
of the orbit, the optics of the simulator, etc. Some of these problems are so vast and 
complex that no attempt at an adequate solution can be made within the brief period 
allotted to the Summer Workshop. 

The spectral distribution of the mercury-xenon lamp is a problem of immediate ur- 
gency. There is no longer any question as to the relative merits of the mercury-xenon 
lamp, the carbon arc, the laser, the plasma jet, and other possible sources; we are com- 
mitted to the mercury-xenon lamp. Many attempts have been made to determine the 


Members of the Summer Workshop Team "A* were Matthew P. Thekaekara, Ralph A. Goodwin, Graham D. Gutsche, 
Anthony Cerkanowicz and Bruce Friedman. 
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spectral-energy distribution of these lamps and the results are highly at variance with 
each other. 

Nor can the question be dismissed as of relatively small importance; the spectral 
distribution is highly significant for the problem of thermal balance. The total energy 
received by the test vehicle in the simulator is an integral of the absorption coefficient 
of the surface coating and of the spectral-energy distribution. If the absorption coeffi- 
cient is uniform throughout the spectral range, spectral matching is of no importance; it 
is sufficient that the total energy incident per unit area be equal to that in outer space. 
However, experiments show that, for practically all types of surface coating, the absorp- 
tion coefficient varies widely from zero to 100 percent and is highly sensitive to the spec- 
tral range. Hence arises the necessity for a spectral match between the solar energy and 
the energy of the artificial source to be used for the simulator. 

The energy due to the sun outside the earth's atmosphere is known to a high degree 
of accuracy. Its spectral distribution is represented, probably within 1 percent, by the 
well-known Johnson curve. The planet radiation due to the earth is relatively small and 
is also sufficiently well-known. The albedo due to the earth is a highly varying factor 
for which no perfect simulation is possible, but whose spectral characteristics are rather 
close to those of the solar radiation. 

The American Institute of Electrical Engineers has laid down certain specifications 
for solar simulators: namely, that throughout the 3,000 - 12,000A range, the total energy 
in any wavelength band shall not deviate from the Johnson solar-distribution curve by 
more than the following percentages in every square centimeter area of the test region. 
For 200A band, 20 percent; for 400A band, 15 percent; for 600A band, 10 percent; and for 
1,000A band, 5 percent. These specifications are very restrictive and cannot be attempted 
except for small test areas. In any attempt for spectral match, a compromise must be 
made between the size of the test area, the cost of the facility, and the spectral match 
specifications. However, it is important to know what degree of mismatch there will be 
in any actual simulator. 

In problems of thermal balance, the major factor is not the energy distribution but 
the total amount of energy absorbed by the spacecraft. By adjusting the input wattage of 
the lamp, it is possible to equalize the energy absorbed per unit area by the test vehicle 
in the simulator and by the spacecraft in outer space. In doing so, it is essential to know, 
first, the variation of absorption coefficient with wavelength, and, second, the energy dis- 
tribution of the mercury-xenon lamp. Extensive studies have been made on the absorp- 
tion coefficients of surface coatings. The techniques are simple; the results are repro- 
ducible. On the other hand, the data concerning the spectral distribution of the mercury- 
xenon lamp leave a great deal to be desired. 
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A search through available literature shows that there are five sets of independent 
data on the spectral-energy distribution of the mercury-xenon lamp. For the sake of 
convenience these data are designated as lamps 1, 2, 3, 4 and 5. 

The most complete data (lamp #1) seem to be those given by the Westinghouse Lamp 
Division, Westinghouse Electric Corporation, Bloomfield, N. J., in a monograph signed 
by J.C.H. (1962). The table gives the spectral-energy distribution of the 2.5-kilowatt Hg- 
Xe short-arc lamp, SAHX-2500B, in watts radiated per 100A band, in the wavelength 
range 0.2 m - 2.0m- The total energy radiated in this whole range is 1485.6 watts. It is 
known that beyond 2.0m the spectral irradiance decreases exponentially, and becomes 
practically zero at 5.0m- Two different methods of extrapolation for the radiant energy 
in the range 2.0m - 5.0m gave the same value, 243 watts. Thus, the total radiant energy 
from the 2.5-kw lamp in the whole spectral range 0.2m - 5.0m is 1729 watts. This is 
taken as a figure to which all other spectral-distribution data are to be normalized. 

A second and fairly complete source of information is a paper published by the Han- 
ovia Company on the Hg-Xe lamps and the Xe lamps. The Hg-Xe data (lamp #2) are on a 
1000-watt lamp, in the range 0.2m - 1.4m, on a percentage basis. The authors state that 
the energy output in this range is 480 watts. Assuming that the luminous efficiency is the 
same for the 1-kw and 2.5-kw lamps, the energy from the 2.5-kw lamp in the range 0.2m 
- 1.4m is found to be 1200 watts = 480 x 2.5. The energy output of lamp #1 in the same 
range is 1217.6 watts. The close agreement between these two values shows that a mean- 
ingful comparison can be made between lamps #1 and #2. Assuming the output of a 2.5-kw 
lamp in the range 0.2 - 1.4m to be 1217.6 watts, the output in watts per 100A band can be 
calculated from the percentage values supplied by Hanovia. 

An independent set of data (lamp #3) is available in a paper published by Mann and 
Dubey. 1 The authors give a histogram comparing solar radiation and the Hg-Xe lamp 
radiation, on a percentage basis, for each 500A band in the range 0.3 - 1.3m. Assuming 
that in this range the total output is the same as that of lamp #1 (namely 1141.8 watts), 
it is possible to prepare a table of energy output in watts per 500 A band for lamp #3. 

Minneapolis-Honeywell performed a series of tests on three Hanovia lamps of 2.5- 
kw rating, operating them under seven different conditions of input wattage and of lamp 
age. The data are presented in a series of tables, on a percentage basis, for bandwidths 
of two or more thousands of angstroms. These values have been averaged (normalized 
to the same output as lamp #2) and values for energy output in watts per 500A band have 
been calculated. These data are referred to as data on lamp #4. 

Minneapolis-Honeywell also made a series of four tests on a 2.5-kw Westinghouse 
lamp #4-21, operated at 2.5 kw, to determine what change, if any, occurs in spectral dis- 
tribution due to aging of the lamps. Though the data are not quite conclusive, they seem 
to indicate that the energy in the range 0.5 - 0.7m increases with the age of the lamp. 
Averages of the four tests were made, the totals were normalized to the same output as 
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lamp #1, and the values of energy output in watts per 500A band were calculated. These 
data are referred to as data on lamp #5. 

In order to compare these data with the solar irradiance, the latter should be nor- 
malized to the same total energy. The solar constant (the energy received per cm 2 out- 
side the earth's atmosphere) is 0.140 watts, of which 99.5 percent is in the wavelength 
range 0.2 - 5^; therefore, to compare the solar constant with anHg-Xe lampof 1729-watt 
energy output, it is necessary to calculate the energy received from the sun over an area 
12,350 cm 2 . 

(Solar irradiance data available in literature 2 give values at very close intervals in 
watts per cm 2 of surface per micron-wavelength band. These values multiplied by 123.5 
give the solar flux in watts per 100A over an area receiving a total energy of 1729 watts.) 

Table 2-1 lists the results of these calculations. The energy values are in watts, 
integrated over 500A bandwidths. These results are also presented in Figure 2-2. A 
more detailed comparison of solar irradiance with that of the Hg-Xe lamp is presented 
in Figure 2-1. The solar irradiance curve in Figure 2-1 is the Johnson curve referred 
to earlier; the curves for the Hg-Xe lamp are based on data from lamps #1 and #2. The 
ordinates are energy values integrated over 100A bandwidths. 


Table 2-1 

Spectral- Energy Distribution of Hg-Xe Lamps: 
Comparison of Data Available in Literature. 


Wavelengths 

(Microns) 

Energy Radiated in 500 Angstrom Interval by Each Lamp (Watts) 

Solar 

#1 

#2 

#3 

#4 

#5 

.3 - .35 

59.1 

442.8 

77.5 

131.3 

33.9 

24.9 

.35 - .40 

75.5 

70.5 

91.0 

154.1 

39.8 

41.9 

.40 - .45 

119.2 

109.5 

97.1 

134.7 

42.5 

63.8 

.45 - .50 

131.0 

26.9 

35.1 

52.9 

13.1 

15.7 

.50 - .55 

120.0 

73.5 

117.9 

94.8 

132.1 

123.0 

.55 - .60 

115.5 

95.3 

146.1 

126.7 

163.7 

159.5 

.60 - .65 

105.9 

36.5 

47.4 

42.2 

53.1 

61.1 

.65 - .70 

95.4 

39.7 

50.6 

50.2 

56.6 

66.4 

.70 - .75 

83.5 

36.7 

49.0 

48.0 

46.4 

22.9 

.75 - .80 

73.8 

37.6 

37.3 

35.4 

35.3 

23.4 

.80 - .85 

65.5 

41.6 

43.0 

28.5 

40.7 

25.9 

.85 - .90 

58.5 

72.1 

45.1 

44.5 

42.7 

44.9 

.90 - .95 

52.5 

76.1 

48.2 

37.7 

45.6 

47.4 

.95 - 1.00 

47.3 

85.0 

49.7 

26.3 

97.0 

53.0 

1.00 - 1.05 

42.9 

94.6 

51.7 

30.8 

85.1 

80.2 

1.05 - 1.10 

39.2 

50.4 

26.9 

25.1 

44.3 

31.8 

1.10 - 1.15 

35.3 

50.6 

34.9 

26.3 

57.5 

42.9 

1.15 - 1.20 

32.5 

34.5 

22.2 

18.3 

36.6 

29.3 

1.20 - 1.25 

29.8 

37.1 

23.9 

14.8 

39.4 

31.4 

1.25 - 1.30 

26.6 

30.9 

18.8 

18.3 

31.0 

26.2 

1.30 - 1.35 

24.0 

30.5 

17.8 


29.3 

25.9 

1.35 - 1.40 

21.5 

30.7 

20.9 


34.4 

26.0 
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Figure 2-2 - Comparison of spectral -distribution data from different Hg-Xe lamp sources 

These five sets of data present considerable variation in the energy values, in par- 
ticular the wide divergence between lamps #1 and #2, and lamps #4 and #5. In the ultra- 
violet range the lamps #4 and #5 indicate considerably higher energy. There is an over- 
all agreement between the data on lamps #1 and #2, as shown by the ratio of the energy 
in the line spectrum to the total energy of the lamp in the range 0.2 - 1.4/u ( 36 percent 
from data on lamp #1, and 37 percent from data on lamp #2). 

As stated earlier, the AIEE specification on spectral matching for solar energy is 
that, for any 200A band, the energy shall not deviate from the Johnson curve by more 
than 20 percent. In ten regions of the spectrum chosen at random, each of 200A width, 
the Johnson curve was compared to the Westinghouse data (lamp #1) and to the Hanovia 
data (lamp #2) . 

The results of this analysis are given in Table 2-2; the average deviation is consid- 
erably greater than the AIEE specification (63% for lamp #1 and 95% for lamp #2). A 
similar comparison made for the 1000A band shows that the average deviation is 41% 
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according to the Westinghouse data and 44% according to the Hanovia data, compared to 
the AIEE specification of 5%. 


Table 2-2 

Comparison of Spectral Irradiance of Hg-Xe Lamps #1 
and #2 with the Johnson Curve 


Wavelength Range 
(Microns) 

Percentage Deviation from Johnson Curve 

Westinghouse Data, 
Lamp #1 
(%) 

Hanovia Data, 
Lamp #2 
(%) 

.21 - .23 

+ 111 

+ 630 

.33 - .35 

- 48 

- 5.3 

.4 - .42 

- 6.8 

- 12 

.58 - .60 

- 47 

- 58 

.76 - .78 

- 48 

- 51 

.87 - .89 

+ 34 

- 8.3 

1.00 - 1.02 

+ 240 

+ 94 

1.12 - 1.14 

+ 24 

+ 46 

1.23 - 1.25 

+ 18 

- 29 

1.36 - 1.38 

+ 55 

- 19 

Average Deviation, 63% 

95% 


The large differences between the data on lamps #1 through #5 raise an important 
question for solar simulation. Not only is there a lack of spectral match; there is no 
reliable information on the extent of the spectral mismatch. The question arises whether 
the mercury-xenon lamp is something that changes its color every time one looks at it. 
At any rate, it is important to know the spectral distribution of the energy of these 127 
lamps which are going to illumine the test floor of the gigantic simulator. 

More than 40 percent of the solar energy is in the visible range from 3500 to 7000A; 
nearly as much of the mercury-xenon lamp energy is also in the same narrow- wavelength 
range, where absorption coefficients of materials vary greatly. It is also the range in 
which the mismatch between the mercury-xenon lamp and the solar radiation is most 
pronounced. In the long infrared range, the spectral mismatch is less striking. In the 
ultraviolet range, the energy is relatively small. 

Team A of the Summer Workshop decided to measure the spectral energy distribu- 
tion in this visible range, and to develop a technique of high-precision intensity measure- 
ment, using a grating spectrograph rather than a prism spectrograph in order to obtain 
the necessary high resolving power and spectral dispersion. Very small slitwidths (20 
microns) were used, as the spectral range viewed by the instrument at any given time 
should be small compared to the half -widths of the lines of mercury and xenon. 
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THEORETICAL CONSIDERATIONS 

The essential information needed about the mercury-xenon lamp is the amount ot 
radiant energy falling per unit time per unit area on the test floor of the simulator. If 
the surface of the spacecraft were an ideal blackbody surface, having an absorption coef- 
ficient unity throughout the wavelength range, or an ideal greybody surface, having an 
absorption coefficient constant and independent of wavelength, it would be sufficient to 
know the total energy incident on the test floor, per unit area per unit time. The total 
energy can be expressed as an integral 

E = (1) 

The limits of integration are strictly speaking zero and infinity, but for all practical pur- 
poses the limits are 0.2^ to 5.0 /x. E is the total energy and is the area under the energy- 
distribution curve of the monochromatic energy E^ versus wavelength k. In wavelength 
ranges below 0.2 /x and above 5.0^, E^ is practically zero for both the solar radiation and 
the mercury-xenon lamp. 

For purposes of solar simulation, not only must the total energy be known, but also 
the energy distribution function E^, since the surface coatings vary considerably from 
the ideal conditions of the perfect black or perfect grey. The absorption coefficient a 
varies between 0 and 1 and is a function of the wavelength. The total energy absorbed 
per unit area of the test body is 

5.0m 

E t = I (2) 

J 0. 2m 

In order to draw meaningful results from the solar-simulation tests, the energy E t 
absorbed by the test surface should be equal to the energy E s absorbed by the surface 
of the spacecraft in outer space. 

E s is the sum of three integrals expressing, respectively, the solar radiation direct- 
ly from the sun, the solar radiation reflected from the earth, and the planet radiation 
from the earth. 

Assume for the present that E s (the radiant energy received by the spacecraft) and a 
(the absorption coefficient of the surface coating) are known to a sufficient degree of 
accuracy. 

A few simple considerations will show how the total energy E t and the monochro- 
matic radiant energy E^ on the test floor differ from the same quantities as emitted by 
the mercury-xenon lamp. The energy emitted by the lamp is collimated and directed to 
the test floor by a complex optical system, consisting of three quartz lenses and three 

aluminized reflecting surfaces. Each of these reduces the total available energy by a 
certain percentage, and also changes the spectral energy distribution to a small extent. 
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There is in particular a noticeable effect in the wavelength range 0.75 - 0.95m, where the 
reflection coefficient of evaporated aluminum drops to about 0.8. In the rest of the spec- 
tral range, the reflection coefficient has a practically constant value 0.95. From the 
known reflection and transmission properties of the optical system, it is possible to cal- 
culate the total energy E t and the monochromatic energy E^ on the test floor, provided 
these quantities are known at the source. 

Methods suggested for measuring the total energy E t emitted by the mercury-xenon 
lamp can be divided into two categories: those which measure the rise in temperature of 
the liquid in a continuous-flow calorimeter, and those which measure the thermoelectric 
electromotive force in a thermocouple. All methods call for fairly elaborate calibration 
procedures for the instruments. The total energy may also be calculated from mono- 
chromatic energy distribution, since the total energy is an integral of the energy- 
distribution function. 


The absolute standard of comparison for spectral energy distribution is the ideal 
blackbody. But an ideal blackbody is impossible to construct, and all practicable approx- 
imations to the ideal blackbody are too complex and require experimental procedures 
which are too elaborate, so that secondary standards of spectral radiance must suffice. 

A secondary standard commonly employed is a tungsten-ribbon lamp. The spectral 
distribution of the tungsten- ribbon lamp can be established in terms of its color- 
temperature, the spectral emissivity of tungsten, and the absolute radiant-energy output 
for selected wavelengths in terms of an NBS standard of thermal radiation. 


The energy distribution within a blackbody enclosure is given by the Planck equation 


87Tch 


^ ~ X S e ch/\kT _ I 


dx. 


(3) 


where 0^dX is the energy density per unit volume in the wavelength range X to X + dX; c 
is the velocity of light; h is the Planck's constant; k is the Boltzman constant; and T is 
the absolute temperature. 


The energy radiated from a small orifice on the side of the blackbody enclosure per 
unit area of the radiating surface, per unit solid angle, in a direction normal to the sur- 
face, is 

ExdX = ^dXxiL. (4) 

The values of these functions for different temperatures and wavelengths are avail- 
able in literature . 3 In using these tables for precision measurement of intensity, it 
should be noted that several authors (in particular, Pivovonsky and Nagel) base the com- 
putations on the following values for the two constants 2hc 2 and hc/k: 

2hc 2 = 1.1909 x 10' 12 watt cm 2 per steradian. 
hc/k = 1.4380 cm °K. 
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If these constants were recalculated from the best available values for h, c and k, a 
slightly different value would be obtained for hc/k. The value computed from constants 
h, c andk, given in "American Institute of Physics Handbook" edited by Dwight E. Gray 
(McGraw Hill Book Company, New York, 1957) p. 7-3, is hc/k = 1.43885 cm K. Another 
good source for fundamental constants is E. R. Cohen and others in Rev. Mod. Phys. 27, 
363 (1955); using this latter source produces very nearly the same value, hc/k = 1.43879 
cm °K. This difference (nearly ten times the estimated experimental error in the deter- 
mination of h, c and k) is probably due to a slightly different value of k assumed by Piv- 
ovonsky and Nagel. 

From the known spectral ir radiance of the standard blackbody, the next step is to 
calculate the spectral irradiance of the secondary standard, the tungsten-ribbon lamp. 

For this, the temperature of the ribbon for a specified current is determined by compar- 
ing the radiation from the tungsten ribbon to that from a blackbody. From the color tem- 
perature of the ribbon (the color temperature being defined as that temperature of a black- 
body which has the same ratio of spectral radiances as the tungsten ribbon in two intervals 
to \ x + d\ x and k 2 to ^ 2 + dk 2 ) may be derived the true temperature. Next, by em- 
ploying blackbody data for this temperature and the spectral emissivities of tungsten, a 
spectral-energy curve for tunsten is calculated. Exhaustive studies are available in 
literature on the spectral emissivity of tungsten. 4 

In actual practice, it is not necessary to go through this elaborate process of deter- 
mining the color temperature and evaluating the Planckian functions; calibration charts 
may be obtained from the Radiometry Section of the National Bureau of Standards for the 
spectral irradiance of a given lamp at wavelength intervals of 500A or less. The values 
are quoted to three significant figures, and may have a maximum error of 5 percent. 

Calibration charts of the two lamps used in the present investigation are given in 
Table 2-3. 

Values of the spectral irradiance of the tungsten ribbon lamp at intermediate wave- 
lengths are determined by linear interpolation. Since l w , the spectral irradiance, changes 
by a factor of 10 4 from one end of the spectral range to the other, it is more convenient 
to plot K = log I w versus k. Let K 1} K 2 be the values of log I w at two wavelengths k v and 

^ 2 - 

Values of K c , where the subscript c denotes that it is a calculated value, are deter- 
mined from the linear equation 

K c -K x K 2 -K x ( 5 ) 

A. “ A. ^ A. 2 — A. i 

where k lt and k 2 are 4000A and 6000A, and k takes successively the values in column 1 
of Table 2-3. 
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Table 2-3 

Spectral Irradiance of Lamps U-187 and U-21 at 35 Amperes 
(in microwatts per steradian per nanometer per mm 2 of source) 



Lamp U-187 

Lamp U-21 

2500 

3.76 x 10" 3 

5.87 x 10 -3 

2600 

8.71 x 10‘ 3 

1.31 x lO" 2 

2700 

1.68 X 10- 2 

2.47 X 10" 2 

2800 

3.09 X lO' 2 

4.49 X 10" 2 

2900 

5.35 x 10* 2 

7.61 x lO" 2 

3000 

8.81 X 10* 2 

1.24 x lO" 1 

3200 

2.11 x 10* 1 

2.96 X 10" 1 

3500 

6.37 X lO" 1 

8.59 x 10' 1 

3700 

1.18 

1.55 

4000 

2.55 

3.30 

4500 

6.84 

8.69 

5000 

1.47 X 10 

1.80 x 10 

5500 

2.54 x 10 

3.05 X 10 

6000 

3.86 x 10 

4.55 x 10 

6500 

5.32 X 10 

6.25 X 10 

7000 

6.82 X 10 

7.86 x 10 

7500 

8.16 x 10 

9.21 x 10 


At each of these values is found the difference ak = k c - K s , where k s = log l w is the 
standard value from the log of i w in Table 2-3. A graph is drawn for AK versus \, which 
permits calculation of log l w at wavelength intervals as close as are needed for the 
problem on hand, by applying the correction AK to the linearly interpolated value of K c as 
given by equation (5). 

The spectral energy distribution of the mercury-xenon lamp can be determined by 
comparing its energy output in a given spectral range to that of the standard tungsten- 
ribbon lamp. It is essential that the energy output from the two sources be measured 
under identical conditions, and that the response of the recording device be strictly 
linear. 

The spectrum of the mercury-xenon lamp shows a large number of lines super- 
posed on the continuum. Since the tube operates under a high pressure of about 20 at- 
mospheres, all the lines show considerable broadening. If the spectral range covered 
by the photomultiplier tube for any given setting of the spectrograph is of a bandwidth 
comparable to the halfwidth of the spectral lines, we are essentially comparing a line 
spectrum of mercury-xenon with a continuous spectrum of tungsten. In such a case, the 
response of the photomultiplier tube is not strictly linear to the intensity, and a correc- 
tion factor must be introduced; this correction factor depends on the wavelength range, 
the slit width, and other constants involving the geometry of the spectrograph. 

On the other hand, if the slit width used is sufficiently small, and the spectral range 
covered at a given setting is therefore very small compared to the halfwidth of the lines, 
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we are essentially comparing one continuous spectrum with another. If the response of 
the recording electronics is linear, the ratio of the two spectral irradiances is equal to 

the ratio of the corresponding signals. Thus the basic equation for absolute intensity 
measurement may be written 

= ( 6 ) 
s w 

where I signifies the spectral irradiance and s a corresponding signal, and the subscripts 
H and W refer respectively to the two sources, mercury-xenon and tungsten. 

It will be helpful to recall a few basic ideas about a plane-grating monochromator in 
Ebert mounting: 

The theoretical resolving power of a grating is 



where n is the total number of rulings on the grating and m is the order of the spectrum. 

In our case we have the first order spectrum; the grating has a ruled area of 52x52 mm 2 , 
and the number of rulings is 30,000 per inch. Hence n = 61,400 lines and for X = 4000A, 
dx = .065A. The theoretical resolving power is such that two lines .065A apart from each 
other will be just resolved, and an infinitely narrow line of wavelength 4000A will have a 
half-width of .065A. 

The theoretical resolving power is rarely achieved because of imperfections in the 
grating, or lack of alignment. 

Another factor that contributes to the observed line width is the finite width of the 
slit. Let the grating be set for a particular wavelength in the neighborhood of X = 4000A 
in the first order. From the grating equation is derived 

X = d(sin a + sin /3). ^ 

where d is the grating spacing, and a and /3 are the angles of incidence and diffraction. 
Assume that the normal to the grating is inclined towards the exit slit and away from the 
entrance slit. 

From the known values of X and d, it is readily seen that (sin a + sin /3) = 0.472. 

Values of a and /3 can be calculated as follows: Let e be the angle which the bisector 
between the incident and diffracted beams makes to either of the two beams, and let y be 
the angle which the same bisector makes with the grating normal. Thus 

a = y + € and /3 = y - e; and, ( 9 ) 

sin a + sin /3 = sin(y + e) + sin(y - e) 

= 2 sin y cos e. 
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The value of e is known from the geometry of the Ebert mounting spectrograph. 
Take for example a pencil of rays which enters the spectrograph through the geometri- 
cal center of the entrance slit and reaches (after reflection, diffraction, and reflection) 
the geometrical center of the exit slit. From the entrance slit to the mirror the rays 

form a divergent cone, and from the mirror to the exit slit they form a convergent cone. 
The axes of these two cones are parallel, and the distance between the points where the 
axes meet the mirror is equal to the distance, D, between the two slits. Thus 

2e = 5 , where f = 50 cm 
F 

Hence e = .1032 radian = 5° - 55'. 

y= 13° - 43'; a = 19° - 21' and /? = 7° - 48’. 

These equations hold strictly for the limiting case where both the slits are geomet- 
rical points. In the actual case of finite slit widths, corresponding to every point on the 
entrance slit there is a point on the exit slit which satisfies the grating equation (8) for 
the same wavelength K. This is the special advantage of the curved slit in the Ebert 
mounting. 5 

Differentiating equation (8) with respect to /3, d\ = d cos /3d/3. The broadening of the 
spectral line due to the exit slit only is given by d\, which for a slit of 25 microns width 
is 

dk = 8.47 » 10' 5 x 0.9441 * 5 x 10 ’ 5 cm = 0.4A. 

The entrance slit contributes about the same linewidth which is superposed on the 
width due to the exit slit, so that the instrumental linewidth due to the slits alone is of 
the order of 0.4A. 

The emission lines of the mercury-xenon lamp have a halfwidth considerably greater 
than this, so that the present set-up is essentially comparing one continuous spectrum 
with another, and equation (6) can be applied without any elaborate correction factors. 


EXPERIMENTAL ARRANGEMENT 

The spectrograph as was stated earlier is a plane-grating monochromator in Ebert 
mounting. The experimental arrangement is shown in Figure 2-3. The light enters the 
instrument through a straight slit 25 microns wide and passes to a concave mirror, where 
it is collimated and reflected to the grating. The grating ruled on a 58 X 58 mm 2 blank, 
has 30,000 grooves per inch on a ruled area 52 X 52 mm 2 . The diffracted light, still par- 
allel, but with separate wavelengths diverging, is reflected back to the mirror, where it 
is again reflected and focussed on the exit slit. 

The wavelength of the light reaching the exit slit is changed by rotating the grating 
about its center. The grating is rotated by action of an 8-inch lever arm mounted on the 
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L.P. LAMP POWER SUPPLY ^ QUARTZ LENS 

P.P. PHOTOMULTIPLIER POWER SUPPLY F FILTER 

L.H. LAMP HOUSING S-l ENTRANCE SLIT 

E. M. EBERT-MOUNTING SPECTROGRAPH M CONCAVE MIRROR 

F. A. FAN-AIR SUPPLY G PLANE GRATING 

L LAMP S-2 EXIT SLIT 

Figure 2-3 - Experimental arrangement for the measurement of intensity 

pivot shaft and driven by a split nut on a precision screw. The linkage is designed to pro- 
vide a linear motion directly proportional to the sine of the angle of rotation, so that the 
screw rotation remains proportional to the wavelength. A simple system of gears per- 
mits the operator to select any of eight different speeds of rotation on the automatic 
drive or four different speeds on the manual drive. Both the manual and automatic drives 
are reversible. A Veeder-Root counter attached to rotating gear reads the wavelength 
correct to 2A. The eight automatic speeds available are 500, 250, 125, 50, 20, 10, 5, and 
2A per minute. 

Light is focussed on the slit by a quartz double-convex lens having an effective aper- 
ture of about F:8.5, the same as that of the spectrograph; in this way the full light- 
gathering power of the spectrograph is utilized, without letting too much stray light into 
the instrument. Care must be taken in aligning the source and the lens along the optic 
axis of the instrument. Once the exact position of the lens was found, it is rigidly fixed; 
as different sources are interchanged one for another, the only precaution to be taken is 
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that a well-focussed image of the source must be formed at the center of the entrance 
slit. The lens, with a diameter of 3.24 cm, is at a distance of 51.3 cm from the slit and 
74.8 cm from the source. 

The light from the exit slit enters a photomultiplier tube. Two different types of 
photomultiplier tubes were used during the experiments, GE tube 1P28 and EMI tube 
9592B, having spectral responses S-5 and S-10, respectively. The 1P28 tube has a rela- 
tively higher gain and narrower spectral-band response than the 9592B tube. Power for 
the phototubes is obtained from a standard phototube power supply. The output of the 
phototube is fed through a 1 -megohm resistor to a recording device. 

Two types of recording device were used. In the earlier phase of the experiment, 
the voltage across the megohm resistor was led to a vacuum-tube digital voltmeter. The 
voltmeter, supplied by the Hewlett Packard Co., had an input impedance of 11 megohms; 
response time, less than 1 second; ranging time, 0.2 to 2 sec.; accuracy, within 0.2 per- 
cent; variable sampling rate, between 5 printouts per second and one printout every 5 
seconds. The instantaneous voltage across the megohm resistor was displayed on a 
lighted panel, and the value was printed out on a strip chart at any desired rate. The 
main advantages of this recording technique over the conventional pen-recorder tech- 
nique were that it eliminated the need for reading out the recorder signals from a roll 
of paper, and that during any single scanning of the spectrum the output signals were re- 
corded over a range of several decades with an accuracy ranging between 1 and 0.1 per- 
cent. The digitizer automatically adjusted the decimal point in the voltage and had a 
range from 999 volts to .001 volt. However, if the relay for changing the decimal point 
tripped at the same time as the relay for the printout tripped, the printout was skipped. 
Thus, quite frequently during the scan of the mercury-xenon spectrum, printouts in the 
vicinity of the strong emission lines were missing from the strip chart and it is pre- 
cisely here that the values are most significant. Another disadvantage was that a great 
deal of information concerning the fine structure of the emission lines was omitted from 
the strip charts. These disadvantages could of course be overcome by scanning the 
spectrum at a very slow speed in the vicinity of the emission lines. There was also some 
uncertainty about the wavelength corresponding to any given printout; this defect was later 
corrected by adding a clock printer unit (561A clock printer) to the digitizer-printout rack. 
On the same line as the printout of the voltage, the printout of the time was also given, in 
hours, minutes and seconds. This constituted a fairly accurate recording of the wave- 
length corresponding to each printout. 

In the later phases of the experiment, the more conventional device was used of apen 
continuously recording the changing voltages over a slowly moving strip of ruled paper. 
The voltage across the megohm resistor was applied to a decade voltage divider (type 
1454A-Serial 3411, General Radio Corp.), which permitted the full-scale deflection of the 
recorder to be set at different values, 50, 25, or 5 millivolts, etc. Through the voltage 
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divider the signal was fed to a Brown Electronic strip chart recorder (model KSY 153 x 
18 (VAH) 163-26). The chart drive speed of 1 inch per minute along with a grating drive 
of 50 angstroms per minute was fully adequate for displaying the profiles of all the strong 
lines. 

The light source used for the initial alignment of the spectrograph was a small high- 
pressure mercury-arc lamp (Central Scientific Company, Serial No. 26830 Type SH). The 
power supply for this lamp consisted of a simple transformer (Hanovia type 7606, 
catalog number 5043). The lamp, which is compact and easily mounted in position, 
proved to be extremely useful for all routine calibration procedures, for testing the 
linearity of response of the phototube and associated electronics, for measuring the 
transmission coefficients of the filters, for checking the accuracy of the wavelength 
counter, etc. 

Two tungsten-ribbon lamps manufactured by the General Electric Company (Type 
30A/T24/3) were used as intensity standards. The lamps are designated U-21 and U-187 
in the calibration charts supplied by the Radiometry Section, Atomic Physics Division of 
the National Bureau of Standards. The values of spectral radiance given in the calibra- 
tion charts can be assumed to hold true for at least 40 hours of operation, provided the 
operating current never exceeds the rated value. Both the lamps were fairly new, and 
the glass walls and quartz window of the lamp showed no signs of darkening through tung- 
sten sputtering, so that the NBS values of spectral radiance could be accepted with 
confidence. 

The power supply for the standard lamp was a variac transformer shown in Figure 
2-4b. The standard lamp operates at 6 volts, 35 amp. The current was read by a sensi- 
tive ammeter connected in series with the transformer output through a calibrated shunt. 
The least division of the ammeter was one ampere and a rough estimate had to be made 
for fractions of an ampere. This was not considered satisfactory; for further work, a 
more precise method of measuring and stabilizing the current should be devised. The 
conventional arrangement is to balance on a precision potentiometer the lamp-input cur- 
rent with a known current, the latter being determined by a standard cell and a standard 
ohm. 


The mercury-xenon lamp used in these experiments was a 2500-watt lamp supplied 
by the Westinghouse Corporation, type SAHX — 2500B, referred to in the records of the 
Test and Evaluation Division, GSFC, as Lamp #R137. A summary of its hours of opera- 
tion is given in Table 2-4. The lamp was enclosed in a large rectangular box of alumi- 
num, 46 x 35 x 27 inches. 

The lamp is held in position by two clamps rigidly attached to the back side of the box. 
The light issues from a circular opening on the front side of the box, which is restricted 
by a diaphragm so that only a narrow pencil of rays emerges from the box, just sufficient 
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RECTIFIER 

Westinghouse HgXe Lamp Osram Igniter Westinghouse Silicon Welder 

R 1 37 SAHX-2500B 2X 200 1 6 Type WS 230/460 Volt. 



220 v n 


Figure 2 -4a - Power supply for Hg-Xe lamp, schematic diagram 



Figure 2 -4b - Power supply for standard lamp, schematic diagram 
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Table 2-4 

Operating History, Hg-Xe Lamp #R137 
(Westinghouse type SAHX-2500B) 


Date 

Average 

Operation 

Hours 

Cumulative 

Operation 

Hours 

Volts 

Amps 

Watts 

7/19/62 

38 

65.5 

2490 

r 1.00 

1.00 

7/31/62 

37-38.7 

65.5 

- 

4.75 

4.75 

8/1/62 

Lamp ignited and s 

.hut off for brief periods during 


the morning 




8/1/62 

40 

64.0 

2520 

3.00 

8.75 

8/10/62 

36.9 

65.5 

2420 

2.00 

10.75 

8/15/62 

- 

65.3 

- 

0.67 

11.42 

8/15/62 

38.0 

65.5 

2490 

1.67 

13.09 

8/15/62 

38.0 

65.5 

2490 

1.22 

14.31 

8/17/62 

37.5 

65.5 

2460 

1.00 

15.31 

8/17/62 

38.2 

65.5 

2500 

2.17 

17.48 

8/20/62 

38 

65.5 

2490 

1.45 

18.93 

8/21/62 

38 

65.5 

2490 

1.08 

20.01 

8/21/62 

38 

65.5 

2490 

0.58 

20.59 

8/23/62 

38 

65.5 

2490 

1.83 

22.42 


27 

56.2 

1515 

1.50 

23.92 

8/28/62 

38 

65.5 

2490 

2.33 

26.25 

9/5/62 

38 

65.5 

2490 

6.63 

32.88 


to cover the collimating lens, reducing to a minimum the hazard due to the strong ultra- 
violet radiation. Excessive heating within the lamp enclosure is avoided by allowing a 
strong air supply from below, provided by a motor-driven air blower, 0.1 hp, 1600 rpm. 
The hot air escapes through an opening above the enclosure. 

The temperature inside the enclosure was measured by two copper-constantan ther- 
mocouples, one located very close to the lamp at the power lead and the other located about 
3 inches behind the lamp. The current output of the thermocouples, calibrated in terms 
of the temperatures within the enclosure, was continuously recorded by a Minneapolis- 
Honeywell strip-chart recorder whenever the lamp was in operation. If the temperature 
showed a tendency to rise or fall, it was stabilized by the simple device of lowering or 
raising a shutter held against the intake of the air blower. By adjusting the rate of cool- 
ing, a steady state was attained, with no tendency to drift either in the temperatures re- 
corded by the thermocouples or in the current and voltage supplied to the lamp. 

The power supply for the mercury-xenon lamp is shown in Figure 2 -4a. The input 
is 440 volts, 3 -phase ac, stepped down and rectified by a Westinghouse silicon welder 
type WS rectifier unit. The rectified current is led to the lamp through an Osram igniter, 
Siemens-Schuckert, type ZX-200-lb. The input wattage of the lamp is determined by 
measuring separately the current and the voltage. 
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In the experimental setup described above, there are several different components, 
each supposed to be adequate for a precision measurement of intensity. However, the 
degree of reliability of each component and the percentage error that might be introduced 
by each had to be separately determined. The spectral radiance of the tungsten-ribbon 
lamp is known at the calibrated wavelength with a percentage error less than 5 percent; 
for most of the values, the error is presumably much less than 5 percent. On the other 
hand, because of the aging of the lamp or due to an inaccurate reading of the current in- 
put, larger and systematic errors might possibly be introduced, whose magnitude would 
be difficult to assess. The current input is a very significant factor, since the C R heat- 
ing which controls the temperature of the ribbon is proportional to the square of the cur- 
rent, and the total radiant-energy output is proportional to the fourth power of the tem- 
perature. Hence, there is a degree of uncertainty about the standard of comparison, and 
the experiment should be set up so that the probable error due to any of the other sources 
is considerably less than that due to the standard. 

The output of the photofnultiplier tube is measured by a vacuum-tube voltmeter in 
the digitizer-printout unit. A calibration of this voltmeter seemed desirable and was 
performed by using a standard voltage-input unit, which can change the input from 1000 
volts to .001 volt in stages of .01 volt, with a rated accuracy of 0.1 percent. In the ranges 
0-1 volt and 10-900 volts, the reading of the digitizer was exactly that of the input. In 
the range 2-10 volts, the readings of the digitizer were about 0.5 percent lower; in the 
range 900-990 volts the values were 0.5 percent higher. In the range 990-999 volts, the 
values were 1 percent higher. As the digitizer readings during the scanning of a spec- 
trum were well below 100 volts, this small difference could be totally ignored. 

A wavelength calibration of the monochromator was made with the aid of a type SH 
mercury-arc lamp. The Veeder-Root counter of the spectrograph had been set accurately 
by the manufacturers to read the wavelength within 1 or 2 angstroms; however, the alum- 
inum casing of the spectrograph and the lever arm controlling the grating drive are both 
subject to small changes in length, caused by fluctuations in ambient temperature, and 
for this reason the calibration would not necessarily remain constant. The mercury-arc 
lamp was focussed on the entrance slit and the reading of the digitizer was noted. The 
grating, manually rotated from the long-wavelength end to the short-wavelength end of 
the spectrum, was set successively on the peaks of the strong mercury lines, and the 
counter reading was noted for each line. The actual wavelength of the line appeared 
higher than the counter reading by about 3A; in no case was the difference greater than 
4A or less than 2A. This correction factor is of some importance in measuring the sig- 
nal for various wavelength settings of the tungsten- ribbon lamp. For the mercury-xenon 
lamp, however, the number of emission lines is sufficiently large to serve as an internal 
calibration of the spectrum chart, so that the counter readings are less significant. 

The phototube signals from the mercury-xenon lamp are considerably stronger 
throughout the spectral range than those from the tungsten lamp. The ratio of the signals 
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varies from about 50 at 6500A to about 5000 at the short-wavelength end of the spectrum, 
3500A, and at the peaks of the mercury lines. Thus, a relatively small error in the data 
from the standard lamp would introduce a very large error in the intensity of the Hg-Xe 
spectrum. The error can be minimized by reducing the Hg-Xe illumination on the en- 
trance slit by a known constant ratio, through the use of a wire-gauze filter, a rotating 
sector, or any other type of achromatic filter. The present experiment used neutral fil- 
ters, manufactured by Bausch and Lomb, which consist of thin aluminum films deposited 
on glass plates. Three grades of filters were available which reduce the intensity by a 
factor of 25, 10, and 5. These filters are respectively marked 1.4, 1.0 and 0.7, which are 
the logarithm to the base 10 of the intensity-reduction factor. By interposing one or more 
of these filters in front of the entrance slit, the signal from the phototube can be brought 
to a value very close to that given by the standard lamp in the same region. Several 
measurements were made on the transmission coefficients of the filters; most of the fil- 
ters were found to have a transmission coefficient very close to the rated value. All 
tests for any dependence of the transmission coefficient on the wavelength were negative. 
In a particular case (filter #1 which was most used in these experiments), the average 
value of the ratio of the voltage output without filter and with filter was 9.93; the value 
was recorded at 9 wavelengths between 4000A and 6500A. The average deviation of the 
values about the mean was 0.18, which corresponds to about 2 percent. 

The aluminum-film filters also gave a ready means of checking whether the voltage 
output on the digitizer or the recorder is directly proportional to the light input at the 
entrance. The linearity of response of the photomultiplier and the associated electronics 
is a necessary condition for absolute-intensity measurements, in terms of the tungsten- 
standard lamp. The tests showed that, within the limits of accuracy of the digitizer unit, 
the response was linear and independent of wavelength, except for very small and very 
large output signals. For very small output signals, the percentage error in the milli- 
volt reading is too high; thus, a reading of 4 millivolts signifies any value between 4.5 
and 3.5 millivolts and has a possible maximum error of 12 percent. At very high volt- 
ages, say 20 volts and above, the current through the megohm resistor is 20 microamps. 
This photoelectric current is generated by light incident on an area 25 microns x 4 mm 
of the cathode surface. If the whole cathode surface of 25 mm x 8 mm were equally il- 
lumined, the total current output would be about 40 milliamps, which is considerably 
above the maximum rated output of the phototube. Hence, at such high values of the sig- 
nal, the linearity law of the phototube response breaks down completely. Below 10 or 15 
volts the output of the phototube is strictly linear with light input. The dark current re- 
mained at zero most of the time and at one millivolt occasionally. 

Several sets of measurements were made on signal strengths from both the standard 
tungsten lamps, U-21 and U-137, using both the phototubes 1P28 and 9592B. The wave- 
lengths were set manually, with due corrections for the counter readings. The set of 
values for signal strengths showed a good consistency even though the measurements 
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were made at intervals of several days. The average deviation between different sets of 
readings was of the order of 1 percent. 


The earlier recordings of the signal strength from the tungsten lamp were made with 
the image of the tungsten ribbon focussed on the full height of the entrance slit. When the 
tungsten lamp was replaced later by the mercury-xenon lamp, it was found that the image 
of the mercury-xenon arc had a height of only about 4 mm, considerably less than the 
height of the slit. The entrance slit was therefore masked to less than 4-mm height. 
Table 2-5 shows how the signal strengths are changed by masking the slit: columns 2 
and 3 give the signal strengths (average of over ten printouts at each wavelength) for the 
full slit and masked slit respectively; and column 4, the ratio of the two values. The 
ratio is fairly constant. The average of the ratio is .2412, and the average deviation is 
.00236, about 1 percent. The logarithm of this value is a correction factor to be sub- 
tracted from the values of log s w . 

For all later recordings of the signal strengths from the standard tungsten lamp, the 
masked entrance slit of height less than 4 mm was used, so that no correction factor was 
needed. 


Table 2-5 

Comparison of Signal Strengths for Masked Slit 
and Full Slit, Using U-21 Tungsten Lamp and 
1P28 Phototube 


Wavelength 

Sj Full Slit 

S 2 Masked Slit 

S 2 /Sj 

6075A 

.6994 

.1656 

.237 

5825 

1.189 

.2816 

.237 

5575 

1.491 

.3575 

.239 

5325 

2.196 

.5355 

.244 

5075 

1.774 

.4290 

.242 

4825 

1.332 

.3250 

.244 

4575 

.9486 

.2275 

.240 

4325 

.5888 

.1428 

.243 

4075 

.3179 

.0775 

.244 

3825 

.1453 

.0352 

.242 


RESULTS 

The spectral distribution of energy of the mercury-xenon lamp as obtained from one 
series of experiments is given in Table 2-6. The intensity values are on the logarithmic 
scale. The value of log i H was obtained from the expression 

log l H = log s H + (log I w - log s w ) , 

Values of log l w , for the standard lamp for wavelengths at 25A intervals, were obtained 
from the NBS calibration table, using the method of linear interpolation. The signals s w 
were measured at intervals of 125A; values log s w for wavelengths at smaller intervals 
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Table 2-6 

Spectral Irradiance of Hg-Xe Lamp U-187 with Photomultiplier Tube 1P28 
(l H in microwatts per steradian, per manometer, per mm 2 of source, 
in a direction normal to the axis of the lamp) 





Ena 

4103 

3.25 

3.45 

4107 

3.32 

3.48 

4111 

3.33 

3.58 

4116 

3.28 

3.42 

4120 

3.26 

3.37 

4124 

3.22 

3.32 

4128 

3.20 

3.27 

4132 

3.17 

3.24 

4136 

3.16 

3.20 

4141 

3.14 

3.18 

4145 

3.12 

3.16 

4149 

3.10 

3.15 

4153 

3.08 

3.15 

4157 

3.07 

3.13 

4161 

3.06 

3.12 

4166 

3.06 

3.12 

4170 

3.06 

3.12 

4174 

3.05 

3.12 

4178 

3.11 

3.13 

4182 

3.10 

3.32 

4186 

3.10 

3.12 

4191 

3.18 

3.12 

4195 

3.10 

3.13 

4199 

3.16 

3.14 

4203 

3.11 

h a 

4217 

3.11 

Hal 

4219 

3.11 


4420 

4424 

4428 

4432 

4436 

4441 

4445 

4449 

4453 

4457 

4461 

4466 

4470 

4474 

4478 

4482 

4486 

4491 

4495 

4499 

4503 

4507 

4511 

14516 

4520 

4524 

4528 

4532 

4536 

4541 

4545 

4549 

4553 

4557 

4561 

4566 

4570 

4574 

4578 

4582 

4586 

14591 

4595 

4599 

4603 

4601 

4611 

4616 

14620 


4624 

4628 

4632 

4636 

4641 

4645 

4649 

4653 

4657 

4661 

4666 

4670 

4674 

4678 

4682 

4686 

4691 

4695 

4699 

4703 

4707 

4711 

4716 

4720 

4724 

4728 

4732 

4736 

4741 

4745 

4749 

4753 

4757 

4761 

4766 

4770 

4774 

4778 

4782 

4786 

4791 

4795 

4799 

4803 

4807 

4811 

4816 

4820 

4824 
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Table 2-6 (Continued) 



bg l H 

K 

bg l H 

K 

bg l H 

K 

bg i H 

K 

bg l H 

K 

bg l H 

K 

bg X H 

5032 

3.18 

5245 

3.17 

5457 

4.58 

5670 

3.32 

5878 

3.49 

6086 

3.20 

6295 

3.28 

5036 

3.19 

5249 

3.16 

5461 

4.79 

5674 

3.34 

5882 

3.58 

6091 

3.19 

6299 

3.29 

5041 

3.03 

5253 

3.16 

5466 

4.86 

5678 

3.34 

5886 

3.56 

6095 

3.19 

6303 

3.29 

5045 

3.27 

5257 

3.17 

5470 

4.72 

5682 

3.35 

5891 

3.62 

6099 

3.19 

6307 

3.30 

5049 

3.13 

5261 

3.16 

5474 

4.34 

5686 

3.37 

5895 

3.46 

6103 

3.19 

6311 

3.30 

5053 

3.19 

5266 

3.17 

5478 

4.36 

5691 

3.38 

5899 

3.31 

6107 

3.19 

6316 

3.30 

5057 

3.16 

5270 

3.17 

5482 

4.22 

5695 

3.40 

5903 

3.82 

6111 

3.19 

6320 

3.30 

5061 

3.15 

5274 

3.17 

5486 

3.98 

5699 

3.42 

5907 

3.37 

6116 

3.21 

6324 

3.29 

5066 

3.15 

5278 

3.17 

5491 

3.84 

5703 

3.44 

5911 

3.35 

6120 

3.29 

6328 

3.30 

5070 

3.14 

5282 

3.17 

5495 

3.71 

5707 

3.46 

5916 

3.34 

6124 

3.26 

6332 

3.31 

5074 

3.05 

5286 

3.17 

5499 

3.64 

5711 

3.52 

5920 

3.33 

6128 

3.19 

6336 

3.30 

5078 

3.05 

5291 

3.18 

5503 

3.58 

5716 

3.57 

5924 

3.33 

6132 

3.18 

6341 

3.31 

5082 

3.05 

5295 

3.18 

5507 

3.59 

5720 

3.61 

5928 

3.32 

6136 

3.19 

6345 

3.31 

5086 

3.05 

5299 

3.18 

5511 

3.50 

5724 

3.66 

5932 

3.32 

6141 

3.19 

6349 

3.31 

5091 

3.05 

5303 

3.18 

5516 

3.50 

5728 

3.72 

5936 

3.30 

6145 

3.19 

6353 

3.31 

5095 

3.04 

5307 

3.18 

5520 

3.46 

5732 

3.79 

5941 

3.30 

6149 

3.20 

6357 

3.32 

5099 

3.04 

5311 

3.12 

5524 

3.43 

5736 

3.86 

5945 

3.29 

6153 

3.20 

6361 

3.33 

5103 

3.04 

5316 

3.11 

5528 

3.41 

5741 

3.96 

5949 

3.29 

6157 

3.20 

6366 

3.33 

5107 

3.05 

5320 

3.11 

5532 

3.40 

5745 

3.80 

5953 

3.28 

6161 

3.10 

6370 

3.32 

5111 

3.05 

5324 

3.12 

5536 

3.37 

5749 

4.18 

5957 

3.27 

6166 

3.21 

6374 

3.33 

5116 

3.05 

5328 

3.13 

5541 

3.36 

5753 

4.34 

5961 

3.20 

6170 

3.21 

6378 

3.34 

5120 

3.04 

5332 

3.13 

5545 

3.34 

5757 

4.51 

5966 

3.26 

6174 

3.21 

6382 

3.34 

5124 

3.15 

5336 

3.19 

5549 

3.34 

5761 

4.71 

5970 

3.26 

6178 

3.23 

6386 

3.34 

5128 

3.26 

5341 

3.15 

5553 

3.21 

5766 

4.87 

5974 

3.25 

6182 

3.23 

6391 

3.34 

5132 

3.15 

5345 

3.16 

5557 

3.21 

5770 

4.87 

5978 

3.24 

6186 

3.23 

6395 

3.33 

5136 

3.17 

5349 

3.18 

5561 

3.30 

5774 

4.70 

5982 

3.24 

6191 

3.24 

6399 

3.35 

5141 

3.16 

5353 

3.20 

5566 

3.29 

5778 

4.52 

5986 

3.24 

6195 

3.24 

6403 

3.37 

5145 

3.16 

5357 

3.26 

5570 

3.28 

5782 

4.61 

5991 

3.23 

6199 

3.25 

6407 

3.37 

5149 

3.16 

5361 

3.42 

5574 

3.27 

5786 

4.84 

5995 

2.95 

6203 

3.24 

6411 

3.38 

5153 

3.15 

5366 

3.47 

5578 

3.27 

5791 

4.84 

5999 

2.13 

6207 

3.24 

6416 

3.38 

5157 

3.15 

5370 

3.34 

5582 

3.26 

5795 

4.67 

6003 

3.21 

6211 

3.25 

6420 

3.37 

5161 

3.15 

5374 

3.26 

5586 

3.26 

5799 

4.49 

6007 

3.20 

6216 

3.24 

6424 

3.38 

5166 

3.15 

5378 

3.26 

5591 

3.27 

5803 

4.41 

6011 

3.20 

6220 

3.26 

6428 

3.38 

5170 

3.16 

5382 

3.26 

5595 

3.26 

5807 

4.37 

6016 

3.20 

6224 

3.26 

6432 

3.39 

5174 

3.02 

5386 

3.26 

5599 

3.59 

5811 

4.30 

6020 

3.20 

6228 

3.27 

6436 

3.39 

5178 

3.17 

5391 

3.28 

5603 

3.26 

5816 

4.21 

6024 

3.19 

6232 

3.27 

6441 

3.30 

5182 

3.17 

5395 

3.39 

5607 

3.26 

5820 

4.13 

6028 

3.19 

6236 

3.28 

6445 

3.41 

5186 

3.02 

5399 

3.33 

5611 

3.26 

5824 

3.99 

6032 

3.19 

6241 

3.27 

6449 

3.41 

5191 

3.15 

5403 

3.35 

5616 

3.26 

5828 

3.94 

6036 

3.20 

6245 

3.17 

6453 

3.41 

5195 

3.15 

5407 

3.34 

5620 

3.27 

5832 

3.88 

6041 

3.21 

6249 

3.25 

6457 

3.42 

5199 

3.18 

5411 

3.34 

5624 

3.26 

5836 

3.82 

6045 

3.21 

6253 

3.26 

6461 

3.43 

5203 

3.03 

5416 

3.45 

5628 

3.26 

5841 

3.76 

6049 

3.21 

6257 

3.26 

6466 

3.43 

5207 

3.14 

5420 

3.49 

5632 

3.26 

5845 

3.71 

6053 

3.21 

6261 

3.27 

6470 

3.44 

5211 

3.15 

5424 

3.44 

5636 

3.36 

5849 

3.66 

6057 

3.21 

6266 

3.26 

6474 

3.44 

5216 

3.14 

5428 

3.47 

5641 

3.27 

5853 

3.62 

6061 

3.21 

6270 

3.26 

6478 

3.46 

5220 

3.14 

5432 

3.53 

5645 

3.27 

5857 

3.58 

6066 

3.21 

6274 

3.27 

6482 

3.44 

5224 

3.15 

5436 

3.62 

5649 

3.28 

5861 

3.55 

6070 

3.23 

6278 

3.27 

6486 

3.54 

5228 

3.15 

5441 

3.74 

5653 

3.28 

5866 

3.53 

6074 

3.24 

6282 

3.28 

6491 

3.48 

5232 

3.17 

5445 

3.91 

5657 

3.28 

5870 

3.51 

6078 

3.21 

6286 

3.28 

6495 

3.46 

5236 

3.15 

5449 

3.98 

5661 

3.29 

5874 

3.48 

6082 

3.21 

6291 

3.28 

6499 

3.46 

1 5241 

3.16 

5453 

4.38 

5666 

3.30 
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of 25A were also obtained by the method of linear interpolation. A graph was plotted for 
log l w - log s w versus wavelength. This graph served as the calibration curve for any given 
set of measurements on the mercury-xenon lamp, for conversion of the digitizer voltages 
into absolute intensity values. 

The values for s R were printed out for approximately every 2A; two successive print- 
outs were averaged on each of two different runs, and the average of both runs was taken. 
Printouts which were skipped because of hunting decimals were supplied by interpolation. 
The wavelengths were calculated from the time indicated by the clock printer. The val- 
ues of s H thus represent the average over a small-wavelength range of 4A in most cases, 
and 5 A in a few cases. 

The units for i H are the same as those of I w , namely, microwatts per steradian per 
nanometer per mm 2 of source. 

The series of experiments whose results are presented in Table 2-6 was made with 
the photomultiplier tube 1P28, and with the standard lamp U-187. Another series of ex- 
periments was made with the same standard lamp, but with the phototube 9592B. Two 
more series of runs were made using the same two phototubes but a different standard 
lamp, U-21. Extensive tabulations have been made of the results of these three addi- 
tional series of runs; these lengthy tables have not been included in the present report, 
because they do not add significantly to the foregoing analysis of the spectrum. Their 
main value is that they serve as a check on the internal consistency of the results. 

Upon comparison of the four series of runs, or different runs in the same series, 
the total light output of the mercury-xenon lamp is seen to change considerably. A small 
increase in the wattage input of the lamp does not necessarily mean an increase in the 
light output. The lamp operates steadily for long intervals of time and there is no indi- 
cation of the energy distribution shifting from the line spectrum to the continuum and 
vice versa, or from one spectral range to another. In other words, the total intensity 
goes up or down without any noticeable change in the relative spectral distribution of 
energy. 

This result is significant in the use of these lamps for solar simulation. Once the 
spectral distribution of the lamp has been carefully determined, for matching absorption 
of energy on the test floor to that in outer space, it will be sufficient to monitor the total 
irradiance of energy on the test floor. On the other hand, the total irradiance should be 
directly measured and cannot be computed as a function of the wattage input of the lamps. 
These two conclusions should, however, be considered tentative for the present, and 
should be tested by more rigorous methods. 

As a check on the internal consistency of the measurements, the data from two ser- 
ies of runs with the two phototubes, and with the same U-21 tungsten lamp, were com- 
pared. The values of log i H from phototube 1P28 were on an average higher than those 
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from phototube 9592B, meaning that the total light output of the lamp was higher in the 
former case. The average of the difference calculated from 83 pairs of values was 0.030. 
Most of these values were taken from the continuum of the spectrum, because the 
sampling of printouts around the peaks did not permit matching of wavelengths. The 
mean deviation of the differences around 0.03 was 0.024. Thus we conclude that during 
the runs with the 1P28 tube the lamp was burning 7.2 percent more brightly (log 1.072 = 
.030), and that the average deviation in the intensity values on a relative scale is 5.7 per- 
cent. Since wavelength assignments for the printouts were highly uncertain, it appears 
that a major part of the 5.7 percent deviation is spurious. Comparison of successive 
runs with the pen recorder shows that on the relative scale the average error is of the 
order of 1 or 2 percent. 

The charts made by the pen recorder showed several interesting features: There is 
a strong continuum whose energy is practically the same for every 100A band throughout 
the spectral range. Superposed on the continuum are the strong emission lines, many of 
which are obviously due to atomic mercury and xenon. A listing of all the spectral lines 
clearly distinguishable from the continuum is given in Table 2-7. Lines whose wave- 
lengths are listed only to four significant figures are unclassified, and the wavelengths 
are estimates from the relative positions of the lines with respect to the known lines of 
mercury and xenon. The unclassified lines are of intensity comparable to the lines of 
xenon. The xenon lines are on the average of intensity about one -tenth that of the five 
groups of strong mercury lines. The intensities of the lines in Table 2-7 are in arbi- 
trary units: , quoted from the "Handbook of Chemistry and Physics," and 1 2 , measured 

approximately from the pen-recorder charts of the mercury-xenon lamp. 

The profiles of the mercury lines are shown in Figure 2-5. All the strong mercury 
lines show clearly the self-reversal effect. The line which originates in the high- 
temperature plasma between the electrodes is very intense. It has also a large halfwidth 
of about 10 or 20A. The profile of the line is highly asymmetric, and is intermediate be- 
tween the Holtsmark profile due to statistical broadening and the Weisskopf profile due 
to collision broadening. There is also a considerable amount of Doppler broadening, 
characteristic of the high temperature of the plasma. Surrounding the plasma is mer- 
cury vapor at a lower temperature which absorbs a part of the energy from the central 
part of the emission line. A dip appears in the middle of the line; the wavelength of the 
lowest point of the dip is very close to that of the sharp emission line of a low-pressure 
discharge. 

A comparison between the results of these measurements and the data reported 
earlier from Westinghouse and Hanovia would be instructive. A direct comparison is 
not possible because the earlier data do not give the spectral distribution within wave- 
length ranges smaller than 100A; the present data, on the other hand, cover a limited 
range (3500 to 6500A) and give the spectral irradiance in one particular direction, and 
of only that part of the source which is focussed on the slit. The earlier data cover the 
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Table 2-7 

Emission Lines in the Spectrum of the Hg-Xe Lamp 
(l x intensity from literature; i 2 , intensity from the 
photoelectric recording. In arbitrary units.) 


k 

Ii 

h 

Element 

4140.38 

200 

10 

Hg 

4193.15 

200 

20 

Xe 

4193.53 

150 

10 

Xe 

4339.23 

150 

400 

Hg 

4347.50 

200 

1000 

Hg 

4358.35 

3000 

6000 

Hg 

4500.98 

500 

150 

Xe 

4524.21 

400 

100 

Xe 

4582.75 

300 

20 

Xe 

4624.28 

1000 

70 

Xe 

4671.23 

2000 

250 

Xe 

4690.97 

100 

10 

Xe 

4697.02 

300 

30 

Xe 

4734.15 

600 

250 

Xe 

4792.62 

150 

10 

Xe 

4807.02 

500 

100 

Xe 

4829.71 

400 

50 

Xe 

4843.29 

300 

20 

Xe 

4916 


1200 


4916.51 

500 


Xe 

4916.04 

50 


Hg 

4960 


700 


5028.28 

200 

40 

Xe 

5047 


30 


5368.07 

100 

150 

Xe II 

5405 


30 


5460.74 

2000 

6000 

Hg 

5769.59 

600 

5000 

Hg 

5789.66 

500 

5000 

Hg 

5790.65 

1000 

5000 

Hg 

5889 


100 


6123 


30 



k 

Ii 

! 2 

Element 

3341.48 

100 

550 

Hg 

3350 


20 


3385.25 

200 

70 

Hg 

3349 


40 


3516 


20 


3542 


300 


3543.08 

40 

50 

Hg 

3559 


250 


3578.75 

40 

30 

Hg 

3604.09 

50 

50 

Hg 

3605.80 

200 

500 

Hg 

3650.15 

500 

3000 

Hg 

3662.88 

400 

1700 

Hg 

3663.28 

400 

1700 

Hg 

3750 


350 


3753 


50 


3769 


100 


3789 


400 


3806.38 

200 

20 

Hg n 

3819 


250 


3858 


120 


3897 


40 


3902 


60 


3908 


100 


3914.29 

100 

20 

Hg 

3918.92 

200 

10 


3983.98 

400 

900 

Hg n 

4008 


20 


4046.56 

200 

3000 

Hg 

4077.81 

150 

2000 

Hg 

4094 


10 


4110 


60 



whole spectral range and give the integrated result of the radiation from the whole source 
in all directions. 

However, a meaningful comparison may be attempted in the following way: The anti- 
logarithm of the values in Table 2-6 gives the intensity i in microwatts per steradian 
per mm 2 of the source for a 10A range. To find the total energy radiated in the spectral 
range covered by the data, multiply I H by 0.4 or 0.5 (depending on whether the wavelength 
intervals are 4A or 5A) and add all the values together. The results from the tube 9592B 
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do not cover the range 3500 to 3600A; therefore, the total for 3600 to 6500A must be con- 
sidered. If the values from the two sets of data (for 1P28 tube and 9592B tube) are aver- 
aged, the lamp is found to emit 1.2295 watts in the spectral range 3600-6500A per mm 2 
of the surface of the source per unit solid angle in a direction normal to the axis of the 
tube. 

Now consider an idealized source, uniformly bright sphere which radiates the same 
amount of energy in all directions. In the actual source, these assumptions are only ap- 
proximately valid, because the irradiance along the axis of the tube is less than that along 
the normal, and the irradiance from the periphery of the plasma is less than that from the 
center. Let r be the radius and A ( = m2 ) the area of cross section of the idealized 
source. The total radiant energy in all directions from the whole source = 1.2295 watts/ 
mm 2 X 477 X A watts. 

From the Westinghouse data, the energy radiated in the range 3500 to 6500A is 412.2 
watts; the Hanovia data for the same range gives the value 534.6 watts. The average of 
the two is 473.4 watts; deducting from this 10 watts (the average of the two sets of data 
for 3500 to 3 600 A), the result is 463.4 watts. 

It seems reasonable to assume that the relative spectral distribution is independent 
of the power level at which the tube is operated. The total radiation in these measure- 
ments could well be made equal to 463.4 watts by a slight adjustment of the power level. 
This produces the equation 

watts „ . 

1. 2295 y * 477 A = 463.4 watts 

Whence A = 29.99 mm 2 or r = 2.9 mm. 

This value is of the right order of magnitude. 

To obtain from the present data a set of values normalized to the same total energy 
as the Hanovia-Westinghouse data, all values of i H should be multiplied by 29.933 x 4 77/ 
10 = 37.69. This would give the energy radiated in microwatts by the whole lamp in all 
directions per spectral range of 1 angstrom. 

This procedure of normalization leads to certain other significant results. 

Of the total energy of 463.4 watts in the range 3600 to 6500A, a major part is con- 
centrated in five short spectral ranges around the strong emission lines of mercury. 
These ranges are shown in Figure 2-5 mentioned earlier. For purposes of comparison 
with the Westinghouse- Hanovia data, a wavelength band of 100A was chosen for each line 
or group of lines, and the total energy emitted in each band was calculated. The results 
are presented in Table 2-8 along with those interpolated from the Westinghouse and Han- 
ovia data. 
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Table 2-8 

Energy Due to the Line Spectrum 


Wavelength Band A 

Strong Hg Lines A 

Spectral Radiance in 

Walls 

Ours 

Westinghouse 

Hanovia 

3600 - 3700 

3650, 3662, 3663 

67.4 

43.8 

46.9 

4012 - 4112 

4046, 4077 

31.8 

37.6 

33.1 

4308 - 4408 

4347, 4358 

54.3 

50.6 

41.2 

5410 - 5510 

5460 

58.2 

44.9 

82.5 

5730 - 5830 

5770, 5790 

105.4 

47.9 

104.7 


The energy in these five wavelength ranges combined is 317.1 watts, about 68 per- 
cent of the total energy in the range 3600 to 6500A. The rest of the energy, 32 percent, 
is distributed over the continuum. Over most of the continuum the energy in each 100A 
range is the same. It varies between 0.8 percent and 1.4 percent of the total energy. 

The average energy per 100A range is 1.1% of the total, that is, 5.1 watts. The spectral 
distribution curve of the solar energy, normalized to the same total output as the 2.5-kw 
mercury-xenon lamp, gives for the region between 4000 and 6500A an average irradiance 
of 23.7 watts per 100A range. Thus it is seen that, if the total radiant energy from the 
mercury-xenon lamp is adjusted to be the same as that due to solar radiation in outer 
space, the energy input on the test floor of a solar simulator is only one-fifth that of the 
solar energy over a large part of the spectral continuum, while it is 25 times greater 
than the solar energy over five narrow wavebands each of about 10 or 20A. This signifi- 
cant difference between the spectral distributions of the solar radiation and the Hg-Xe 
lamp radiation clearly shows the necessity for accurate determination of the spectral 
quality of the radiation incident on the test floor of the simulator. 
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PART III 


HEAT TRANSFER AND CRYOGENIC PUMPING 
IN SPACE SIMULATION 


PRELIMINARY STATEMENT 

(J. H. Boeckel, Associate Chief, Test and Evaluation Division, GSFC) 

While there is a good deal of glamor attached to the space program and our contri- 
bution to it, our work in the end is largely the application and extension of basic engi- 
neering diciplines. I would like to discuss two such areas in which I feel the Summer 
Workshop could make a fruitful contribution. 

First, simulation of the space environment requires that the test object "see" the 
blackness of space while subjected to a very low pressure. To simulate black space, the 
temperature of the chamber walls needs to be below 100°K; Liquid Nitrogen (LN 2 ) pro- 
vides an excellent heat-transfer medium at these temperatures. Once an extended area 
LN 2 surface has been introduced into the chamber, one soon realizes that water vapor, 
the bete noir of high vacuum, will condense on these panels and will in effect be pumped 
Lom the system. Further, with LN 2 available as a thermal shield, it becomes practical 
to operate a gaseous helium (GHe) system at a temperature of about 20°K to condense 
nearly all the remaining gases. 

Cryogenic pumping in this manner has received a good deal of attention. However, 
logical design of a system depends heavily on the actual gases which must be removed. I 
feel that identification of the constituents of the atmosphere within a space simulator is a 
matter of great importance. When this has been accomplished, the study can then be ex- 
tended to cover the manner in which these gases are best pumped. 

Second, simulation of the space environment requires the provision of a source of 
thermal radiation which simulates the input from the sun. We now have what we believe 
to be a good optical design for a solar simulator; however, the engineering practicality 
of the system remains in question. In particular, the energy balance in a system which 
dissipates some 800 kilowatts in a confined area deserves careful analysis. Optical tol- 
erances, vacuum seals, material properties may all depend heavily on how the excess 
energy is channeled. This analysis clearly leads to many other problems in the engi- 
neering design of the solar simulator. 
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Although I have cited two specific questions whose study I believe would be reward- 
ing, we hope for a greater benefit from this program than the generation of two isolated 
pieces of information. In our group we are very close to the turmoil of last-minute op- 
erations in getting a satellite ready for launching. Thus we are somewhat concerned 
that, in our evaluation of our needs, we have in the press of circumstances overlooked 
important unresolved problems. We believe that the Summer Workshop may well uncover 
some of these for us. 

Last, as Dr. Goett said earlier, we feel that contact with the local scientific commun- 
ity can be of great mutual benefit over a long period of time. The Summer Workshop 
gives us an excellent opportunity to get to know each other, our capabilities, and our 
goals, so that both groups can expect quick and effective response in the future when 
either has a problem and turns to the other for help. 
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SECTION A 


CRYOGENIC PUMPING AND VACUUM TECHNOLOGY 

by 

J. L. Shereshefsky, Murray Geller, 

George Pick, and J. B. Crooks 


ABSTRACT 

This report includes a study of the vapor pressures of ice, 
water, solid and liquid C0 2 , liquid Hg, and liquid CO at moder- 
ate and low temperatures. It is shown that Antoine’s equation 
log P = A - B/(C + T) fits the data exceedingly well, and lends it- 
self to more confident extrapolation. 

The report contains an account of a cryosorption experi- 
ment at liquid- nitrogen and helium-vapor temperatures, with 
coconut charcoal as adsorbent. The pressures obtained at 
liquid-nitrogen temperature were 20 times lower than without 
the adsorbent, and about 200 times lower at the helium-vapor 
temperature. 

It also includes an account of a new method of calibration 
of an ionization gauge; a discussion of several problems related 
to cryogenic pumping; and a review and bibliography of adsorp- 
tion of common gases on steel and related metals. 


INTRODUCTION 

Within the broad frame of reference of "heat transfer and cryogenic pumping" which 
was the area of its concern. Team B has chosen to devote itself mainly to study and in- 
vestigation of the following subjects: 

• Vapor pressure and temperature relations of residual gases 

• Cryosorption as a method of cryogenic pumping 

• A method of calibration of ionization gauges in high vacuum 

• A survey of the literature on adsorption of common gases on steel and other metals 
It has also considered various problems related to analysis of residual gases, measure- 
ment of contamination on the surfaces of spacecraft, location of leaks in cryoliners at 
liquid-nitrogen temperature, and orientation of vacuum sensors in a space simulator. 
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VAPOR PRESSURE AND TEMPERATURE 
RELATIONS FOR RESIDUAL GASES 


The purpose of this study was to assess the validity of existing equations for vapor 
pressure at low temperatures of the common gases and vapors occurring in the atmos- 
phere and evolved by the various components of a vacuum chamber. Time limitations did 
not allow an extensive investigation. However, as shown in the following analysis, the 
Antoine equation when applied to liquid and solid H 2 0 and C0 2 , liquid CO, and liquid Hg 
gives an extremely good fit over a wide range of temperatures and can be extrapolated 
with greater confidence than can the equations containing five constants which are pres- 
ently in use. 

Theoretically, one can start from the Clapeyron equation (which is exact) 

dP = ah (1) 

dT TAV 

where 

P = pressure 
T = temperature 

AH = enthalpy of phase change from solid — gas or liquid— gas 
AV = change in volume from solid— gas or liquid— gas , 


or 


and since V con « v gas 


dP _ AH (2) 

dT ~ T (V gas - V con ) con. = condensed phase 


dP _ AH V 1 ( Vcon \ ‘ 

dT - TV gas \ V gas ) 


For the V gas , we can use a virial expansion 


PV 


gas 



RT 


where r = gas constant, = virial coefficient. 
Substituting (4) in (3) for V gas we obtain: 


00 



L 

i ~ n 



(3) 


(4) 


( 5 ) 
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A usual and extremely valid (except near critical point) assumption is to let 



so that the only major stumbling block to integration of Equation (5) is to find an expres- 
sion (analytical or otherwise) for AH. A first approximation is to let AH = constant and 
bj = 0 (i= 1 , . . ., <x) and bo = 1 so that we can integrate to obtain 


InP = - 


AH 

RT 


+ Constant 


( 6 ) 


This is the celebrated Clausius - Clapeyron equation: 
Instead of assuming AH= Constant, an exact expression is 


AH = AHo + f (Cp gas - Cp con ) dT ( 7 ) 

J o 

where 

Cp = heat capacity 

AHo = Constant of integration-enthalpy at absolute zero. 

This produces (within the approximation of V con ./v gas = o) an exact expression forP vs T: 


InP + bjP + b 2 p 2 + 



Thus, an accurate determination of vapor pressures is possible at all temperatures, if 
the heat capacities are known well. 


Using empirical expressions for heat capacities (see any standard thermodynamic 
text) as 


Cp = A + BT + CT 2 (9) 

Cp = A 1 + B X T + CVT 2 

or others (and this is a dangerous procedure for the range of validity is limited), an 
empirical expression for the vapor pressure vs temperature can be obtained in any of 
the forms: 


log P = A - | (10) 

log P = A - | + c log T 

log P = A - | + CT 
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log P = A - | + C log T + DT + FT 2 + 


The disadvantage of most (if not all) of these equations is that they are strictly limited 
to the range in which the empirical constants are evaluated, and cannot be extrapolated 
to lower temperatures (which is the case of interest to us!) 

An empirical equation with the advantage that it can be extrapolated out of the region 
of determination of the constants is the Antoine equation 

1obP = A "cTt (U) 

Rossini* claims that it is the best empirical equation for vapor pressure and is the most 
reliably extrapolated. What follows is an attempt to evaluate the constants for H 2 0, C0 2 , 
and CO, and to compare the calculated pressures with experimental values in other re- 
gions of the scale. 

H 2 0: Liquid-vapor pressure 

log P = 5.12146 -1696,7 (12) 

T -42 


T(°K) 

p (atm) predicted 

p (atm) - exper. 

275 

.006909 

.006883 

300 

.035081 

.034844 

325 

.13361 

.13329 

350 

.40988 

.41069 

375 

1.0623 

1.0675 

400 

2.4102 

2.4232 

425 

4.9137 

4.9337 

450 

9.1804 

9.1968 

475 

15.958 

15.947 

500 

26.113 

26.050 

525 

40.607 

40.490 

550 

60.461 

60.387 

575 

86.724 

87.014 

600 

120.44 

121.88 

625 

162.62 

166.93 


> error 

_ 38 
.68 
.24 
.20 
.49 
.54 
.41 
.18 
.07 
.24 
.29 
.12 
.33 
1.18 
2.58 


.53 average % 
deviation 


(Experimental data NBS Circular 564 (1955)) 


Solid - vapor pressure 

27Q2 5 

1 og P — 7.83158 - L 

T + 5 


( 13 ) 


*Rossini-Chemical Thermodynamics 
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T (°K) 

p (atm) predicted 

p (atm) - exper. 

270 

.0 4754 

.0 2 4633 

260 

.0 2 1967 

.0 2 1933 

250 

.0 3 7596 

.0 3 751 

240 

.0 3 2714 

.0 3 2696 

230 

.0+8884 

.0+883 

220 

.0+2633 

.0+262 

210 

.0 5 697 

.0 5 687 

200 

.0 S 162 

•0 s 158 

190 

.0 6 324 

-0 6 32 

180 

.0 7 546 

.0 7 5 

170 

.0s749 

.0 8 7 

160 

.0+808 

.0+8 

155 

.0+239 

.0+3 


Empirical constants evaluated from 210-260; 


rest of values are extrapolated 
(Experimental data NBS Circular 564 (1955)) 

Using the same equation and changing units to mm 

log P mm = 10.71239 - 2792 5 

t+ 278.1 ( 14 ) 


T (°C) 

P mm predicted 

P exper. 

P. K. Kelley 

-20 

.7815 

.776 

.779 

-30 

.2863 

.2859 

.2931 

-40 

.09641 

.0966 

.0994 

-50 

.02951 

.02955 

.02985 

-60 

.008103 

.00808 

.00840 

-70 

.001965 

.00194 

.00198 

-80 

.000413 

.00040 

.000413 

-90 

.0 4 736 

.0+70 

.0 745 

-100 

.0+108 


.0+110 

-110 

•0 S 126 


.0 S 125 

-120 

.0 6 112 


•0 6 113 

-130 

.0 8 719 


•0 8 698 

-140 

.0 9 3 10 


.0 9 293 

-150 

.0 n 819 


.0i+74 

-183 

.0 lg 223 


•0 2 1 16 


The experimental values are taken from Handbook of Chemistry and Physics (&I.C.T.). 
Experimental values go down to -98°C, and in that range the predicted values agree ex- 
ceedingly well (using Antoine-type equation). The column marked Kelley* is not experi- 
mental but based on the accepted empirical equation 


log P mm = 1.207+3.857 logT +3.41.10' 3 T +4.875.10 -8 T 2 2461 

T 


(15) 


Dushman -Vacuum Tech. 
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derived from Af expressions. Thus, although Kelley uses a five-constants expression, 
the values are less sure than the three-term empirical Antoine equation. Though both 
empirical equations may not be highly valid below -120°C, the Kelley expression yields 
much smaller values than the Antoine equation. 

Preliminary calculations on solid C0 2 seem to indicate the same general results. 
The Antoine equation (16) gives an extremely good fit over a wide range and can be ex- 
trapolated, whereas the Kelley empirical equation (17) 


log P mm = 9.4364 - 


1217.4 
t + 264.1 


(16) 


log Pmm = 8.882 + 0.8702 log T - 3.891. 10' 3 T - 


1408 

T 


(17) 


gives much too small values at low temperatures. 
The results for solid C0 2 are as follows: 


log 


P mm = 9. 4364 - 


1217.4 
T - 9 


(18a) 


T K 

P(pred.) mm 

P(exp.) mm 

% error 

135 

.595 

.59 

_ 

145 

3.054 

3.05 

- 

155 

12.53 

12.53 

- 

165 

42.90 

42.91 

- 

175 

126.67 

126.62 

- 

185 

330.65 

330.8 

- 

195 

778.47 

782.4 

.5 

200 

1155.0 

1165.8 

1.18 

205 

1679.5 

1706.4 

1.58 

210 

2397.1 

2457.9 

2.47 

215 

3362.8 

3491.4 

3.68 


Data from NBS-Circular 564 (1955). 


Agreement is very good until the critical temperature 217 K is approached, 
same equation, but putting it in the form of °C, produces 


log P mm = 9.4364 - 


1217.4 

t + 264.2 


Using the 

(18b) 


Its application to the data of Meyers and Van Dusen is given in the following table. 
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T°C 

p(pred.) mm 

p (exp.) mm 

% error 

P (pred.) Kelley 

-60 

2982.6 

3073.1 

3.03 


-70 

1471.0 

1486.1 

1.03 


-80 

671.87 

672.2 

.05 


-90 

280.48 

279.5 

.36 


-100 

105.26 

104.81 

.43 


-110 

34.791 

34.63 

.46 


-120 

9.862 

9.81 

.53 


-130 

2.317 

2.31 

.30 


-140 

.4310 

.431 



-150 

.05973 

.0605 

1.29 

.06194 

-160 

.0 5664 

.0 2 59 

4.00 

,0 2 730 

-170 

.0 3 326 

.0 3 37 


.0 3 382 

-180 

.0 4 0951 

.0 4 13 



-185 

.0 5 1162 

.0s 17 



-189 

•0 6 18 

.0 6 3 




(Data from Meyers and Van Dusen, J. Res. NBS 10, 381 (1933)) 


Kelley uses the expression 


log P mm = 8.882 + .8702 log T - 3.891.10" 3 T - 1408 

T 

derived from af data, but this equation gives a very poor fit. 

For liquid C0 2 , the best Antoine constant is actually zero, which produces: 


(19) 


log P mm = 7.57368 - 861 7 
T 


( 20 ) 


This equation gives an average deviation of ±.33 percent over the liquid range from 
216°K to 304°K, as shown in the table below. 


t°k 

p mm (pred) 

p mm (exp.) 

% error 

220 

4,538 

4,500 

.84 

230 

6,717 

6,705 

.18 

240 

9,622 

9,625 

.03 

250 

13,393 

13,393 


260 

18,174 

18,143 

.17 

270 

24,110 

24,023 

.36 

280 

31,346 

31,196 

.48 

290 

40,022 

39,865 

.39 

300 

50,274 

50,329 

.11 


(Data from NBS Circular 564 (1955)) 


The data for Hg (liquid) can best be represented by an Antoine equation: 

log P mm. = 7.70536 - ? 975 3 

t + 260 


( 21 ) 
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575.2 

960.86 

557.86 

305.87 

156.87 
74.377 
32.123 
12.405 

4.1813 

1.1923 

.2758 


157- 

961 

55' 

30 

15i 

7 

3 

1 


40 

,0 2 6133 

.0 2 6079 

.90 

20 

,0 2 1200 

.0 2 1201 

- 

0 

.0 3 1828 

•° 3 185 

1.19 

-20 

.0 4 203 

.0 4 181 


-38.9 

.0sl77 

.0 s 124 



Kelley gives an expression for liquid Hg: 


3285 / 

log P mm = 10.377 - .8254 log T - — - — 

This expression gives very poor results, in comparison with the Antoine equation 
which holds extremely well over 9 decades of pressure. 

The results for liquid CO, using the Antoine equation, 


log P mm = 6.74265 


303.57 
T - 3 


are given below: 



T (°K) 

P (pred) mm 

P (exp) mm 

% error 

Triple point 

68.09 

11.99 

11.53 

3.99 

70 

162.85 

158 

3.07 


75 

336.07 

332 

1.23 


80 

631.26 

629 

.36 

Boiling point 

81.62 

761.1 

760 

.15 

85 

1,098.0 

1,098 

- 


90 

1,792 

1,796 

.22 


95 

2,773 

2,780 

.25 


100 

4,103 

4,110 

.17 


105 

5,841 

5,840 

.02 


110 

8,046 

8,040 

.08 


115 

10,770 

10,770 

- 


120 

14,062 

14,110 

.34 


125 

17,963 

18,140 

.97 


130 

22,509 

22,990 

2.14 

Critical point 

132.88 

25,430 

26,242 

3.10 
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Summary of constants in the Antoine equation for the several substances investigated: 


log P mm = 


A - 


B 

T + C 


(23a) 


System 

A 

B 

C 

H 2 0 - liquid 

5.12146 

1696.7 

-42 

H 2 0 - ice 

7.83158 

2792.5 

+5 

C0 2 - liquid 

7.57368 

861.7 

0 

C0 2 - solid 

9.4364 

1217.4 

-9 

Hg - liquid 

7.70536 

2975.3 

-13.2 

CO - liquid 

6.74265 

303.57 

- 3 


3-11 




SUMMER WORKSHOP REPORT 


SECTION B 


CRYOSORPTION AS A METHOD OF CRYOGENIC PUMPING 

INTRODUCTION 

The cryopanels presently used in vacuum chambers and simulators are of limited 
surface area, adequate only for condensing the vapors in the residual gases. For con- 
densing the noncondensible gases, the chambers must be equipped with active adsorbents 
having large specific surfaces. At low temperatues, these materials also tend to lower 
the vapor pressure of the condensible gases. The condensation process, whether of va- 
pors by lowering the temperature, or of gases by adsorption, is nearly instantaneous 
under proper conditions and is in any case much faster than a diffusion pump. These 
materials can be easily incorporated in existing cryopanels, or cryopanels with surface 
properties of active adsorbents can be developed. 

To demonstrate the applicability and effectiveness of cryosorption in cryogenic pump- 
ing, an inner shroud was designed and built into an existing vacuum chamber, and an ex- 
periment carried out as described below. It should be pointed out that the design and 
form of the shroud are far from ideal, and were conditioned by considerations of time, 
availability of material, and ready fabrication facilities. 


DESCRIPTION OF EXPERIMENT 

EXPERIMENTAL FACILITY. The experimental apparatus is shown in Figures 3-1 
and 3-2. Testing was conducted in a stainless steel cylindrical vacuum chamber 22 
inches in diameter, and 30 inches in height. It was equipped with an oil-vapor diffusion 
pump, type PMC-4100, backed by a mechanical pump, and included an aluminum cryo- 
panel. The oil diffusion pump was equipped with a liquid nitrogen- cooled trap to prevent 
the oil vapors from diffusing into the test chamber. 

The cryogenic fluid used in the cryopanel was liquid nitrogen, which entered under 
3-4 lb/in 2 pressure from a pressurized container. The liquid nitrogen evaporated in the 
cryopanel and the vapor was vented out. No attempt was made to recover the nitrogen 
vapor. 

TEST ARRANGEMENT. The experimental shroud assembly, mounted within the 
aluminum cryopanel as shown in Figure 3-3, was made of copper and consisted of 18 
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Figure 3-2 - Experimental apparatus, stainless 
steel vacuum test chamber 


3-14 



HEAT-SINKING TECHNIQUES 



Figure 3-3 - Assembly of experimental shroud, diagram 

vertical 3/8" O.D. tubes. The tubes formed five sections connected at the bottom of the 
assembly, approximating a cylindrical array. Total length of the tubes was approxi- 
mately 35 feet. 

Each tube was equipped with 29 cone-shaped cups which were silver-soldered at 
equal intervals along the length of the tube surface. Total number of cups was 532. The 
cups (13/16-inch deep, and 2 inches base diameter) were fabricated from 0.020-inch 
thick copper sheet material, and were formed on the lathe, using a die. Figure 3-4 shows 
the experimental shroud before installation, and Figure 3-5 is a detailed drawing of one 
section. 

Deposited in the cups was a total amount of 2300 gm of activated coconut charcoal 
having a potential specific surface of 890 m 2 /gm. To obtain low surface temperatures, 
liquid nitrogen was introduced in the tubes as a coolant (77°K); helium vapor was also 
used, providing an even lower temperature (~25°K). No attempt was made to recover the 
nitrogen or helium vapors. 


INSTRUMENTATION. Pressure and temperature sensors were installed to provide 
information about the pressure and temperature distribution. Five pressure gauges pro- 
vided pressure data in different locations. One thermocouple gauge measured the fore- 
pressure in the mechanical pump, and a similar gauge registered the pressure in the 
diffusion pump. In the test chamber, one thermocouple gauge registered the pressure 
from 1 mm to 1 x 10 3 mm/Hg. Below this pressure range, a Veeco tubular ionization 
gauge (in connection with a Veeco type RG-31A control panel) registered the pressure in 
the vicinity of the vacuum-chamber wall. In the middle of the test chamber, a nude gauge 
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Figure 3-4 - Experimental shroud 
before installation 


INNER SHROUD FOR CRYOGENIC PUMPING 


SIU/EF SOLDER 


SINGLE UNIT WITH 
&Tu™ PLANE 


EXHAUST TUBE — ' c 
WT«£ FEEDER TUBE 


SEAL WITH SILVER SOLDER 


'SILVER SOLDER 


Figure 3-5 - Section of inner shroud, detailed drawing 
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was installed in connection with a Cooke ionization-gauge control panel. Both gauges 
were capable of measuring pressures between 10 ' 4 to 10 ~ 9 mm/Hg. 

The temperature distribution of the experimental shroud and of the cryopanel was 
observed by using 15 copper-constantan thermocouples (30 gauge) in connection with a 
switchboard, an ice junction, and a Minneapolis -Honeywell Model No. 2733 potentiometer. 
The measuring range of the potentiometer was between 0.005 to 100 mv, with the smallest 
division being 0.005 mv. 

The thermocouples registering surface temperatures were fastened to the surfaces 
by means of a special Scotch brand aluminum tape with adhesive on one side.* The ther- 
mocouple wires were protected by radiation shields. The thermocouple location chart is 
shown in Figure 3-6. 


Symbols: • center pipe 
if cup 

No. 13 Outside of outer shroud 
No. 11 Inside of outer shroud 
No. 7 Outside of outer shroud 

(Opposite to No. 3 shroud panel) 



Figure 3-6 - Thermocouple location chart 


EXPERIMENTAL PROCEDURE AND RESULTS 

The pump-down procedure was designed to minimize oil contamination and build-up 
of condensate upon the surfaces of the test assembly. Immediately after the vacuum 
reached 2.5 X 10 ' 2 mm/Hg pressure, obtained with the mechanical pump, the cooldown of 
the liquid nitrogen trap was begun and continued concurrently with the warm-up of the dif- 
fusion pump. When the diffusion pump became effective, chamber pressure dropped to 
7 X 10 -s mm/Hg, and the bakeout of the coconut charcoal was begun. 

The baking of the charcoal was done by means of hot dry nitrogen gas. The gas, 
pressurized at about 2 psi, was forced through an electrically heated coil where it 
warmed up to 220°C. The hot gas was then introduced into the inner shroud to heat the 
charcoal. 


The adhesive material did not endure the severe thermal cyclings sufficiently, and three thermocouples were lost in the 
course of the experiment; however, the tape material did not show evidence of extensive outgassing. 
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Table 3-1 contains selected parameters related to the baking period. 


Table 3-1 


Temperature drop of nitrogen gas 

70 °C 

Temperature drop between gas and tube-wall surface 

20° C 

Average temperature drop in a vertical tube 

47°C 

Average temperature difference between tube and cup surface 

3.5°C 

Temperature drop between the shroud and the cryopanel 

53° C 

Temperature uniformity at the same elevation 

4% 

Pressure ratio between inner and outer gauge 

3.4 


Figure 3-7 contains the pressure- and temperature-time history of the bake-out per- 
iod. The bake-out period lasted 4 hours, of which the first 35 minutes took to reach a 
temperature of 98°C. After this, the temperature dropped slightly. Parallel with this 
process, the pressure increased sharply in the first hour and then dropped very sharply 
in the next hour. The gas flow was adjusted and a higher gas temperature level was ob- 
tained; concurrently, the pressure increased again and reached in the next hour a peak 
which was very much smaller than the first peak:. The pressure dropped rapidly after- 
ward. It is interesting to note the parallel behavior of the pressure-time curve of the 
inner and outer gauges; the general behavior of the curves is largely the same, but the 
ordinate differs by a nearly constant factor. 

After the baking was completed, 30 minutes was required to cool down the assembly 
to room temperature. Dry cold nitrogen gas was used for this purpose. 

Upon completion of the pumping, liquid nitrogen was introduced into the shroud, and 
in one -half hour the pressure in the middle of the vacuum chamber dropped 2 decades; 
in the next hour it dropped nearly another decade. The difference of pressure readings 
between the inner and the outer gauges increased roughly to a decade. 

On reaching a steady state after 8 hours, the cryopanel was cooled with liquid nitro- 
gen and the inner shroud was warmed with gaseous nitrogen to clear it of liquid nitrogen. 
After this, liquid helium was introduced at a pressure of about 4 ounces and with an 
approximate rate of 1/2 liter/minute. In 10 minutes, the pressure in the inner portion of 
the chamber dropped 2 decades to 3.2 X 10‘ 9 mm/Hg. In the next 10 minutes, the pres- 
sure dropped by a factor of 2 to its lowest value. Doubling the liquid helium rate did not 
affect the pressure. 

The pressure registered by the outer gauge after reaching its lowest value remained 
approximately 2 decades above the pressure registered by the inner gauge. 

It should be noted that, in the middle of the experiment, a change of helium tanks was 
made which took 5 minutes; during this period, the inner pressure gauge registered a 
three-decade increase in pressure. This phenomenon shows that the pressure is ex- 
tremely sensitive to temperature variations, and that the response is fast. 
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Figure 3-7 - Pressure- and temperature-time history of bakeout period 

Figure 3-8 is the pressure-time history of the cooling period of the experiment, and 
Figure 3-9 is the temperature-time history. Table 3-2 contains selected data pertinent 
to the cooling period with liquid nitrogen used as coolant; Table 3-3 contains selected in- 
formation pertinent to the cooling period when liquid helium was used as coolant. 
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2 4 6 8 10 12 14 16 

TIME (hours) 

Figure 3-9 - Temperature -ti me history 
Table 3-2 

1 

18 

Average temperature of the vertical tube 

-175°C 

Average temperature rise in vertical tube 

23°C 

Maximum temperature rise between inlet and outlet 

30°C 

Temperature rise between the LN 2 and the tube wall 

3-5°C 

Temperature difference between the inner shroud and the cryopanel 

143°C 

Average temperature difference between the cup surface and the tube 

21°C 

Temperature uniformity of the same elevation 

2% 

Temperature fluctuations in the same location 

1% 

Corr. pressure ratio between the inner and outer gauge 

10.6 


Table 3-3 


Average temperature in the vertical tube 

-253°C 

Average temperature rise in a verticle tube 

9°C 

Maximum temperature rise between inlet and outlet 

22°C 

Temperature difference between the cup surface and tube 

20° C 

Temperature uniformity at the same elevation 

8% 

Temperature fluctuations in the same location 

1% 

Corr. pressure ratio between the inner and outer gauge 

85 

Temperature of cryopanel 

-186°C 


INTERPRETATION OF RESULTS 

THE ANNULAR AND CENTRAL CYLINDERS. Assuming that the ionization gauges 
reflect the pressures prevailing in the portions of the vacuum chamber with which they 
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directly communicate-and there are no valid reasons for assuming otherwise-the strik- 
ing facts of this experiment are: in the liquid-nitrogen experiment, the corrected steady- 
state pressure registered by the outer gauge is 10 times higher than the pressure regis- 
tered by the inner gauge; in the helium-vapor experiment, the disparity between the two 
gauges is even greater, namely 85 times. 

This is readily explained by the following facts: a) at low surface coverage the ad- 
sorption process is almost instantaneous and the equilibrium pressure is practically con- 
stant; b) the transport of matter from a region of higher pressure to one of lower pres- 
sure by molecular flow is much slower than by mass flow, and is equivalent to a restrictive 
connection between the two regions; c) evolution of gas from the chamber walls and in- 
leakage through O-rings represent a source of gas supply in the annular cylinder of the 
vacuum chamber. 

It thus seems that the vacuum system consists of an annular cylinder delineated by 
the chamber walls and the cryopanel, characterized by a high pressure, and an inner cyl- 
inder delineated by the inner walls of the cryopanel, including the inner shroud with the 
adsorbent. The pressure in the latter is determined by the "vapor pressure" of the ad- 
sorbed gases on the charcoal at the prevailing temperature, while the pressure in the 
former is determined by the efficiency of the pumping system. 

TRANSPORT OF MATTER FROM ANNULAR CYLINDER TO INNER CYLINDER. 
The steady state observed in this experiment cannot be explained on the basis of equili- 
brium established by transpiration, because of the temperature difference between the 
two compartments; such an equilibrium would require that the ratio of the temperatures 
of the two sections of the chamber be about 110 for the liquid -nitrogen experiment, to 
about 6400 in the helium-vapor experiment. The steady state must, therefore, be ex- 
plained on the basis of a process involving mass transport, due to net molecular flow 
from the high-pressure to the low-pressure compartments, and immediate uptake of the 
transported gas by the adsorbent. As will be shown, because of the large area of adsorb- 
ent, the surface coverage remains very low, even after hours of gas flow. The low cov- 
erage accounts for the constancy of the pressure in the inner cylinder. 

The area through which the mass transport takes place is confined to two circular 
openings in the cryopanel, 3 inches in diameter and at right angles to each other, and an 
area between the cryopanel and the base plate described approximately by the surface of 
a cylinder 10 inches in diameter and 1.5 inches in height, the total being approximately 
400 cm 2 . The opening on top of the cryopanel is neglected. 

The rate of transport can be calculated from the equation for the rate of molecular 
collisions with a surface 


Z 



( 1 ) 
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where A is the area,n is the molecular density and c is the average molecular velocity. 
Substituting for n and c the equivalent kinetic terms, z is given by 


z A_ P_ 

(27rmk ' T' A ’ 


( 2 ) 


where P is the pressure, T is the absolute temperature, m is the molecular mass, and k 
is the Boltzmann constant. This equation is also applicable to the evaluation of the back- 
ward molecular flow from the low- to the high-pressure compartments, the net flow being 
the difference between the flows in the opposite directions. Thus the net flow in the di- 
rection of low pressure is given by 


Az 


A /Pi P 2 \ 

(277mk l T 2 !4 j ' 


(3) 


where Pj and Tj are the pressure and temperature of the annular cylinder, and p 2 and t 2 
are the pressure and temperature of the inner cylinder. 

As shown in Table 3-4, the net rate of molecular flow is approximately 9 x 10 16 
molecules, or 3.36 X 10 3 cc of gas at STP, per second, in the liquid-nitrogen experiment, 
and 3.44 x 10 16 molecules, or 1.29 x 10* 3 cc at STP per second, in the helium -vapor ex- 
periment. The backward molecular flow is appreciable, about 15 percent in the latter 
case and about 10 percent in the former case. The composition of the gas is assumed 
to be equivalent to that of air. 


Table 3-4 

Parameters of the Experiment Involved in the Calculations 


Volume of chamber 

197 liters 

Volume of cryopanel 

134 liters 

Volume of annular cylinder 

63 liters 

Surface area of chamber walls 

1.83 sq. m. 

Transport area 

396 sq. cm. 

Specific surface of charcoal (assumed) 

400 sq. m./gm 

Temperature of charcoal in LN 2 experiment 

120 °K 

Temperature of charcoal in He-vapor experiment 

40°K 

Temperature of cryopanel in LN 2 experiment 

241°K 

Temperature of cryopanel in He-vapor experiment 

9l°K 

Time average pressure in annular cylinder in LN 2 experiment 

5.95 X 10" 7 Torr 

Pressure in annular cylinder in He-vapor experiment 

1.5 x 10- 7 Torr 

Temperature of tubular gauge 

350°K 

Rate of transport of gas in LN 2 experiment 

8.97 X 10 16 molec 
sec 

Rate of transport of gas in He-vapor experiment 

3.44 X 10 16 molec 
sec 
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The pressure values for the annular cylinder used in these calculations were cor- 
rected and time-averaged, the corrections being made on the basis of transpiration 
equilibrium between the tubular gauge at an estimated temperature of 77°C and the meas- 
ured temperature of the cryopanel facing the gauge. This corrected pressure reflects 
the pressure of gases at the cryopanel. The equation used is given by 

Pj =p 2 (T 1 /t 2 >k, - (4) 

where Pi is the corrected pressure, p 2 is the pressure registered by the tubular gauge, 

T is the temperature of the cryopanel, and T 2 is the temperature of the glass envelope 
of the gauge. 

FRACTION OF ADSORBENT SURFACE COVERED WITH ADSORBED GAS 
MOLECULES. The calculations of coverage are based on an assumed value of specific 
surface area for the coconut charcoal, of 400 m 2 per gram. This value is about 45 per- 
cent of the value given for the material, when activated at about 350 C. Also involved is 
the value of 13.47 sq.A for the area of an "air-molecule," obtained from the given values 
of nitrogen and oxygen at -253°C. 

The coverage, that is, the fraction of surface covered, after nine hours (the duration 
of the liquid-nitrogen phase of the experiment) is 4.3 X 10 ' 4 ; at the end of the helium- 
vapor phase of the experiment, (after 10 hours and 45 minutes) the coverage is about 
4.6 x 10 ' 4 . It is of interest to point out that, at the average rate of coverage of 4 x 10' 5 
per hour, it would take about 2.5 x 10 4 hours, or 2 years and 9 months, to form a mono- 
layer. This leads to the possibility of either decreasing the quantity of adsorbent, or 
lowering the temperature of activation, which is preferable, or both. It should be pointed 
out, however, that the equilibrium pressure tends to rise with increase of coverage. A 
desirable limiting coverage can be obtained from a study of adsorption vs equilibrium 
pressure. As pointed out elsewhere in this report, the adsorption literature seems not 
to include information of this type in high -vacuum systems. 

SOURCES OF ADSORBED GASES. According to the calculations in section 2, at 
least 3.36 microliters of gas, or 9 x 10 16 molecules, had to be supplied from somewhere 
to the annular cylinder of the chamber each second. If the walls of the chamber (which 
are at a temperature not much below that of the room) are the source of this gas, it 
seems that 5 x 10 12 molecules, or approximately 2 X 10 4 microliters/ cm 2 , would be 
evolved each second, and a monolayer of adsorbed gas would be depleted in 2.5 minutes. 
Although many monolayers of gas would be supplied by the chamber wall, which was not 
baked out, it is quite certain, in view of the duration of the experiment, that the major 
source of the gas supply is due to in-leakage. 
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CONCLUSIONS AND RECOMMENDATIONS 

This experiment and the interpretation of its results show that cryosorption (a) is a 
highly efficient method of cryopumping; (b) is at least 8 times as fast as the 10-inch dif- 
fusion pump used in the experiment; (c) has a high potential for coping with inleakage; and 
(d) produces higher vacua than the diffusion pump capability, 20 times in LN 2 experiment 
and about 200 times in the He experiment. 

It is recommended that further studies be made with respect to: 

a. C ryosorption efficiency of different adsorbents 

b. Cryosorption efficiency of given adsorbents, with respect to different gases 
and to binary gas mixtures 

c. Cryosorption efficiency of adsorbent activated at different temperatures 

d. Design and fabrication of cryopanels that would incorporate efficient 
cryosorbing surfaces 
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SECTION C 


A METHOD OF CALIBRATION OF IONIZATION GAUGES 

INTRODUCTION 

The calibration of ionization gauges is frequently made by establishing the value of 
the constant, K, in equation 


P = K 



( 1 ) 


by comparing the ratio of the ion current i + to the electron current I. of the ionization 
gauge with the pressure reading, p, of a McLeod gauge in the pressure range of 10 4 
to 10 - 5 mm/Hg. This procedure assumes the constancy of K in the higher decades of 
vacuum. However, it has been shown by Nottingham and Torneyi that in the regions of 
pressure 10 - 3 to 10 4 mm Hg, that the constant, K , depends on the pressure and the elec- 
tron current l_. To test the constancy of K at pressures less than 10 5 with the McLeod 
gauge is not reliable, for reasons connected with the reproducibility of McLeod gauge 
readings. 

THEORY OF NEW METHOD 

Since the McLeod gauge is not precise at pressures less than lx 10 5 mm/Hg, and 
the ionization gauges are most useful at these very pressures, it was considered feasible 
to compare the ratio i + /L of an ionization gauge with the pressure calculated from a 
known quantity of gas in a known volume at constant temperature. 

This method based on Boyle's law could be made very precise, if it were not compli- 
cated by the phenomenon of adsorption and desorption. With the expansion of the gas into 
a baked-out chamber which is provided with an ionization gauge, at pressure less than 
10‘ 6 mm/Hg, the amount of gas adsorbed on the surface of the chamber can be an appre- 
ciable fraction of the total. 

To account for the gas removed from the gas phase by adsorption, a small vessel 
containing chamber-wall material in a finely divided state and of known geometric surface 
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area may be brought into contact with the gas in the expansion chamber. From the ratio 
of the readings of the ionization gauge, the adsorption coefficient may be calculated. 

Let Pl be the pressure of the gas confined in the small vessel of volume Vl at the 
constant temperature T, and let V be the volume and S the geometric surface of the cham- 
ber at temperature t . The quantity of gas in moles, n , is given by 

n = PjV^/RT (2) 

This quantity of gas on expanding into the chamber of volume V is given by 


n 


PV 

RT 


+ A, 


(3) 


where P is the pressure and A is the number of moles of gas adsorbed on surface S. On 
expanding the surface to S' , the gas quantity is now given by 


where p' is the gas pressure and a' is the quantity of gas adsorbed. If A' differs from 
A, then p' is different from p. 

At low pressure and in a narrow region of pressure, it is reasonable to assume (on 
the basis of the Langmuir adsorption isotherm) that 


A = KSP 


(5a) 


and 


a' = ks'p' (5b) 

where s and s' are the respective geometric surfaces, P and P' are the gas pressures, 
and K is the adsorption coefficient per unit surface and unit pressure. 

Equating (2) and (3) and combining with (5a) and (5b), we obtain 

^ (P - P') = K(S'P' -SP) (6a) 

or 

5 (*-£)-* («■*-») 
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Now, assuming that in a narrow region of pressure the constant K in equation (1) is 
independent of pressure, therefore at a given emission current i_ , the ratio P'/P is 
given by the ratio of ion corresponding currents of the ionization gauge 



I + 


Combining (7) with (6b), and substituting r for s'/S, produces 


(7) 


V 

RT 




( 8 ) 


From this equation, the value of K , the adsorption coefficient, can be evaluated, when 
it is combined with equations (3) and (5a). In order to compute the gas pressure Pj, 
equation (8) may be combined with equations (4) and (5b). 


APPARATUS 

The apparatus made of pyrex glass, shown in Figures 3-10 and 3-11, includes three 
pressure gauges, a mercury manometer for measuring pressures of 20 mm/Hg and 
above, a small McLeod gauge for pressures of 10' 4 to 1 mm/Hg, and a large Consoli- 
dated Vacuum Corporation McLeod gauge for pressures of 10' 6 to 10 " 4 mm/Hg. It also 
includes a calibrated capillary for introducing desired quantities of gas, and a flask of 
calibrated volume containing a known number of glass spheres of known diameter. The 
flask, prior to the introduction of the glass spheres, was used for the calibration of the 
volume of the system. Two ionization gauges, a Veeco and a NRC, were sealed to the 
apparatus for calibration. The apparatus was evacuated by means of a mechanical pump 
and a Varian ion pump. The ion pump and the ionization gauges were protected from Hg 
vapor by liquid-nitrogen traps. 

The parameters of the apparatus are given in Table 3-5. 
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Figure 3-10 - Ionization-gauge calibration apparatus, diagram 



Figure 3-11 - Ionization-gauge calibration apparatus 
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Table 3-5 


Capillary diameter 

~ 3 mm 

Capillary volume 

0.1187 cc. 

Flask volume 

116.54 cc. 

Diameter of capillary in small McLeod Gauge 

2.635 mm. 

Diameter of glass sphere, average 

2.962 mm. 

Weight of spheres in flask 

131.25 gm. 

Number of spheres 

3806. 

Surface area of spheres 

1057 sq. cm. 

Surface area of apparatus 

1594 sq. cm. 

Volume of apparatus (exclusive of ion pump) 

3.473 liters 


PROCEDURES 

CALIBRATION OF VOLUME OF APPARATUS. Helium is introduced into the appa- 
ratus through the intake stopcock to a desired pressure, read on the manometer; the 
stopcock to the flask is then closed, and the apparatus evacuated to at least 10" 5 mm/Hg. 
The gas is allowed to expand into the apparatus (exclusive of the ion pump) and the pres- 
sure is read again, either on the manometer or on the small McLeod gauge. The volume 
is calculated from the known volume of the flask and the pressure readings, using Boyle T s 
law. 


CALIBRATION OF ION GAUGE. The ionization gauges and the glass spheres are 
baked (the former at about 400°C, and the latter at about 27 5° C) for six to ten hours, 
while the ion pump is evacuating the system to about 2 X 10 8 Torr as indicated by the 
ion vac pump gauge, or until a sticking vacuum is reached on the CVC McLeod gauge. 

The connecting tubes and the McLeod gauges are flamed while the baking and pumping 
is in process. The ionization gauge to be calibrated is allowed to degas and stabilize to 
a constant reading. 

The gas taken for calibration purposes and stored prior to the baking process in the 
calibrated capillary at a known pressure (from 0.01 to 0.2 mm/Hg) is allowed to expand 
into the apparatus, exclusive of the flask of spheres. Readings at different emission set- 
tings are taken on the ionization gauge, and also on the CVC McLeod gauge, if the pres- 
sure is within the reliable range. 

The stopcock to the flask of spheres is now opened, and readings are taken periodi- 
cally on the ionization gauge at different settings of emission current, and if possible on 
the McLeod gauge, until constancy of reading is obtained. 


RESULTS OF DISCUSSION 

Regretably, the time remaining after completing construction of apparatus, and elimina- 
ting unforeseen difficulties, was inadequate for a full test of the method and its limitations. 
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Only one experiment was carried out, on which it is not possible to base any 
evaluation. 


The data are as follows: 

Pressure of nitrogen in capillary: 12.6 x 10" 2 mm Hg (1) 

PV of nitrogen: 14.95 X 10 3 cc - mm Hg (2) 


Veeco reading prior to introduction of gas: 

emission, ma 10 8 6 4 2 

Reading, Torr X 10 7 10.5 8.2 5.4 3.4 1.6 

Veeco reading after introduction of nitrogen: 

emission, ma 10 8 6 4 2 

Reading, Torr x 10 6 2.5 2.0 1.5 1.0 0.57 

Veeco reading after opening stopcock to spheres: 

emission, ma 10 8 6 4 2 

Reading, Torr x 10 6 2.5 2.0 1.5 1.0 0.62 


(3) 

(4) 

(5) 


The background pressure recorded by the Veeco gauge and given under (3) is too 
high, and is evidently due to insufficient degassing of the apparatus. The difficulty of 
baking the several sections of the apparatus simultaneously is the weakest aspect of the 
apparatus design; the sections of the apparatus which require baking should be brought 
together to make possible enclosing it by an oven. 


The reading in (4) and (5) show that no adsorption by the walls has taken place. The 
pressure calculated from Boyle ! s law, 4.3 X 10 6 mm/Hg, is higher than the values indi- 
cated by the Veeco gauge at the several emissions. The low reading on the Veeco is in 
agreement with observations of Barnes^ ascribed to the high ratio of positive ions to 
neutral molecules in the space between the grid and the tube walls. This coincidence 
could be fortuitous, and no affirmative statement should be based on a single observation. 
However, the method of approach shows as possessing potentialities. 
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SECTION D 


A SURVEY OF THE LITERATURE 
ON ADSORPTION OF COMMON GASES 
ON STEEL AND OTHER METALS 


INTRODUCTION 

Originally, in collecting these references, information was sought on the adsorption 
of gases at ultra-high vacuum, and at room temperatures or below, on steel and other 
structural metal surfaces. As satisfactory information of this type was not available, the 
scope of the search has been broadened. References in the appended bibliography will 
be discussed in relation to the following topics: 

The magnitude of adsorption of common gases on steel and on metals going into 
steel; the nature of the "irreversible" adsorption observed on metals; experiments on the 
adsorption of gases at very low pressures. 


ADSORPTION OF COMMON GASES ON STEEL AND SOME OF ITS 
COMPONENT ELEMENTS 

Some direct measurements have been made of adsorption of gases on steel surfaces. 
M. H. Armbruster has contributed several articles on this subject (31-34). Measure- 
ments were made at -195 to 20°C at pressures ranging from 10 3 to 1 Torr, on both re- 
duced and unreduced surfaces. Two types of adsorption were observed: a reversible 
adsorption which could be removed by pumping in a moderate (10 6 Torr) vacuum, and 
an ’'irreversible” adsorption which could not be thus removed without raising the temp- 
erature. The reversible adsorption was recognized to be held by physical (e.g., van der 
Waals) forces, and approached a limit at 0.1 to 1 Torr, comparable in magnitude to a 
monolayer. 

Nitrogen and argon on reduced surfaces at -195°C, carbon dioxide at -78°C, and ethyl 
iodide at 20°C, were found to be held entirely by physical forces. Except when observed 
on unreduced surfaces, the "irreversible” adsorption was thought to be held by chemical 
forces (i.e., chemisorbed). Carbon monoxide, water vapor, argon, and nitrogen on un- 
reduced surfaces were all observed to have at least partial irreversibility of about a 
monolayer or less. Oxygen on a reduced surface was found to be chemisorbed to the ex- 
tent of at least 8 monolayers. It should be noted that temperatures used in reduction and 
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outgassing were no greater than 600°C, which may not have been great enough to produce 
clean surfaces. 

Chessick et al (27) made a similar study on the physical adsorption of common gases 
at pressures greater than 10 5 Torr on molybdenum powder. Kavtaradze (20) has recent - 
cently made adsorption measurements of H 2 , CO, 0 2 , and C 2 H 4 on metal films at pres- 
sures of 10' 7 to 10‘ 2 Torr at temperatures ranging from -195 to 250°C; he found adsorp- 
tion to be either large and irreversible, or small and reversible. This is in agreement 
with what one would expect from Armbruster 1 s data. 

Several articles have been published on the adsorption of oxygen on iron and steel, 
especially in relation to passivity. One such article is that by Burshtein and Levina (77) 
on the adsorption of highly alloyed stainless steels. At 20°C, the magnitude is l/5th the 
amount observed by Armbruster (surface measured by B.E.T. method), i.e., about 1.3 
monolayers. Zettlemoyer and Chessick (76) found that water vapor was completely 
reversibly adsorbed on nickel and steel powders reduced at 25°C, indicating that the oxide 
layer was thick enough to prevent diffusion of water vapor to the metal surface. When 
elevated temperatures were used to "activate" the surface, 10 to 13 layers of water were 
adsorbed at room temperature. This is far in excess of the adsorption observed by 
Armbruster. 

If temperatures are sufficiently high, dissociation will occur on adsorption. While 
Armbruster found only physical adsorption of C0 2 at -78° C, Takao Kwan (81) found evi- 
dence that C0 2 was dissociated into CO and O at 150°-200°C on reduced nickel. Wright 
and Ashmore (43) found that both ethane and methane were dissociated on adsorption by 
evaporated tungsten surfaces at 0-100°C. Ethane was found to adsorb in the same manner 
on nickel and on iron at temperatures greater than 70°C. 

The adsorption on evaporated metal films seem qualitatively different from adsorp- 
tion on powders or bulk metal surfaces; in particular, chemisorption appears to be more 
prevalent and pronounced on the evaporated films. A. S. Porter and F. C. Tomkins (28) 
measured the adsorption of CO and H 2 at 10 5 Torr on iron films. Most of the gas ad- 
sorbed at high pressures was found to remain to the low-pressure limit. The heat of ad- 
sorption for both carbon monoxide and hydrogen was calculated to be about 20 kcal/mole. 
Armbruster, in contrast, did not observe any adsorption of H 2 , and calculated the average 
heat of adsorption for CO to be about 3.5 kcal/mole. While Armbruster found only physi- 
cal adsorption for C0 2 on reduced steel at -78°C, Collins and Trapnell (65) found strong 
chemisorption of C0 2 on Fe, W, Mo, and Ni films at -78°C to 70°C. Mazumdar (50) found 
that ethane and methane were chemisorbed on Ni films at -195°C. 

IRREVERSIBLE ADSORPTION AT LOW PRESSURE 

No reference has been found which gives information on the evolution of irreversible 
adsorbed gases at pressures below 10" 7 Torr. The data of Porter and Tomkins (28) 
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suggest that this evolution could be appreciable for CO and H 2 on evaporated iron films. 
Also, Armbruster thought that the irreversible adsorption of N 2 and A on unreduced sur- 
faces of steel at -183°C might have been due to strong physical forces between the ad- 
sorbed molecules and an oxide layer. The reasons put forth for suspecting physical ad- 
sorption were (1) no adsorption for N 2 and A was observed at -78°C, indicating the weak 
nature of the forces involved; (2) argon, usually chemically inert, was irreversibly ad- 
sorbed to about the same extent as nitrogen. One might expect adsorption of this charac- 
ter to desorb at outgassing pressures lower than those used by Armbruster (10 6 Torr). 

The heat of adsorption is one indication of the strength of the bonds holding an ad- 
sorbed substance. Heats of physical adsorption are usually less than 10 keal/mole; heats 
of chemisorption, on the other hand, may be as great as 100 kcal/mole. Molecules ad- 
sorbed in a state characteristic of such high heat (adsorbed atoms, for example) have a 
negligible equilibrium pressure at room temperature or below; however, materials with 
intermediate heats of adsorption might evolve at a rate appreciable at very low pres- 
sures. Porter and Tomkin f s data suggested this might be 20 kcal/mole. 

SOME EXPERIMENTS ON ADSORPTION ON METALS IN HIGH VACUUM 

Three references have been found which deal with the adsorption of gases at pres- 
sures of about 10' 8 Torr. Ehrlich and Hudda (22) made a study of the adsorption of the 
rare gases on a tungsten field emitter at low temperatures, and found no detectable ad- 
sorption of He, Ne, or A below 10" 7 Torr above 20°K. R. J. Hill (21) made a field- 
emission study of the adsorption of water vapor on tungsten at 80°K and 400°K; but the 
adsorption is only indicated qualitatively by photographs of the resulting patterns. 

P. Kisliuk (51) has determined the adsorption of nitrogen on tungsten to be from 10 8 to 
10- 5 Torr at 1100° to 1700°K. As might be expected, the nitrogen is dissociated into 
atoms, and adsorption is proportional to the square root of pressure. 

A complete search of the literature on this topic is expected to yield more informa- 
tion of this nature. However, it seems that the literature is lacking in systematic equili- 
brium studies of adsorption in the pressure region below 10 6 Torr. 
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APPENDIX A 

INTERCHANGE OF INFORMATION 

From the interchange of information and ideas between 
members of the Summer Workshop and the permanent staff 
have crystallized the investigations and studies pursued by 
the Workshop. In the course of this interchange, a number 
of questions and problems have arisen which merit further 
discussion. 

ANALYSIS OF RESIDUAL GASES 

To improve cryogenic pumping and space simulation, it is essential to know the 
nature and partial quantities of the vapors and gases evolved from the components of the 
vacuum chamber and the spacecraft. 

The mass spectrometer is an instrument that lends itself to analysis of the gases 
being evolved in a cryogenic pumping system at high vacuum. However, the gas analyzer 
capable of measuring minute quantities of gases or vapors over a range of molecules 
from 2 to 80 was constantly and continuously in a state of disrepair since the beginning 
of the Summer Workshop. Besides, the instrument is of limited molecular range, and 
some of the components given off by elastomers and plastics are of such high molecular 
weights that the instrument is not capable of detecting them. 

It is therefore desirable to consider other methods of analysis. One such method is 
gas chromotography, which requires larger quantities of gas samples than can be obtained 
in vacuum test-chambers at high vacuum; this, however, can be overcome by collecting 
samples on the fore-pressure side, where the gases evolved are compressed from 10 4 to 
10 6 times. It is essential to point out that samples of gases thus collected will also con- 
tain vapors of the diffusion oil, for which corrections could be determined. The applica- 
tion of this method to analysis of residual gases in high-vacuum chambers merits fur- 
ther investigation. 


CONTAMINATION OF SPACECRAFT SURFACE 

In addition to the common gases occurring in the atmosphere, which are given off in 
the course of evacuation by the different components of the chamber during a thermal- 
vacuum evaluation, gases and vapors are also being evolved which are characteristic of 
the components of the spacecraft and their histories. Steel, for instance, evolves-in 
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addition to the gases adsorbed from the atmosphere, such as water vapor, oxygen, nitro- 
gen and carbon dioxide-also carbon monoxide and hydrogen formed and absorbed in the 
course of its manufacture. The elastomers and plastics evolve a variety of different or- 
ganic substances related to the polymeric materials of which they are composed. Some 
of these gases, including hydrogen and carbon monoxide, may be irreversibly adsorbed 
on the metallized surface of the spacecraft and may change its properties. 

In this connection, it is of interest to point out that surface potentials are sensitive 
to surface composition, and lend themselves to measurement under conditions prevailing 
during simulation testing. Surface potentials with respect to an inert electrode can be 
measured by using the vibrating-electrode technique, or by means of a radioactive elec- 
trode that is an alpha-emitter. Both techniques are described in the literature. 

DETECTION OF LEAKS IN CRYOPANELS 

There have been instances of leakage in cryopanels only when liquid nitrogen is cir- 
culating in them. This was the case with the cryoliner in the cryosorption experiment 
described elsewhere in this report. 

One approach to detecting a low-temperature leak is to find a material that will dis- 
solve in, or mix with, liquid nitrogen and will produce an imprint on the surface of the 
cryoliner at the point of leakage. If iodine could be dissolved to any appreciable extent 
in liquid nitrogen, it would leave an imprint on a starch-painted cryoliner. 

There are perhaps to be found materials that fluoresce in ultraviolet. 

Another more promising approach is the use of radio tracers. A radioactive isotope 
of nitrogen, oxygen, or tritium would be injected into the circulating liquid nitrogen, and 
the area of leakage detected with a Geiger-Muller counter or scintillator. 


ORIENTATION OF SENSORS IN SPACE SIMULATORS 

In a gaseous system where the molecular mean-free-path is very much less than the 
dimensions of the vessel, intermolecular collisions are the major factors in the consider- 
ation of the physical properties of the gas. In such a gas, the vacuum sensor would not 
depend on its orientation, because the number of collisions per unit volume in a given 
position in the vessel would be the same in all orientations. 

However, in a region of pressure where the molecular mean-free-path is of the same 
magnitude as the dimensions of the vessel, or greater (the so-called region of free- 
molecular flow), the molecule-surface interaction is the governing factor. This condition 
prevails in space simulators. The kinetic energy of the molecules is characteristic of the 
temperature of the surface with which they interacted. Molecules originating from the 
surface of the spacecraft will generally be endowed with energies different than those 
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originating from the cryoliner, and will possess kinetic energies characteristic of the 
temperature of the spacecraft. For this reason, the orientation of the vacuum sensor is 
important . 

To derive a relationship between the pressures and temperatures of free-flowing 
molecules in these two orientations, assume equilibrium between the sensor and the 
spacecraft at which it is looking. The number of collisions per second with unit surface 
is given by 

1 

v = — nv, 

4 

where n is the number of molecules per unit volume and v is the average velocity. The 
mass striking unit surface per second is 

1 

mv - — pv 

4 

where m is the mass of a molecule and^> is the mass density. For the condition of equil- 
ibrium, the mass striking unit surface in either direction is the same, and therefore, 


Pl v l = P 2 V 2 

where the subscript 1 refers to the molecules originating from the sensor and the sub- 
script 2 refers to the molecules coming from the surface of the spacecraft. Substituting 
the kinetic expressions for_£ and v, we obtain 

Pi/^ p 2 /V^ 


where and are the temperatures of the sensor and of the surface of the spacecraft, 
respectively, and Pj and P 2 are the corresponding pressures. 
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CORONA AND DISCHARGE EFFECTS INSIDE SPACECRAFT 
AT LOW EXTERNAL PRESSURE 

PRELIMINARY STATEMENT 

(F. C. Yagerhofer, Space Power Technology Branch, Spacecraft Technology Division, GSFC) 

The problem in corona will be discussed in three parts: The first will be a review 
of the background of the corona problem; the second will be a summary of present space- 
craft corona problems; and the third, methods and techniques for the detection and sup- 
pression of corona. 

First, let us glance into the background of the corona phenomenon. Corona can be 
defined as that ionization phenomenon which is produced when the ionization potential of 
a gas or air is exceeded. Corona forms small quantities of ozone, which attacks insu- 
lation materials, and oxides of nitrogen which corrode metals. This problem was first 
encountered years ago in the electrical distribution industry. Even nowadays, the engi- 
neers who design overhead-distribution lines will allow a value of about 1000 watts of 
power loss per mile for a fair-weather day. If the weather gets foul, in rainy or humid 
weather, for example, this power loss can readily go up by a full order of magnitude. 

Another problem encountered years ago was in the aircraft electrical industry. 
During the war, it was found that dynamotors used aboard the aircraft to change 28 volts 
to over 400 volts as necessary to power the navigational and surveillance equipment on 
the airplane had a tendency to generate corona. This shortened the life of the insula- 
tion, and at the same time generated a high-frequency component which was reflected 
into the electronics aboard the aircraft. 

Although a lot of vile things have been said about the corona, it can nevertheless be 
employed to a useful end. At the Naval Research Laboratory, we tried to use the corona- 
starting voltage as a method for studying the quality of insulation materials. We wanted 
to establish a nondestructive test to depict adequately the true condition of the insulation 
in, let's say, an aircraft generator. We found in our tests that the breakdown potential 
was approximately twice the corona- starting voltage (CSV). We determined that this 
corona-starting voltage test would be a valid nondestructive test for periodically testing 
the quality of insulation, if we could find some method or means to precisely detect the 
existence of corona. A number of major problem areas remain in the study of corona; 
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in fact, as far as can be determined, there really is no universally accepted definition for 
corona even today. 

Item 2, a summary of our present spacecraft problems: A number of experimenters 
here at Goddard have encountered corona. Leo Davis, for example, in his rocket probes, 
encountered corona in a 2000-volt power supply. He figured, though, since this occurred 
at altitudes corresponding to 60,000 kilometers, the corona did not exist for a sufficient 
length of time to create a carbonized path, or anything of that sort. So, since it didn't 
occur at sea level, nor at the extremely high altitudes, the corona was tolerated. But in 
his equipment used aboard spacecraft, he encapsulated these critical components in a 
hard epoxy in order to minimize or eliminate this corona. 

Another incident occurred during the thermal vacuum test in the UK-1 satellite, in a 
power supply that supplied 40,000 volts for the X-ray ionization chambers. This power 
supply used a plug-in connector to attach the external wiring into the power supply, and 
arc -over occurred right at the connector. Now, somehow or other, this gaseous dis- 
charge got into some of the transistors located in close proximity and damaged them. A 
puncture occurred from the emitter to the base in one of the transistors. The fix in this 
instance consisted in removing the entire connector and installing a solid conductor with 
a heavy insulation, plus an additional insulating sleeving, which eliminated the problem 
from this particular experiment. 

Mr. Flanick encountered corona in his high-voltage power supplies for the Geiger- 
Muller as well as for the photomultiplier tubes. He employed the standard technique of 
encapsulating the critical components in a hard epoxy; in some instances, he used a rub- 
berized silastic compound to eliminate corona. 

In our own case we designed a 6,000-volt power supply to be used for the Redhead 
gauge in the S-6 atmospheric structure satellite. During vacuum tests, corona appeared 
in the transformer core. We increased the thicknesses of the Teflon insulation between 
the layers of this saturable core and potted it with a foster-cast polyester, which elimi- 
nated the corona. 

This brings us now to our third item, the methods and techniques for corona suppres- 
sion. There are three basic methods: one, by pressurization or evacuation. Pressuri- 
zation, as you probably know, is employed for the wave guides used on high-powered 
radar equipment. A second way is to increase the radius of curvature of the current- 
carrying conductor. The voltage gradient at any point on the surface of a conductor 
varies approximately inversely with the radius of curvature. Thus, to eliminate dis- 
charge, increase the curvature of all points at high voltage. This system, by the way, is 
used in our high-voltage power-distribution lines. You find in a 500-kv line a 2-inch OD 
conductor; this is really a thin-walled tubing, with about a quarter -inch or an eighth of 
an inch wall. Now, this large diameter isn't necessarily for carrying heavy currents - 
it's to eliminate, or at least minimize, the corona power loss. 
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The third method which I mentioned earlier, is to encapsulate or to immerse the 
critical components in a high dielectric material such as an epoxy, rubberized silastic, 
or oil. Any of these methods can be used, but there must be no contaminants on the sur- 
face of the conductor, and no voids in which ionization can form between the wire and 
insulation. 

I have briefly covered a few of the problems that we have had in the past. A number 
of problems remain to be solved, particularly an accurate means of detecting corona, as 
well as improved dielectric-strength materials capable of withstanding satellite 
environment. 
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ON THE BREAKDOWN VOLTAGES OF SOME 
ELECTRONEGATIVE GASES AT LOW PRESSURES 


by 

Stefan Schreier 

ABSTRACT 

The breakdown voltages of five Freon gases and SF g were 
determined at low pressures between parallel plates, with a view 
to determining their insulating properties under these conditions, 
and especially to find their minimum breakdown voltages. It was 
found that the advantages of these gases over air is much less at 
low pressures than at atmospheric pressure or higher, and that 
the use of these gases as insulators at lower pressures is there- 
fore limited. Preliminary investigations reveal, however, that 
the vapors of compounds which are liquid at normal temperature 
and pressure, such as FC-75, maintain their dielectric strength 
much better than do electronegative gases as the pressure is 
reduced, and thus seem to hold promise as insulators at lower 
pressure. 

INTRODUCTION 

Certain artificial earth satellites, having direct-current power supplies in the range 
of 1,600 to 2,000 volts, were found to have difficulties caused by corona-type discharges 
when passing through altitudes in the vicinity of 100,000 feet. These difficulties occurred 
because the equivalent pd (p equals pressure and d equals distance between electrodes) 
in the satellites passes through the range of minimum breakdown voltage, which in air is 
about 350 volts. No difficulties occurred at either higher or lower altitudes, as might be 
expected from Paschen's Law, which states that the breakdown voltage between two elec- 
trodes is a function of the pressure times the distance between the electroces and de- 
creases with this product down to a certain minimum, after which it rises again as pd 
decreases further. 

The exact shape and location of the curve of pd versus breakdown voltage will depend 
on the shape of the electrodes, the material of which they are made, and the nature of the 
dielectric between them. 

Various suggestions have been made for preventing corona-type discharges in satel- 
lites. The solution presently in use consists of imbedding all electrodes (charged con- 
ductors) in a solid dielectric, thus making gaseous discharges impossible. This has the 
disadvantage that, once conductors are so imbedded, they are difficult to remove. Other 
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suggestions include the use of electronegative gases as insulators (the subject of the 
present investigation), pressurizing high-voltage components, immersing them in oil, 
evacuation of the entire system, keeping the power off until the critical values of pd have 
been passed (delayed switching), and artificial suppression of free electrons in the gas. 
Some of these suggestions may well bear further investigation. Other possible avenues 
of investigation are the seeking of improved design criteria to delay corona, and the 
further study of the effects of electrode materials and configurations with a view not only 
to prevent, if possible, the onset of corona, but also to provide the designer with im- 
proved methods of predicting when corona may occur. 

In order to solve the satellite problem, it was proposed to enclose the high-voltage 
components in an electronegative gas which would be allowed to leak out as the satellite 
rose, maintaining the pressure inside the satellite about the same as the pressure out- 
side. This would eliminate the need for pressurization and the danger of unintentional 
leaks, and would prevent corona-type discharges as long as the minimum breakdown 
voltage of the gas was higher than the maximum voltage expected in the satellite. For 
this reason, a study was undertaken to find the minimum breakdown voltages of certain 
electronegative gases. 


PREVIOUS WORK ON ELECTRONEGATIVE GASES 

As is well-known, electronegative gases obtain their dielectric strength from the 
fact that their molecules have a tendency to absorb loose electrons, thus inhibiting the 
electron avalanches necessary for breakdown. This is already discussed briefly in Meek 
and Craggs. 1 A somewhat longer summary has been given by Devins and Sharbaugh. 2 
The effect of pressure on the positive-point-to-plane discharge for some electronegative 
gases was investigated in 1939 by Pollack and Cooper. 3 After the war, further investi- 
gations were undertaken by Camilli and Chapman 4 and Wilson Simons and Brice. 5 Fur- 
ther work was done on fluorine-containing gases by Camilli and Plump in 1952. 6 Addi- 
tional theoretical considerations were given by Geballe and Reeves in 19 53. 7 The electric 
breakdown of perfluorocarbon vapors and their mixtures with nitrogen were discussed by 
Berberich, Works, and Lindsay in 1955; 8 the dielectric behavior of some fluorogases 
and their mixtures was also discussed in the same year by Camilli, Liao, and Plump. 9 
Additional discussion of factors controlling electric strength of gaseous insulation was 
given by Narbut, Berg, Works and Dakin in 1959. 10 The advantages of gas insulated 
power transformers were discussed by Camilli in 1959. 11 In the same year, Blodgett 
evaluated some of the dielectric properties of octafluorocyclobutane. 12 A further discus- 
sion of fluorocarbons as electrical insulators was given by Reuther in 1961.' 3 

The above is by no means a complete listing, but is merely given as an illustration 
of some of the things which had been done in this field. Further references may be found 
in AIEE Special Publication S-97. 13 The difficulty with most of the previous work is that 
it is concerned chiefly with the insulating properties of gases at atmospheric or higher 
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pressures; very little had been said about breakdown at lower pressures, and especially 
about minimum breakdown voltages. 


GASES CONSIDERED IN THE PRESENT STUDY 

The gases considered in the present study were Freons 14, 114, 115, 116, and C318. 
SF 6 was chosen because it is perhaps the best known and most used of the electronega- 
tive gases. The newer Freons were evaluated because little is as yet known about their 
properties, and it was hoped that they would prove superior to their older relatives. 

FREON 14. Freon 14 is perfluoromethane, CF 4 . It has a molecular weight of 88.01, 
boils at -128°C, and freezes at -168°C at one atmosphere pressure. Its dielectric strength 
at one atmosphere and 23°C relative to nitrogen is 1, and its dielectric constant is 1.0006 
at 24. 5° C and one atmosphere. 

FREON 114. Freon 114 is CC1F 2 -CC1F 2 . It has a molecular weight of 170.93, 
boils at 3.55° C and freezes at -94°C at one atmosphere pressure. Its dielectric strength 
relative to nitrogen is 2.8 at one atmosphere and 23°C, and its dielectric constant i^ 

1.0021 at 26.8°C and one atmosphere. 

FREON 115. Freon 115 is chloropentafluoroethane, CC1F 2 -CF 3 . It has a molecular 
weight of 154.48, boils at -38°C, and freezes at -106°C at one atmosphere pressure. Its 
dielectric strength relative to nitrogen is 2.8 at one atmosphere and 23°C, and its dielec- 
tric constant is 1.0018 at 27.4°C and one atmosphere. 

FREON 116. Freon 116 is hexafluoroethane, C 2 F 6 . It has a molecular weight of 
138.0, boils at -78.2°C and freezes at -100.6°C at one atmosphere pressure. Its dielec- 
tric strength relative to nitrogen is 1.96 at one atmosphere, and its dielectric constant 
is 1.00197 at 23°C and 711 mm/Hg. 

FREON C318. Freon C318 is octafluorocyclobutane, C 4 F g (cyclic). It has a molecu- 
lar weight of 200, boils at -6.0°C and freezes at -41.4°C at one atmosphere pressure. Its 
dielectric strength relative to nitrogen is 2.63 at one atmosphere, and its dielectric con- 
stant is 1.0034 at 10°C and 760 mm/Hg. 

The properties of sulfurhexafluoride, SF fi , are well-known. Its molecular weight is 
146, it boils at -63.8 C at one atmosphere, and its dielectric strength relative to nitrogen 
is 2.5 at one atmosphere. 


THE EOUIPMENT USED 

The equipment used was a variant of that used by Dakin and Lim 15 of Westinghouse 
and that of Brancato* of the Naval Research Laboratory (Figures 4-1 and 4-2). Two 
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(a) Total setup, showing scope, variacs, vacuum 
gauge, voltmeter, bell jar, and pressure gauge 


ectrodes 


Figure 4-1 - experimental equipme 
determining breakdown voltage of e 
negative gaser 
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Bell Jar 



Figure 4-2 - Experimental equipment, block diagram 


parallel circular brass plates, 4 inches in diameter and 1/4-inch thick with rounded 
edges of 1/8-inch radius, mounted on a teflon base, were used as electrodes. These were 
connected to a 3,000-volt transformer connected to an ordinary 60-dycle 110-volt power 
supply. 60-cycle alternating current has the advantage over direct current that, while it 
yields the same type of corona as direct current, because of its relatively slow oscilla- 
tion, it triggers corona at each peak of the cycle on a regular basis at the proper voltage 
and thus makes it much easier to study the phenomenon (Figure 4-3). Filters were used 
both to suppress high-frequency noise from the power supply before it reached the elec- 
trodes, and to filter out the 60-cycle sine wave before feeding the signal to the oscillo- 
scope. The presence of corona was determined by observing the signal from the elec- 
trodes on a Tektronix type 545A oscilloscope, using a type L plug in unit preamplifier 
with a fast rise time. The electrodes themselves were placed in a bell jar which was 
evacuated by means of a mechanical pump capable of reducing the pressure in the bell 
jar to 45 microns of mercury. Voltages across the electrodes were read from a care- 
fully calibrated Simpson vacuum-tube voltmeter. The pressure inside the bell jar was 
read from a Wallace and Tierney absolute-pressure gauge calibrated to read NACA 1956 
standard atmosphere to 200,000 feet; from 1 mm/Hg on down, the pressure was also read 
from a Consolidated Vacuum Corporation thermocouple vacuum gauge, giving a double 
check in the range of equivalent altitude from about 150,000 to 200,000 feet. The filters 
were immersed in transformer oil to prevent the possibility of local discharges occuring 
in the circuit. 


THE EXPERIMENTAL PROCEDURE 

The equipment was first tested in air (Figure 4-4) to make sure it was operating 
properly and to provide a basis for comparison with subsequent readings, as it is well- 
known that the exact shape of the Paschen curve tends to depend on the particular con- 
figuration and equipment being used. To test other gases, the bell jar was first evacuated 
to a pressure of 60 microns of mercury; enough gas was then introduced to raise the 
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(a) Corona discharges being triggered 
by 60 -cycle sine wave 



(b) Single corona pulse, showing initial spike 
(double exposure) followed by reverberation 



(c) Detail of reverberation following corona 
pulse (about 1 millisecond) 

Figure 4-3 - Phenomenon of corona discharge 
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pressure to 60 mm/Hg or more, thus guaranteeing a purity of 99.9% in the gas being 
tested. After the pressure was recorded, the voltage was raised across the electrodes 
until breakdown was observed on the oscilloscope. The voltage was then backed off until 
the discharge ceased and then raised again in small increments until corona was again 
observed. This procedure was followed to prevent errors due to time-lag effects. For 
each reading, the entire procedure was followed twice to make sure the reading was 
correct. Once the reading was obtained for a given pressure, enough gas was pumped out 
to lower the pressure by about 10,000 feet equivalent altitude (around 5,000 feet near the 
minimum voltage), and the next breakdown voltage was obtained. Three runs were made 
for each gas, using a gapwidth between the electrodes of 1 mm, 2/3 mm, and 1/3 mm. 


CONCLUSION 

The results show that the electronegative gases tested lose most of their advantage 
over air as insulators at lower pressures. Freon 114, which has one of the highest diel- 
ectric strengths of any of the Freons at atmospheric pressure, was found to have the 
lowest minimum breakdown voltage of any of the gases tested (about 435 volts). SF 6 was 
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found to have the highest minimum breakdown voltage (520 volts) as compared to a mini- 
mum for air of 353 volts, in the present configuration. Note that the relation between 
breakdown voltage and molecular weight mentioned by Camilli 11 ceases to hold for the 
minimum breakdown voltage. 

The results of these tests, shown in Figures 4-5 through 4-10, indicate that gases of 
the type tested are not too promising as insulators at lower pressures in satellites. On 
the other hand, preliminary tests indicate that vapors of such liquids as FC-75 seem to 
keep their dielectric properties much better than gases at lower pressures. Two new 
Freon compounds currently under development, which are liquid at normal temperature 
and pressure, also hold promise; they are 1,1,1,3-tetrachlorofluoropropane, and 1,1,1- 
trichloropentafluoropropane. The former especially is claimed by the manufacturer to 
have a dielectric strength in vapor form 4.7 times that of FC-75. These vapors would 
seem to merit further investigation. 



4-12 




CORONA AND DISCHARGE EFFECTS 



Figure 4-6 - Paschen curve. Freon 114 
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Figure 4-10 - Paschen curve, SF 6 
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MAGNETIC FIELD SIMULATION 

PRELIMINARY STATEMENT 

(By Dr. J. Heppner, Fields and Particles Branch, Space Sciences Division, GSFC) 

As is well-known, the earth’s field is roughly a dipole, and some of the principal 
areas of investigation in space science are the magnetic field, the radiation belts, and 
so forth, so that quite a bit of emphasis goes into measuring the magnetic field. If we 
consider the requirements, say, of a satellite instrumented with magnetometer, we’ll hit 
on problems encountered in almost all other satellites. 

Under the category of magnetometer test, it is desired to map the magnetic field at 
the surface of the earth, where it ranges from around 31,000 gammas (1 gamma = 10 5 
gauss) at the equator to about 65,000 gammas at the poles. On the other hand, at great 
distances from the earth, the fields may interfere with and cancel each other. Hence, in 
order to be sure that magnetometers will operate under such conditions, magnetic -field 
simulation must be capable of ranging anywhere from one gamma up to, say, 65,000 
gammas. 

To calibrate these various types of magnetometers and to determine their frequency 
response, they must have a very homogeneous field in which to operate, certainly much 
better than the field existing in the air around us. 

Since a magnetometer generally must be located some distance away from the space- 
craft, as on a boom, so that the fields may not be severely distorted by the spacecraft 
itself, additional problems arise. At this distance, it will have thermal problems all its 
own; it is not part of the overall heat balance of the system, and its thermal condition 
fluctuates wildly due to its very small heat capacity. This would require that it be tested 
in a nonmagnetic thermal -vacuum chamber. 

After a magnetometer has been placed in a certain position, it must stay there. 
Components of the rest of the spacecraft (and even the whole package) must be tested to 
make certain they do not interfere with the operation of the magnetometer. If interfer- 
ence does occur, the magnetic materials and current loops within the spacecraft must be 
eliminated. 
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As the spacecrafts grow larger, a very large volume must be sought over which 
there is a controlled homogeneous magnetic field; this leads to rather large coil sys- 
tems, with very precise controls and very fine tolerances in the coils. 

These same considerations apply to a lesser degree to practically all payloads, one 
reason being that many satellites -sounding rockets especially-use small magnetometers 
as an inexpensive method of determining orientation up to a few degrees of accuracy. 

This isn't too bad, as long as the rocket is going up only a few thousand kilometers or 
less. Nevertheless, magnetometers still must be located properly for proper operation. 

There is another category, namely, satellite orientation control, that is, using the 
magnetic field to control the orientation of the vehicle. When you put in a large magnet, 
which then locks into the magnetic field, such as in TRANSIT, one will want to know how 
well it's going to lock in. One must also know beforehand the moment of the payload, 
etc. This is control in the passive sense. 

Semiactive systems, like TIROS, are designed to change the orientation of the 
spacecraft. A large current is put through a coil and the torque between that field and 
the earth’s field produces a transient change in the position of the satellite. In this way 
it can be reoriented, and pictures can be taken over more of the orbit. In these last two 
cases one must be able to take measurements on the variability of the magnetic torque. 

Then there are systems, such as the OAO, where the magnetic field is used in an 
active sense; measurements of the magnetic field are made, and that information is used 
to operate gas jets, etc., to orient the vehicle. 

One of the important things is to know the variability of the field and how this affects 
the payload. The magnetic properties of the materials inside haven't changed, but their 
position on a B-H curve has changed, and measurements must be made to find the whole 
range of variation to determine its limits. 

The spinning of a satellite, especially if it is small and metallic, introduces new 
problems. When a satellite spins in a magnetic field, the spin damps out rather fast if 
there is much of a surface. In fact, fields of 5,000 gammas or so can be generated in- 
side a metal ball like a typical Vanguard satellite, just due to the fact that there is a 
conducting surface spinning in a magnetic field. So spin damping becomes another of 
these problems. 
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MAGNETIC FIELD SIMULATION 

by 

D. L. Waidelich and M. E. Pittman 


ABSTRACT 

Several new and, in some respects, better solutions have been 
obtained relating to the production of homogeneous magnetic fields 
by the use of circular coil systems. 

Previous work exists in many forms, mainly giving particular 
solutions in which the author is interested, such as setting all cur- 
rents equal or setting the radii equal. 

The work deals with the calculation leading to two particular 
types of optimal systems: the optimum for the largest coil radius 
and the optimum for the mean coil radius. 

One calculation is presented on the tolerances necessary 
when one coil is moved from its optimum position. 


INTRODUCTION 

In testing the magnetic properties and instruments of spacecraft, there is a need to 
cancel the earth's magnetic field and then to produce a very homogeneous controlled 
magnetic field. It is desired to produce a magnetic field homogeneous to 1 part in 10 5 , 
such that there is ample access to the working volume. The electrical design and con- 
struction should be simple, and the largest possible mechanical tolerances should exist. 
It is the purpose of this report to review briefly some of the previous work on circular 
coil systems, and to present the results of calculations made on these systems during 
the Summer Workshop. In addition, some work on tolerances will be given. 

THEORY OF IDEAL COIL SYSTEMS 


In his paper, McKeehan 1 reviews many of the more important particular solutions 
for producing a homogeneous magnetic field using two, three, four, six and eight loops, 
and also gives several new solutions. Our equations and terminology are the same as 
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his. In some cases the distances to the periphery of the coils are kept constant. This 
was done by Maxwell 2 for 3 coils, by McKeehan 1 and Garrett 3 for 4 coils and by 
Williams and Cain 4 for six and eight coils. Fanselau 5 ’ 6 and Braunbek 7 held the currents 
equal. McKeehan 1 used equal currents in a six-coil solution. Barker 8 ’ 9 solved for 3 and 
4 identical coils. Additional work on these systems is given in Sauter. 10 

Two methods of attack were considered; one 11 ’ 12 used the idea of making the field 
the same at a number of points along the axis, but did not control the magnetic gradients 
that exist. The second 1 used the idea of making as many terms zero as possible in the 
series for the magnetic field on the axis of the system. This method, which afforded 
some control of the magnetic gradients, was the one used here. Both circular 1 and rec- 
tangular 11 coil systems were considered, but it was decided to work on the circular coils 
because more work had been done on these and there is more symmetry present; in addi- 
tion, rectangular systems are more difficult to calculate. The problem of calculating the 
circular coil system may be worked by the use of elliptic functions, or by an infinite 
series of Legendre functions; the use of the elliptic functions proved to be quite difficult, 
and the method using Legendre polynomial expansions was adopted. 

One other method of obtaining a uniform field was considered: this was the use of a 
spherical current sheet 13 ’ 14 in which the current density varies as the sine of the angle 
that the radius vector makes with the axis. This method, while giving a theoretically 
perfectly uniform field, would allow an extremely limited access or would require the 
use of many connections to permit access. It was felt that the disadvantage of poor ac- 
cess ruled out further consideration of this method. 

The scalar potential due to a single coil, whose periphery is described in polar 
coordinates as r p and at a point r and o is given by: 


V - 

p C 


r 00 

'l+u -(l-u 2 ) / -A. (r / r ) n p' (u ) P n 

| P v P 7 Z— i n K P n v p 7 n 


(u) 


n - 1 


FOR r < r p . 


(i) > 


Here V p is expressed in Gaussian units, i p is the current in the pth coil, c is the 
velocity of light, and u and Up are used in place of cos <•' and cos P , respectively. P n (x) 
and P' n (x) are the nth Legendre polynomial and its first derivative, respectively. This 
is illustrated in Figure 5-1. The total potential v equals I. v p . 

Since |p n (u)| <1 and equals one if and only if u = 1 , it is felt that, if the main in- 
homogeneities on the axis are eliminated, then the rest of the field will be at least as 
well off as points on the axis. It now becomes easier to work with the field itself on the 
axis. If x is the distance (on axis) from the origin, then 
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H(X) = ^Z/ n {2]p (1 " U p 2)P n + l (U P ) / r p n + 1 } (2) 

n “ 0 L P J 

where H is the component of the field parallel to the axis. 

As long as the coils occur in pairs, symmetric about the origin, then the parity of 
the Legendre polynomials will eliminate the inner brackets for n odd, and the sum need 
be taken only over even terms. The field will be of the form: 


00 

H(x) = ^ a 2n x2n = a o + a 2 x 2 + a 4 x4 + •••• 

n - 0 


( 3 ) 


Perfect homogeneity is obtained if a 2 , a v etc., are all identically zero, as in the case of 
the so-called sine-windings on a sphere. 

Approximations can be made here so that a finite number of terms, a 2 , a 4 , . . . a n 2 
are zero. To find the volume of homogeneity, rewrite the equation as: 


H(x) = a Q 


HSer*-)- 


where r = mean coil radius = « r x • r 2 


( 4 ) 
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Neglecting higher order terms, we see that (for .001%) 




S 10-5, or*< 


r* n — • 10' 5 
a 


1/n 


( 5 ) 


Such an expression as x/r must be used rather than x/r 1? etc., to prevent the ratio 
from going from infinity to zero, as r x goes from zero to infinity leaving the other coils 
constant. 

The specification for which x/r is a maximum for a particular number of coils, is 
called the optimum for the mean-coil radius. 

Another useful quantity to consider is x/ rmax where r max is the largest r p of the sys- 
tem. The point at which this quantity is a maximum is called the optimum for the 
largest-coil radius. If one is limited to a certain largest volume in which to operate, 
then obviously this would be the better of the two to use. Unfortunately, this involves 
changing functions as the ratio of radii crosses unity. It is possible to define other op- 
tima as well, but it is believed that the two mentioned above are the most useful and 
practical. 

PARTICULAR CASES 

For specifications of all coils mentioned, consult Table 5-1. 


EXPLANATION OF TABLE 5-1 

The name ’’efficiency optimum" represents the optimum for the mean-coil radius. 

The ratios of the currents and radii are given in the table, whereas the magnitudes are 
determined by the amount of field desired and the volume of workroom desired respec- 
tively. The field at the center is given by: 

«<.>■?£*. <‘-.Vv 

P 

As mentioned previously, 7 = ( ri . r 2 -..r n ) 1/n , noting especially that, for three coils, 
r = (j-j • r 2 2 ) 17 5 Furthermore, i p is actually equal to the number of turns in each wind- 
ing time the current in each turn, so that if the same current flows in all wires, i 2 ' i j = 
n 2 /Nj, etc. 

To determine H(x) at any point on the axis, the table is sufficient. As an example, 
consider the Helmholtz coil, with a 2 equal to zero. H(x) is then proportional to 
(1 - 1.8 (x/r) 4 ). 
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Table 5-1 

Specifications of Various Coil Systems 


Name and Number of Coils 

Ampere (1) 

Helmholtz (2) 

Barker (3)* 

Maxwel 1 (3)* 

Ef f ic iency Opt. (3)* 

Assumpt ions 

none 

none 

Z 1 = z 2 

r l = r 2 

none 

Which a’s are zero 

none 

a 2 

a 2 /a 4 

a 2 ' a 4 

a 2 ' a 4 

Coefficient of next derivative 

-1.5 

- 1.8 

-5.0695298 

-2.0428571 

-2.036426 

U 1 

0 

±.447214 

0 

0 

0 

u 2 

- 

- 

±.760051 

±.6546537 

±.61627 

r 2 /r J 

- 

- 

.796143 

1 . 

1.19697 

1 2 /l 1 

- 

- 

2. 2604 

1.53125 

3.075846 

x/7 

.26% t 

4.86% 

11.30% 

13.031% 

13.037% 

x/r max 

.26% 

4.86% 

9.71% 

13.031%* 

11. 564% 

Name and Number of Coils 

Fanselau (4) 

Braunbek (4) 

Barker (4) 

McKeehan (4) 

Ef f ic iency Opt. (4) 

Assumpt ions 

ij = i 2 

1 1 = x 2 

Z 1 =z 2 

r l = r 2 

none 

Which a’s are zero 

a 2' a 4 

a 2 ' a 4 * a 6 

a 2 ,a 4 ' a 6 

3 2 ' a 4 ' a 6 

a 2' a 4 ,a 6 

Coefficient of next derivative 

-1.289 

-2.23448 

-2.4725 

-2.255102 

-2.2219679 

U 1 

±.285232 

±.26789 

±.243186 

±.2852315 

±.2723547 

u 2 

+ .765055 

±.74213 

±.940731 

±.7650553 

±.7484183 

r 2 /r l 

1.136009 

1.098 

1.336 

1 . 

1.0712777 

1 2 Z 1 1 

1 . 

1 . 

2.2604 

.6821109 

.9040608 

x/7 

14.07% 

21.682% 

21.480% 

21.607% 

21.677% 

x/r „ 

max 

13.20% 

20.682% 

18.583% 

21.607* 

20.944% 

Name and Number of Coils 

Neumann ( 6 ) 

Braunbek- 
McKeehan ( 6 ) 

Cain- 

Wi 1 1 iams ( 6 ) 

Neumann 
McKeehan ( 8 ) 

Cain- 

Wi 1 1 iams ( 8 ) 

Which a’s are zero 

a 2 ,a 4' a 6 

a 2 to a io 

a 2 to a io 

a 2 to a g 

3 2 ^ ® a 14 

U 1 

±.20929922* 

±.190655 

±.2 0 9 2 9 922 

±.1652754* 

±.1652754 

u 2 

+.59170018 

±.550274 

±.59170018 

±.4779250 

±.4779250 

u 3 

±.87174003 

±.843307 

±.87174003 

±.7387739 

±.7387739 

u 4 

- 

- 

- 

±.9195342 

±.9195342 

r 2 /r l 

1.071723 

1.046147 

1 . 

1.0222398 

1 . 

r 3 /r 1 

1.242359 

1.157907 

1 . 

1.1827288 

1 . 

r 4 /r 1 

- 

- 

- 

1.2382935 

1 . 

l 2 /l 1 

1 . 

1 . 

.8270469 

1 . 

.891626 

i 3 /i l 

1 . 

1 . 

.5108492 

1 . 

.686604 

Vm 

- 

- 

- 

1 . 

.406992 

x/r 

max 

31% 

34% 

36% 

16% 

47% 


♦Current in middle coil = 2ij. 
fWhere inhomogeneity = 10 - 5 . 

❖Volume optimum, the optimum for the largest coil radius. 
❖Zeroes of P 7 ' (u) . 

♦ Zeroes of P 9 ' (u) . 
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A convenient way to find x/F (or x/r max ) for inhomogeneities other than .001% is to 
multiply the figures given by io (5 + a)/b where a^is the logarithm of the desired inohom- 
geneity, and b is the subscript of the first nonvanishing a. 

To find where x/r crosses .1% for the Helmholtz coil, a = log 10" 3 = -3, and b = 4 . 
Therefore, x/F = 4.86% • 10 1/2 = 15.37%. 

For convenience: 

10 12 = 3.1622777 
10 1/4 = 1.7782794 
10 1/6 = 1.4677990 
10 1/8 = 1.3335214 
10 1712 = 1.2115275 

TWO COILS. Since only one variable is free, only one term can be made to vanish, 
and that is a 2 . The criterion, due to Helmholtz, is that u p 2 = 0.2. For such a system, 
x/F = 4.9%. 

THREE COILS. The third coil is placed at u = 0 to maintain symmetry. a 2 and a 4 
can be set equal to zero, as does Maxwell 2 and Barker. 5 If, however, the term |a 6 |F 6 /a o 
is determined in terms of one of the variables and minimized analytically, a slightly 
better system results. 


f(y) = I a 6 | r 6 /a Q oc (15y 3 - 9y 2 + 5y + 5)/(y + 1). 


( 6 ) 


where y = u 2 2 . 

Then f' (y) = 0, yields u 2 2 = .4l£ - .6, 
or u 2 ~ .61627 

x/F then equals 13.037% as compared with 13.031% of Maxwell's coil and 11.30% as cal- 
culated from Barker's. 

In studying the optimum for the largest coil radius, x/r M , it was found to coincide 
with Maxwell's solution to within six decimal places. It should be mentioned that the 
knowledge gained in solving the three-coil problem was very useful in planning the four- 
coil solution. 

FOUR COILS. With the aid of an IBM 7090 computer, again slightly better results 
were obtained than previously found. a 2 , a 4 , and a 6 all vanish. After an unsuccessful 
attempt to make a 8 vanish, the following scheme was undertaken: 
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a 2 = 0 = P' 3 (u 1 )-kP' 3 (u 2 ) 
a 4 = 0 = P' 5 (u 1 )-k/r 2 P' 5 (u 2 ) 
a 6 = 0=P' 7 (u 1 )-k/r 4 P' 7 (u 3 ) 

where r = r 2 /r j and k = i 2 (1 - u 2 2 )/i j r 3 (1 - u j 2 ). (7) 

Eliminating k and r yields: 

P' 3 (Ui)P' 7 (Ui) P’ 3 (u j) P' 7 (u 2 ) 

_ k (u 2 ) 

[P'gCuj )] 2 [P' 5 (u 2 )] 2 


or 


p 3 (u i) P' 7 ( Ul ) - F(u 2 ) [P'sCuj)] 2 = o. (8) 

Over the entire range of u 2 , this equation was solved for u l , and then k and r were 
calculated. Literally hundreds of different solutions were obtained. These solutions in- 
cluded the old ones of McKeehan, Braungek, and Barker, as well as many new ones. It 
was then decided to minimize the next two terms in the series for added accuracy. 


a o C 1 + r ' 8a 


8 /a o (x/r) : 


8 + F 10 a 


10 


/a„ (x/F) 10 ) 


and 

(x/F)’ 8 = F 8 /a 0 |a g +F 2 a 10 (x/F) 2 |. 


( 9 ) 


Iteration procedures yield that 


x/F = 21.677%. 

This compares to 21.607% as calculated from McKeehan's 1 data; 21.671% of Braunbek ; 7 
21.480% of Barker ; 8 and 14.07% of Fanselau . 5 Fanselau's solution makes a 2 and a 4 zero, 
a 6 is not zero. The other solutions make all three a’s zero. 

Again the optimum for the largest-coil radius coincides with the solution for equal 
radii, that of McKeehan. 
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Another interesting result is that of Fanselau 6 in which the two coils to the one side 
of the center of the system are placed in one plane, while the other two coils are placed 
in another parallel plane. This arrangement may have considerable advantages in struc- 
tural design and accessibility to the object being tested. The disadvantage would be a 
smaller x/7 or x/r max than that of the other solution of Fanselau 5 given in Table 5-1. The 
results for this solution are u 1 = ±0.285232, u 2 = ±0.765045, r 2 / ri =±0.372820, i 2 /i 1 = 
28.2897. Scott 17 may have a similar solution. The solutions of Fanselau 6 and Scott 17 
appear to make a 2 and a 4 zero. It is believed that, among the many four-coil solutions 
obtained by use of the computer, there is one that will place two coils in the same plane 
and the other two coils in another parallel plane. This solution would have the advantage 
of making a 2 , a 4 , and a 6 zero, and thus should have a larger volume of uniformity than 
the solutions of Fanselau 6 or Scott 17 , although smaller volume than the optimum four- 
coil solution given in Table 5-1. It is hoped to pick out this solution and add it to Table 
5-1 later. 

SIX COILS. A computer program was also set up to solve the six-coil problem. 

Since seven variables are available (three angles, two current ratios, and two radii 
ratios), an attempt was made to set six of the derivatives equal to zero. Previous ef- 
forts (l,4,16) have made a maximum of five vanish. 

Three of the variables can be easily eliminated from the first three equations, and 
substituted into the remaining derivatives. This leaves the problem as three equations 
in three unknowns, plus one more which could be set at an arbitrary value. If u x 2 is less 
than u 2 2 and u 3 2 , then its range is limited between 0 and .2. u x 2 was to be taken over 
this range in steps of .001. 

Upon meeting an equation such as: f( x , y, z) = 0, and guessing approximate solutions: 
x o , y Q > z o ; a correction to the variable x could be obtained by letting f(x o + Ax, y D , z o )=0 

0 = f (x + Ax, y , z ) = f(x ,y , z ) + f'(x ,y ♦ z ) A x + • • • » 
x o oo ooo ooo 


where the prime denotes derivatives with respect to x. 
This gives us the first order correction to x o : 


Ax = -f(x o , y o , Z 0 )/f'(x 0 - y o - z o> 

The next approximation to be used is x x = x o + Ax. Various loops were designed 
into the program, for all three variables to zero in on the correct values. When the 
criterion was met that the corrections to the three variables were all less than 10' 7 , the 
computer would calculate the seventh derivative with the eventual hope of minimizing it. 
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RESULTS 

In this particular program, numbers frequently arose which were too large or too 
small for the computer to handle. As a result, traps were built into the program to elim- 
inate these particular values of u t 2 . This happened in 44 of the 200 values. In another 
21 cases, the ratios of the radii went imaginary, so these values were eliminated. 132 
values did not register, because they were apparently diverging. This was previously 
taken care of by setting a maximum to the number of times each variable could be incre- 
mented. In the remaining three cases, "solutions" were obtained. Here current ratios 
were on the order of 10 2S , and radii ratios were about 10 3 . 

An example is the following set of data: 


”1 2 = 

.02900 

u 2 2 = 

.67269 

”3 2 = 

.71167 

r l /r 2= 

.0003897 

r l /r 3= 

.0004642 

i 2 /i l = 

1.05-10 24 

i 3 /i l = 

-1.64-10 24 


As a result, these also should be thrown out, due purely to the ridiculous orders of 
magnitude. 

Despite the very pessimistic results obtained, it is still felt that a much more care- 
fully planned program could be devised without the need for such traps, possibly by some 
sort of compensation for these large orders of magnitude. 

The available solutions of Neumann, 16 and Braunbek-McKeehan 1 for equal currents, 
and of Cain and Williams 4 for equal radii are listed in the Table. 

EIGHT COILS. Solutions for eight coils have also been made. The Neumann- 
McKeehan 1 solution held the currents equal, and only made four coefficients vanish. Cain 
and Williams 4 were able to set seven of them zero by holding the radii equal. These 
solutions are also in Table 1. 

TOLERANCES. Raff 15 has discussed the effect of displacing one coil of a particular 
system, first parallel to the coil system axis and then perpendicular to the system axis. 

A similar solution was made for the optimum four-coil system, assuming that one of the 
inner coils is moved slightly parallel to the axis of the system. For the optimum system 


H(x) = a o + a 1 x + a 2 x2+... _ 


( 10 ) 
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where the odd a’s are zero because of the symmetry, and a 2 , a 4 and a 6 are zero because 
of the choice of dimensions and current ratios. Suppose now that one of the inner coils 
is moved a distance in the direction of its axis. Then a given a^ which is initially 
zero has a Taylor* s series, two terms of which are 




3£i 


~daf> 

AC, = ^ 


AC j = o 


*C, 


AC i = o 


AC, 


(ID 


where 


Ba^ 


1,(1 - Ui 2 )p ^ +2 (uj) a + 1 ) 


2r 


1*2 


( 12 ) 


From (11) and (12) 


-a + 1) i x (i - Ul 2 )p; (u,) 

H = r~z A Ci • 


2r 


1*2 


(13) 


Now 


= 7^- C(1 - Uj 2 ) + i (1 - u 2 2 )/r] 


where i = i 2 /ij and r = t 2 /t 1 . 
From (13) and (14), 


Using (10): 


: i H 


/a 


(^ + 1)P^ +2 (u 2 ) 


2 


1 + 


i(l - u 2 2 ) 

r(l -u x 2 ) 



H = a [l + r a/a (x/r ) + 
o 1 1 o v V 


+ t l a^/a o (x/r^ + ■••] 


(14) 


(15) 


(16) 
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If the results for the optimum four coil system are used 


- 1.69 (x/ Tj) 8 + 2.58 (x/r,) 10 + . • .J . 


The terms involving AC X , are caused by moving one coil from its optimum position, while 
the terms involving (x/rj ) 8 and (x/i-j ) 10 are the residual terms after the optimum po- 
sitions have been assumed. The latter two terms are used to obtain the optimum value 
of x/F = 21.677% given previously. The effect of the terms involving AC t is to shift the 
center of uniformity a distance D where 


D = - 


/aC,\i 

0,674 Vw 7 " 1 


AC 

1.81 I— At 


— = -0.187 r. 


(18) 


If the coil which is shifted is to the right of the center of the system, the center of 
uniformity is shifted to the left of the center of the system by an amount given by (18). 
Furthermore, if AC x is positive (the coil is shifted to the right), H increases as the dis- 
tance from the center of uniformity increases. Similarly, if the coil is shifted to the left, 
H decreases as the distance increases. The length L of the uniform region is given 
approximately by 


L = 




0.5548 

(ACj'rj)* 


(19) 


where 8 = Id' 5 for 0.001% uniformity. 

Suppose that AC X = 1/32 inch and r x = 10 feet, then ACj/rj = 2.6 X 10"*, L = 2.92 
feet and D = 1.87 feet. Hence, the effect of moving a coil 1/32 inch from its optimum 
position is to produce a region of uniformity approximately that of a sphere 2.92 feet in 
diameter, and to shift the center of this sphere approximately 1.87 feet to the left. If the 
field along the axis of an actual coil system could be measured, it might be possible to 
determine which coil must be shifted, and how much and in which direction it should be 
shifted. 
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CONCLUSIONS 

The optimum four-coil solution has been found, while other solutions obtained include 
those of various earlier investogators. The best six-coil solutions known at this time are 
also given. The effect of moving one coil out of position in a four-coil system is calcu- 
lated. This one calculation seems to indicate that the position must be held well within 
1/32 inch if 0.001% uniformity in a three-foot sphere is required. 

It is recommended that: 

1. Additional work on tolerances be done, particularly on moving the coil perpendi- 
cular to the system axis and on deforming the coil. Work on the coil-winding depth and 
breadth would also be useful. 

2. The six-coil system solution should be finished. 

3. A study should be made of rectangular coil systems, because these might be 
easier to construct than circular ones. 

4. A model of a four-coil system, say five feet in diameter, would be helpful in 
checking the conclusions, particularly those of the tolerance study. 

NOTE: Further work on the calculations of the six-coil system is now under way. 

The complete manuscript with all computations will be submitted for pub- 
lication as a Goddard Technical Note. 
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RESPONSE OF VEHICLE STRUCTURES TO ACOUSTIC FIELDS 

PRELIMINARY STATEMENT 

(Dr. J. E. Greenspon, Contract Consultant) 

During certain phases of spacecraft flight, the structure is subjected to intense jet 
noise and thrust reactions which can cause violent vibrations. The payload can receive 
vibrations through many paths: Vibrations of the engines can be transmitted directly 
through the structure to the payload attachment points; also, a certain portion of the 
acoustic load can be transmitted to the payload through the air. Part of the acoustic 
energy will excite vibrations of the casing. One aspect of Goddard's program in this 
field is an attempt to predict the response of the structural shell and the payload when 
subjected to this intense acoustic environment. Internal pressurization, temperature of 
the fuel inside, and structural damping of the casing are a few of the factors affecting the 
response of the shell. 

The response of structures to these intense acoustic fields has been studied by a 
number of investigators. In particular, the writer has developed certain expressions for 
the response of unstiffened cylindrical shells in terms of the sum of complex frequency- 
response functions for sinusoidal-loading conditions. These response functions have been 
extended to cover random loading, and are programmed for the electronic computer. 
Results of this analysis have been correlated with experiment. 

Thus far, only the response of unstiffened shells has been computed; the next logical 
step is to consider the shell with ring stiffeners. That is the purpose of the present work 
in the Summer Workshop Program. 
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EFFECTS OF RESIDUAL STRESS ON THE 
VIBRATION OF RING-STIFFENED THIN 
CYLINDRICAL SHELLS 


by 

D. F. Haskell 


ABSTRACT 

In this investigation, expressions are presented for the total 
potential and kinetic energies of one bay of a ring- stiffened, cir- 
cularly cylindrical, thin shell subjected to internal hydrostatic 
pressure. A procedure is explained for determination of the 
natural frequencies and modes of vibration, from the equations of 
motion which may be formulated by means of Lagrange's equa- 
tions. In addition, a method is presented which may be used to 
determine the effect of residual stress upon the vibration charac- 
teristics of the system. 

The equations presented in this report are general equations 
which may be employed for analysis of axial excitation of the 
structure as well as of the radial excitation considered here. An 
extension of the present theories is required before a comprehen- 
sive method can be obtained for theoretical prediction of launch 
vehicle and satellite response to the various vibration excitations. 
The response excitation is not yet available. The vibration char- 
acteristics of longitudinally stiffened shells should be studied as 
well as those of cross-stiffened shells. The response of a bulk- 
head to excitation from various shell configurations, the effect of 
residual stress on longitudinally and cross-stiffened shells, the 
response of a satellite to bulkhead- and shroud-vibration inputs, 
and the response of satellites to various vibrations of solid-fuel 
launch vehicles are some of the areas that require future study. 


INTRODUCTION 

A satellite launch vehicle system is a complicated structure. During various stages 
in its flight, this structure is subjected to severe vibrations which can cause direct phy- 
sical damage to the system. Two examples of deleterious effects of the vibration environ- 
ment are malfunction of electronic and mechanical equipment due to induced oscillatory 
motions, and generation of oscillatory stresses leading to fatigue failure (actual fracture 
or fretting fracture wear)! Such effects can damage both the launch vehicle and the satellite. 
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The vibration arises from various sources, some of which are combustion oscilla- 
tions in rocket motors, aerodynamic flutter, and noise from the rocket exhaust. The 
combustion oscillations in the rocket motor can be transmitted up through the fuel -storage 
section or casing of the launch vehicle to the satellite in the form of axial vibrations. On 
the other hand, noise from the exhaust transmitted through the surrounding air can im- 
pinge on the external surface of the launch vehicle and satellite shroud in the form of 
radial vibratory motion. 

Of central importance in the test and evaluation of satellites is the development of a 
theory to enable prediction of satellite response to this external acoustical excitation. 
Preliminary studies in this area have been conducted by A. Powell, 2 I. Dyer, 3 S. H. Cran- 
dall and A. Yildiz, 4 and J. E. Greenspon. 5 Greenspon has presented theoretical analysis 
of the random excitation of unstiffened cylindrical shells due to jet noise, and has com- 
pared the results of his theoretical analysis with experiment. The next step in the devel- 
opment of a theory to predict satellite response to external loading is to consider the 
random response of stiffened cylindrical shells. Therefore, the present investigation is 
concerned with the development of a theory to predict vibration characteristics of ring- 
stiffened, thin cylindrical shells. The energy method used in the analysis is extended to 
include the effects of residual stresses in the shell which would arise if the stiffened 
rings were welded to the shell. The latter study may be useful in design because, if the 
effects of residual stress on the characteristics of the vehicle could be predicted, then 
stiffened rings could be welded to the shell, rather than chemically milled, which is the 
present procedure. This would serve to reduce the fabrication costs of the vehicle. 


THEORETICAL ANALYSIS 

RING-STIFFENED THIN SHELL WITHOUT CONSIDERATION OF RESIDUAL 
STRESS. For purposes of this analysis, the structure of a launch vehicle is considered 
to be a multiple-bay, ring-stiffened, thin cylindrical shell which is subjected to internal* 
hydrostatic pressure. The method used is formulation of the equation of motion of one 
bay of the system, by means of Lagrange’s equations for continuous motions. 

The total potential of one bay of the system is taken as the sum of the change in ex- 
tensional bending, plus the shear-strain energies of the shell, total potential of the finite 
rigidity stiffening rings, and total work done by external forces acting on the shell. The 
sum of these expressions, equations (1)** was developed by Nash 6 and is written in terms 
of displacement components of a point in the middle surface of the shell. The kinetic 
energy of the system consists of the kinetic energy of the shell and of the stiffening rings. 
The shell kinetic energy, T s , is given by equation (2), and the kinetic energy of the stiffen- 
ing rings, T sr , is given by equation (3); these expressions are also written in terms of 


• This presentation will hold either for internal or external pressure. The equations presented in this report are for 
positive external pressure. 

* * Equations are given in the Appendices 
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the displacement of a point in the midsurface of the shell. The elasticity of the stiffening 
rings is taken into account by assumption of enough terms in the displacement functions, 
equation (4), to allow the rings to translate, dilate, bend in their plane, and deform in the 
circumferential direction. 

By substitution of the displacement configuration (4) into equations (1), (2), and (3), 
the general expressions for total potential of one bay, shell kinetic energy, and stiffening- 
ring kinetic energy -equations (7), (5), and (6) respectively-are obtained in terms of the 

system constants, generalized coordinates, A, B N, and the time functions U, V, and 

W. To obtain the equations of motion, equations (5), (6), and (7) must next be substituted 
into Lagrange's equations (12). The natural frequencies of vibration may then be calcu- 
lated from the equations of motion by assuming suitable time functions. 

RESIDUAL-STRESS CONSIDERATIONS. The previous development may be extended 
to include the effect of residual stress upon the vibration characteristics, by adding an 
equivalent residual-strain energy term to the total potential energy expression. To do 
this, it is assumed that an "equivalent residual-stress load" may be determined. By this 
procedure may be calculated the external hydrostatic pressure required to produce the 
same deflection as that which would be caused when stiffeners are welded to a cylindrical 
shell. 

It has been found 7 that a cylindrical shell acquires an initial measurable deflection 
between stiffening rings where the rings are welded to the shell; if a static external- 
pressure distribution can be determined which produces the same deflected shape as 
would the actual residual-stress distribution, then the change in potential due to the 
actual residual-stress distribution may be replaced by the change in potential due to this 
"equivalent residual-stress load." 

In Appendix B, the static external-pressure distribution and initial displacement 
configuration are assumed to have the same form as the second-degree curve measured 
by Lunchick and Short. 8 These distributions are substituted into Flugge's equations of 
motion, 9 and the "equivalent residual-stress load" is determined in terms of the shell 
parameters and assumed maximum initial shell deflection-this maximum initial shell 
deflection may be either assumed or taken from experimental evidence. This calculated 
"equivalent residual-stress load" (po) is then added to the applied hydrostatic stress (p) 
as a sort of correction factor for residual stress. In this manner, the effect of residual 
stress is introduced to the change in total potential of the system through the pressure, P. 
After this is accomplished, the equations for kinetic and total potential energies presented 
in Appendix A (which excludes the effect of residual stress) may be extended to include 
the effects of residual stress, by including in the total pressure term, P (which only 
included the applied hydrostatic pressure, p, in Appendix A), the "equivalent residual- 
stress load," po ; i.e., instead of P = p, P is set equal to the sum of po and p for the case 
where residual stress is considered. 
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APPENDIX A 

KINETIC AND TOTAL POTENTIAL ENERGY EXPRESSIONS 


The total potential of one bay of the system, which includes elastic extension, bending, 
and shear strain energies, work done by the external load and the ring strain energy, may 
takes the form 


2?rR L r 

n n -I 


+ K, 




* 2 - - -> 



+ 



2 Bfo. 

R 2 3 S 2 


Cl) 


dxds 


2ttR l 

+k 3 J J ® + S ) 2 dxds 

o J o 

2itR L _ 

*™j / [- 1 (if ) a * i (sr - i - 1 (H) 2 ] — 

0 0 


bh E 

o 

+ 2 


277R1 r , _ . 

f f 3y M \ 2 2 3v 

J |U» R/ + (rj "RB-s-J x=0 


ds 



( 1 ) 


6-9 


SUMMER WORKSHOP REPORT 


The kinetic energy of the shell, T s , and of the stiffening rings, T SR , may be given by 


T s = 


2itR L 

f f (u 2 + v 2 + <£> 2 ) dx ds 

Jo Jo 


( 2 ) 


and 


t sr “ 


b(h n -t)p 


2.R2 

I (u 2 + V 2 + cl) 

Jo 


ds - 


(3) 


The displacements during vibration to be used in these kinetic energy and total po- 
tential expressions are taken as 


u = j^A s in -jp sin28x+B s in ^ cosSx+Ccos ^ sin 26x + D cos ^cos 6x + F -j-J U 


/ v = ^sin— - (G + H sin 28x + J cos 2Sx)J V 


(4) 


= sin + Mcos^j (1 - cos 2Sx) +nJw. 


After the displacement configuration (4) is substituted into equations (2), (3), and (1), 
the shell kinetic energy, ring kinetic energy, and total potential of one bay of the system 
take the following forms: 


t s = k 4 (k 10 u 2 +k u v 2 +k 12 w 2 ), 
t sr =K s (k 13 U 2 + k 14 V 2 + k 15 W 2 ), 

bf = [(k 38 + kjy + kjg)U 2 + (k 39 + k2Q + k2^)V 2 + (k 34 + Pk 33 )W j 
+ t( k 36 + k 37 ) UV + (k 3g - k 3g )VW + (k 40 + Pk 41 )WU] , 

where 


(5) 

( 6 ) 

( 7 ) 


K 4 = wRLpt K s yobR(h o - t) 

k 34 = ( k 22 + ^23 + ^24 + ^25 + ^26 + *^27 “ ^28 + ^32^ 
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k 3S - (-k 29 + k 30 - k 31 + k 33 ) 
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APPENDIX B 

DETERMINATION OF THE EQUIVALENT RESIDUAL-STRESS LOAD 


The displacement configuration caused by residual stress and external pressure dis- 
tribution is assumed to be of the form 



and are to be substituted into Flugge's equations of motion, (5) which are 
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. S = — £l_ . 

(i -v 2 ) 

If the maximum deflection caused by residual stress is written in terms of a fraction, n, 
of the shell radius 


A = 7ja 


and it is assumed that 


Px = Ptf> = o 


Pr =Pi, 

then, after substitution of (8) and (9) into Flugge’s equations (10), it is found that the 
"equivalent residual-stress load" is given by 


Po = 


Et(12a 2 + t 2 )' 
,12a 3 (l-v 2 ) y' 


(ID 


As indicated, this "equivalent residual-stress load" is a direct function of the magnitude 
of the residual stress, which is represented by the factor n. 
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APPENDIX C 

EQUATIONS OF MOTION FOR RING-STIFFENED THIN SHELLS 


The equations of motion for the ring-stiffened thin shell, with or without the consid- 
eration of residual stress, are formulated by means of Lagrange's equation of motion: 


_d ( 3T \ _9T_ _BU_ 

dt \Bq J + Bq t 


Q it (i = 1,2,3. . . ), 


( 12 ) 


where q, are the generalized coordinates andQj are the generalized forces. The sum 
of the kinetic energy terms, equations (5) and (6), and the total potential energy, 

(7), are substituted into (12), and the number of equations obtained is equal to the number 
of generalized coordinates used in the kinetic-energy and potential -energy expressions. 
To determine the natural frequencies of motion, is set equal to zero. 

For those cases in which the effect of residual stress is to be considered, the ex- 
ternal pressure, p, is set equal to the sum of the "equivalent residual-stress load" and 
the actual applied hydrostatic pressure: 

P = P 0 + P. 

For those cases in which residual stress is not a factor, the external pressure term is 
only equal to the actual applied pressure: 


P = P . 
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RELATING DAMAGE AND RELIABILITY TO ENVIRONMENTAL DATA 

PRELIMINARY STATEMENT 

(R. E. Dorrell, Quality Assurance Branch, Test and Evaluation Division, GSFC) 

The environmental effects on systems elements, be they minimal or extreme, may 
cause severe damage, resulting in catastrophic failure to the entire spacecraft system. 
This is the type of "hard-nosed" problem which Goddard is facing and which Goddard 
must solve. The Quality Assurance Branch is responsible for making sure that, for 
every ground and spacecraft application, the proper material is used and that this mate- 
rial is used properly. Two recent occurrences will illustrate some of the problems that 
the Quality Assurance Branch encounters. 

The first will exemplify the case of minimal environment and maximum damage. A 
failure was reported involving the IN 461 and IN 252 diodes made by the Sylvania Elec- 
tric Company. The first indication of failure was a complaint that several of the diodes 
broke apart as their leads were straightened out upon being taken from the card in which 
they were shipped. Also, when they were put on printed circuit boards or otherwise 
soldered in place, a good many open diodes were encountered. 

In order to determine just what was causing these failures, the opaque coating cov- 
ering the diode to prevent any penetration of light photons was removed with a chemical. 
The first thing noticed was that, as soon as the leads were straightened out, failure of 
the glass resulted. This was no doubt due to a lack of annealing of the glass. During 
manufacture, although the semiconductor wafer itself must be kept relatively cool, the 
glass should be hot in order to provide a good seal. If the glass were not annealed, failure 
due to cracking of the glass would and evidently did result. Next, as soon as the diodes 
were put in circuits and cutting pliers were used to trim off the leads, open diodes re- 
sulted. Apparently the slightest of mechanical shocks, such as the action of cutting the 
lead wires, was sufficient to open-circuit the diode. Finally, the diode was examined 
under a microscope. To briefly describe the construction, an aluminum rod is pressed 
against the silicon wafer, heat is applied, the aluminum rod is cut off, and a pressure 
contact made with an S-shaped ribbon. However, under a microscope, the contact ap- 
peared to be almost completely severed, a factor which again might easily have led to 
failure. 
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The ideal method of relating all this to environmental parameters would be a vibra- 
tion or shock test to determine whether or not failures of this type are present. However, 
when these poorly constructed diodes were placed in competition with others superior in 
design, on the same vibration test, they showed up every bit as well. In this vibration 
test, the forward characteristic was measured first; then the vibration test was conducted ; 
then the forward characteristic was measured again-the two coincided very well. How- 
ever, if this diode were used in a compter, any transient that the diode generated would 
give false indication or false signals. To perform a valid test, therefore, a current is 
passed through the diode and the forward current characteristic is displayed on an oscil- 
lograph; a vibration test is then conducted, and any discontinuities in this curve indicate 
a poor diode. 

This was a case involving minimal environment because the diodes were not mis- 
treated at all. Now, we shall discuss a problem in which the environment is more ex- 
treme. It is this problem toward which the Summer Workshop people will be contributing. 

While the UK-1 satellite was being tested, particularly in thermal vacuum, a prob- 
lem arose concerning a transistor used in a regulator circuit to govern the amount of 
energy put into the battery from the solar paddles. The regulator circuit consisted of a 
detector and power transistor, the transistor being used solely to dissipate excess power 
from the solar paddles. The system worked well in room temperature, but in thermal 
vacuum the transistor junction and case temperatures rose tremendously and the transis- 
tor eventually failed due to "thermal runaway." The transistor was rated for 50 watts 
while on a heat sink, yet it was only handling a maximum of 5 watts for this application, 
so the question arose of why this thermal runaway occurred. 

The transistor without a heat sink was rated for only 3 watts, so one of the first 
things to be looked into was the matter of heat sinking. A power transistor generally is 
constructed by mounting the wafer on top of a stud or to the header that forms the base 
or lower end of the case. Now, the energy to be dissipated in the form of heat must be 
conducted away from the stud at the bottom to a metallic sheet or heat sink. However, it 
is very inconvenient to electrically connect the transistor collector power supply to the 
spacecraft ground or the frame. Therefore, an electrical insulator must be interposed 
between the transistor base and the heat sink, but must nevertheless provide a good heat- 
conductive path between transistor and heat sink. The interposed material takes the form 
of a washer, made of beryllium oxide or aluminum oxide or some other material which 
has high electrical resistance combined with high thermal conductivity. 

Now, while the regulator was working in air, apparently there was enough convective 
cooling to alleviate the difficulty. However, since no convection is possible in a vacuum, 
the system must depend on cooling by conduction to the heat sink. The difficulty here, 
and probably the source of the problem, is the number of interfaces through which the 
heat must go to be conducted away from the transistor: first, between the wafer and the 
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stud; then between stud and can; another between can and washer; and another between 
washer and heat sink. Each interface contributes a certain resistance to heat flow, so 
that the sum of all the interface resistances may be sufficient to reduce heat dissipation 
to a dangerously low level. 

One possible solution, tried by some transistor manufacturers, is to electrically in- 
sulate the transistor with a beryllium oxide or aluminum oxide tab placed inside the 
transistor case. This would help the situation by removing the necessity for so many 
external interfaces. However, this unit must undergo vibration exposure and thermal 
analysis before its effectiveness can be determined. 

In summary, I think we have pointed out the area where we hope that the Summer 
Workshop will be able to give us some assistance. We are looking forward to seeing 
what can be evolved. 
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HEAT-SINKING TECHNIQUES FOR POWER TRANSISTORS 
IN A SPACE ENVIROMENT 

by 

J. John and J. Hilliard 


ABSTRACT 

An investigation was made of the cooling of power transistors 
in a space environment, where the only available mode of heat 
transfer is that of conduction to a heat sink and radiation from the 
heat sink to space. An attempt was made to minimize the thermal 
resistance between transistor case and heat sink, allowing the 
transistor to dissipate as much power as possible while maintain- 
ing its temperature within the maximum tolerable level to prevent 
thermal runaway. Further, it was necessary to electrically insu- 
late the transistor from the heat sink. The use of beryllium oxide 
washers provided electrical insulation, while adding very little to 
the thermal resistance between case and sink, the Be O being a 
good heat conductor. However, the problem of contact thermal 
resistance at each interface arose, especially in vacuum-this 
contact resistance providing practically all the thermal resist- 
ance between case and sink. The effect on the contact resistance 
of surface pressure, insertion of foils, and soldering was ex- 
amined. It was concluded that, for most efficient heat sinking, 
indium foil should be inserted at each interface, the indium foil 
having the effect of reducing the contact resistance in vacuum 
by a factor of 10. 


INTRODUCTION 

The problem of maintaining the junction temperature of a transistor within certain 
limits to prevent runaway and failure is well-known. The work performed to date has 
been concerned with transistor cooling in an atmospheric environment, either by free or 
forced convection. Typical of the work in this area is that performed by Abel 1 and 
others. 2 The investigation performed here, however, had the purpose of extending pres- 
ent knowledge to cover transistor cooling in a space environment. 

The only mode of heat transfer available in space is radiation, either directly from 
the transistor or else from a heat sink to which the transistor is thermally coupled. Cal- 
culations show that the surface area of the transistor is too small to provide more than 
a few milliwatts of direct radiative heat dissipation, whereas several watts may have to 
be dissipated. So, the problem resolves itself into an investigation of the thermal path 
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between the transistor junction and the heat sink with the goal of making the thermal con- 
ductance as high as possible. 

Figures 7-1 and 7-2 show how a typical power transistor is mounted on a heat sink. 
Two thermal paths are available: one from the case, through the top washer to the heat 
sink; and another from the case, through the stud and nut, through the bottom washer and 
to the sink. Essentially, then, we have the case of thermal resistances in parallel. 



Figure 7-1 - Mounting of 2N1724 
Transistor on anodized aluminum 
heat sink, using two BeO washers 



Figure 7-2 - Mounting of a typical pawn transistor, diagram 


It is desired that the net thermal resistance be as low as possible. However, it is 
often necessary that the transistor also be electrically insulated from the sink. This 
means that some sort of device such as a washer, must be used that will provide good 
thermal conductivity while serving as an electrical insulator. Further, if washers are 
to be used, they must have good mechanical properties to resist cracking when the nut 
on the transistor is tightened during mounting. Materials satisfying these specifications 
to a greater or lesser degree are beryllium oxide, 3 mica, and aluminum oxide. 
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Figure 7-3 - Thermal network 
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Figure 7-3 is a schematic diagram of the thermal network. It can be seen that all 
quantities are readily determinable with the exception of contact resistance (the thermal 
resistance at the interface of two different materials). 


Contact resistance is due to improper thermal contact between the two materials. 
Since no surfaces are perfectly smooth, they only touch at a limited number of points; the 
remainder of the space being filled with a nonconducting vacuum. Variables affecting 
the contact resistance would include the pressure between the two surfaces, the smooth- 
ness of the surfaces, the materials themselves, and the possible addition of greases, sol- 
ders, or some soft material to fill the empty regions between the surfaces. 

Several theoretical and experimental investigations have been made in the general 
area of thermal contact resistance. Theoretical considerations must necessarily assume 
an idealized shape of contact point and a distribution of contact points. An analysis car- 
ried out by Fenech and Rohsenow 4 attempted to predict, with some degree of success, the 
thermal conductance of metallic surfaces in contact. It is felt, however, that an analysis 
of this type can only be at best a fair approximation, and that experimental values should 
be determined. Several investigators have worked on the experimental side of the prob- 
lem. 5 ’ 6 - 7 Their results are not considered applicable to this problem, because their data 
were obtained in air, not in vacuum; also, their data necessarily depended on the configura- 
tions of the surfaces they employed, somewhat different from the ones to be dealt with here. 

Tests were run to determine the thermal resistance from transistor case to heat 
sink as a function of the previously mentioned variables. It was hoped to minimize this 
thermal resistance so as to provide a practical means of heat sinking in a space 
environment. 


EXPERIMENTAL APPARATUS AND PROCEDURE 

In this experiment, the junction temperature, case temperature, sink temperature, 
and energy dissipated per unit time were measured for each test. Variations in contact 
resistance were introduced by applying different torques on the transistor nut with a 
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torque wrench, and by employing various interface materials such as foils and solders. 
Variation in the smoothness of the surfaces was not undertaken because of the difficulty 
of measuring and controlling this parameter; in any case, it was felt that this could offer 
only limited improvement. 

Case and sink temperatures were measured with copper-constantan thermocouples 
and a potentiometer. The thermocouples were fastened to the surfaces with aluminum 
tape, as shown in Figure 7-4. Two thermocouples were used on the heat sink, one at the 
base of the transistor and the other at the extreme edge of the heat sink. 



Figure 7-4 - Location of Thermocouples 
on anodized aluminum heat sink 

As the collector junction was inaccessible to a thermocouple, its temperature had to 
be measured indirectly. Junction temperature (t . ) was determined by measuring the for- 
ward voltage drop (V CBC ) from collector to base with emitter open, since this voltage is 
directly dependent on junction temperature. The appendix contains an explanation of the 
technique employed and a schematic of the circuit used to make this measurement. 

Figure 7-1 shows a typical power transistor (2N1724, which was used for all tests) 
mounted on a 10 x 10 x 1/16-inch black anodized aluminum plate which served as a heat 
sink. In this configuration, beryllia washers were placed between the transistor case and 
the plate, and also between the nut and the bottom of the plate. The stud of the transistor 
was 1/4 in. in diameter and the hole through the plate was 5/16 in. The base of the tran- 
sistor had a diameter of 3/4 in. All the beryllia washers employed were of the same 
thickness, 1/16 in. 
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An environmental pressure of approximately 10 5 mm Hg was obtained by the use of 
a bell-jar vacuum system employing a mechanical fore pump and an oil diffusion pump. 

In the tests made under both atmospheric and vacuum conditions, the measuring proce- 
dure was the same. Power was applied to the transistor from a constant-voltage power 
supply and all temperatures were allowed to stabilize for approximately 30 minutes. The 
junction, case, and heat-sink temperatures were then measured, and the temperature dif- 
ference between the case and sink was calculated. The thermal resistance was obtained 
by dividing the temperature difference by the input power. The input power was set at 
several different levels and the procedure was repeated at each level. 


RESULTS 

Table 7-1 is a compilation of the experimental data obtained in this investigation, 
both in air and in vacuum. 


TABLE 7-1 

Experimental Data Obtained in Air and in Vacuum 


Configuration 

Vac. or Atm. 

Torque 
(in. lb.) 

(° C/ watt) 

2 BeO washers 

Vac. 

6 

4.79 

2 BeO washers 

Air 

6 

1.19 

2 BeO washers + 3 indium foil 

Vac. 

6 

.61 

2 BeO washers + 3 indium foil 

Air 

6 

.50 

2 BeO washers + 3 aluminum foil 

Vac. 

6 

1.59 

2 BeO washers + 3 aluminum foil 

Air 

6 

.81 

2 BeO washers + 3 copper foil 

Vac. 

6 

3.15 

2 BeO washers + 3 copper foil 

Air 

6 

1.04 

2 BeO metallized washers + 3 indium foil 

Vac. 

3 

.71 

2 BeO metallized washers + 3 indium foil 

Air 

3 

.54 

2 BeO metallized washers + 3 indium foil 

Vac. 

6 

.66 

2 BeO metallized washers + 3 indium foil 

Air 

6 

.50 

2 BeO metallized washers + 3 indium foil 

Vac. 

9 

.61 

2 BeO metallized washers + 3 indium foil 

Air 

9 

.43 

2 BeO metallized washers + 3 indium foil 

Vac. 

12 

.55 

2 BeO metallized washers + 3 indium foil 

Air 

12 

.42 

2 BeO metallized washers 

Vac. 

6 

4.40 

2 BeO metallized washers 

Air 

6 

1.12 
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Table 7-1 (Continued) 


Configuration 

Vac. or Atm. 



2 BeO metallized washers + 2 indium foil + 
transistor soldered to washer (indium solder) 

Vac. 

6 

.49 

2 BeO metallized washers + 2 indium foil + 
transistor soldered to washer (indium solder) 

Air 

6 

.40 

Transistor soldered to metallized washer 
(indium solder), washer soldered to plate; 1 
BeO washer + indium foil on bottom 

Vac. 

6 

.56 

Transistor soldered to metallized washer 
(indium solder), washer soldered to plate; 
1 BeO washer + indium foil on bottom 

Air 

6 

.43 

1 BeO washer + 2 indium foil; stud insulated 

Vac. 

6 

.70 

1 BeO washer + 2 indium foil; stud insulated 

Air 

6 

.49 

2 BeO washers + 2 indium foil + Apiezon 
grease on threads 

Vac. 

6 

.57 

2 BeO washers + 2 indium foil + Apiezon 
grease on threads 

Air 

6 

.47 

No washers 

Vac. 

6 

5.16 

No washers 

Air 

6 

1.36 


•Washer dimensions - O.D. - .90" 
I. D. - .26" 


Early in the program it became apparent that the thermal resistance of the beryllium 
oxide washer ( R = L /Ak = .02°C/watt) was much lower than the thermal contact resis- 
tances. Hence, instead of testing the effects of different washer materials, it was decided 
to use beryllium oxide exclusively, since little improvement over it could be expected. 

EFFECT OF ENVIRONMENT. Figure 7-5 illustrates the increase in thermal re- 
sistance between the transistor case and sink due to the vacuum environment, as com- 
pared to the resistance in air. In a vacuum, the only mode of heat flow between the sur- 
faces is conduction through the few discrete contact points; in air, there is also the 
possibility of heat convection between the surfaces or heat conduction across the narrow 
air layer. The presence of a vacuum thus tends to amplify the problem of heat sinking. 

USE OF FOILS. Figure 7-6 illustrates the effect of using foils along the interfaces 
of the beryllia washers in a vacuum environment. In each case, a foil washer with the 
same surface dimensions as the beryllia washer was placed between the transistor case 
and top washer, between the top washer and the heat sink, and between the bottom wash- 
er and heat sink. The graph shows that the use of any of the three interface materials 
lowered the thermal resistance from case to heat sink. Indium foil, however, proved most 
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effective, lowering the thermal resistance to .63 ° C/watt, approximately 1/8 of its value 
without interface material. 

The effect of the foil is to fill the void between the surfaces with heat-conducting 
material. The softness of the foil seems very critical in determining its effectiveness. 
Indium, by far the softest of the foils, in most cases was found to adhere to the contact 
surfaces due to the penetration of contact points into the foil. 

SURFACE PRESSURE. The effect of surface pressure on thermal resistance is 
shown in Figure 7-7. The transistor was mounted with two beryllia washers and three 
indium washers as interface material, and torque on the nut was varied from 3 to 12 
in. lb. in 3 in. lb. steps. As the graph illustrates, thermal resistance decreases linearly 
with an increase in surface pressure. Greater pressures probably would result in a 
still smaller thermal resistance; however, the cracking of the washer limits further 
tightening of the nut. 



Figure 7-7 - Effect of surface pressure 
on thermal contact resistance in a vac- 
uum environment 


This increase in conductivity is thought to be due to further penetration of the indium 
into the voids as the surfaces are further compressed. 

EFFECT OF SOLDERING ON CONTACT RESISTANCE. One way of reducing ther- 
mal contact resistance is to solder the surfaces together, in this case, soldering tran- 
sistor to washer and then washer to heat sink. As the ceramic washer material cannot 
be directly soldered, the beryllium oxide must be metallized and then copper- or silver- 
plated. The metallizing is accomplished by depositing on the washer surface a 
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molybdenum manganese compound, a process carried out by the National Beryllia Cor- 
poration or the Brush Beryllia Company. The first tests made were for the purpose of 
investigating any change in thermal contact resistance due to the plated washers over 
those not plated. Using the same configurations, plated and unplated washers gave the 
same value for contact resistance. Next, the plated washer was soldered to the transis- 
tor and also soldered to the aluminum heat sink, the heat sink being nickel -plated to hold 
the solder more easily. A low-melting-point indium-alloy solder was used on the tran- 
sistor to avoid any possible damage to it. Results showed that the soldering reduced the 
thermal resistance to roughly the same level as that achieved using indium foil in the 
interfaces. 

To determine the magnitude of heat conduction through the stud, nut, and bottom 
washer in relation to that path through the top washer, the stud was insulated with teflon 
and paper, and conduction was allowed only through the upper path. The change in ther- 
mal resistance was so small that the lower path might almost be neglected when using a 
good conductive configuration such as the indium-beryllia washer method for the upper 
path. The high resistance of the lower path is thought to be caused mostly by ineffective 
contact between the screw threads and the nut, and between the nut and beryllia washer. 
An attempt was made to improve thermal conduction through the lower path by applying 
a film of Apiezon vacuum grease along the threads and putting indium between the nut and 
the beryllia washer; however, this made no appreciable change in overall thermal 
resistance. 

Several power transistor manufacturers have attempted to combine the necessary 
thermal conduction and electrical insulation by insulating the collector from the transis- 
tor case internally, as in the 2N1724/I (which is identical to the 2N1724 except for this 
modification): this transistor therefore may be mounted directly on the heat sink. A test 
run in vacuum showed that direct contact between the two metal surfaces still gave a 
high thermal resistance, 5.16 °C/watt. This could be lowered considerably by the addi- 
tion of indium foil to the interface. However, measurements made with the 2N1724/I 
transistor showed it to have a thermal resistance from collector junction to case of 
1.17 ° C/watt, while that of the 2N1724 transistor was .36° C/watt. If the latter transis- 
tor is provided with the indium-beryllia washer method of heat sinking, the total resis- 
tance from junction to heat sink is .98° C/watt, which is less than the junction-to-case 
resistance alone of the modified 2N1724/I transistor. 

CONCLUSIONS 

In summarizing, several facts stand out as important in the problem of heat-sinking 
a power transistor in a space environment. The use of a soft interface foil is highly 
effective. Indium foil, used as an interface material, reduces interface resistance to 
almost 1/10 of its normal value, and the ease with which it may be shaped to the required 
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geometry makes it very desirable. Aluminum and copper foils which offer some reduc- 
tion in interface resistance, could be used depending on the amount of thermal conduc- 
tivity desired. 

The effect of surface pressure, although critical when no interface material is used, 
is not of too much importance when indium is employed. At very high pressures, the 
effect could probably be made appreciable; however, the cracking of the beryllia washers 
prevents this. It might be noted that the use of indium along the faces of the washers 
allows a sizeable increase in the amount of torque that may be applied to the transistor 
nut before the washers crack. 

The process of soldering the transistor to the BeO washer and the washer to the 
heat sink provided effective heat sinking. The indium eutectic solder used between the 
case and washer proved very satisfactory, not only because of its conductive properties, 
but also because of its low melting point, which decreases the danger of harming the 
transistor. 

RECOMMENDATIONS 

As a result of this investigation, recommendations can be made concerning the most 
efficient methods of heat-sinking a power transistor in a vacuum. First, the insertion of 
indium foil at the interface between beryllium oxide washer and heat sink and between 
the beryllia washer and transistor base is sufficient to reduce the thermal resistance 
between case and heat sink to a very low level, and thus to aid in maintaining the junction 
temperature within its maximum permissible value, while the transistor is able to dis- 
sipate a fairly large amount of power. The use of indium in space will necessarily de- 
pend on its rate of sublimation and consequent deterioration. One reference, 8 however, 
indicates that, even at a temperature of 400°C, the rate of sublimation of indium in 
space is only 10 - 5 centimeters per year, which is less that that of lead or zinc, for ex- 
ample. Thus, since the temperature encountered by the indium in this application is less 
than 100°C, it is felt that sublimation will be no problem. 

A second recommended procedure for heat sinking in vacuum is to metallize the 
BeO washer and to solder the transistor to washer and the washer to the case. This 
method provided roughly the same thermal resistance from case to sink as the previous 
method involving indium foil. 

The first procedure recommended appears more advantageous, mainly because of its 
ease in assembly and adaptability. Once the joints have been soldered, it would be im- 
possible to easily remove the transistor from the heat sink. Further, it is felt that the 
former method gives a more reliable joint, which would better tolerate vibration and 
shock. 
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APPENDIX 

MEASUREMENT OF THE JUNCTION TEMPERATURES 

The transistor junction temperature, j can be determined by measuring electrical 
parameters of the transistors which are functions of T j# The parameter used in this 
experiment was the forward voltage drop^^, from collector to base with emitter open. 
A schematic of the circuit to make this measurement is included as Figure 7-8 (also see 
Reference 9). 



Figure 7-8 - Circuit for measuring heat dissipation 
and Vqbq of a transistor, schematic diagram 

To determine the relationship betweenV CBO and T for each transistor, the transistor, 
in a nonoperating condition, was placed in an oven and heated to temperatures in the 
transistor’s operating region (between room temperature and 100°C). At each tempera- 
ture, the transistor was allowed to reach a constant temperature throughout its structure. 
The switch (S-l of Figure 7-8) was then closed, causing the relay, D-l, to energize mo- 
mentarily at set intervals, permitting v CBO to be measured on the oscilloscope and the 
calibration curve to be plotted (Figure 7-9). 

In the actual tests, the transistor was mounted on a metallic plate serving as a heat 
sink, and placed in a bell jar, which was then evacuated to a vacuum of 10 * 5 mm Hg (Fig- 
ure 7-10). Electrical connection was provided between the transistor and the circuit of 
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Figure 7-9 - Junction temperature vs , 

2N 1724 transistor 

Figure 7-3. Power into the transistor was varied by means of the rheostat in the base 
lead of the transistor; the emitter-to-collector voltage, V EC , was kept constant at 14 volts 
while the emitter current, I E , was varied. At each power level, the transistor tempera- 
ture was allowed to stabilize, assuming that the input power would be equal to the heat 
dissipated per unit time. For each configuration of heat sinking, three or four different 
power levels were used. Switch S-l was then closed causing the multivibrator to be acti- 
vated; subsequently, power was switched from the power circuit to the V CB0 measuring 
circuit and back. The transistor remained in an operating condition for intervals of 2.25 
seconds and was in the measuring circuit for 100 milliseconds. Since the switching time 
was in the order of 2 milliseconds, it was possible to read v CB0 before the junction tem- 
perature had decreased from its operating value. v CB0 was then measured, converted to 
degrees centigrade, and recorded. 
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Figure 7-10 - Experimental equipment for measurement 
of junction temperatures 
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